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■ Pages in the original (foolscap paper) which had text longer than A4 were cut and 
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0.1 bottom margin on A4 paper. The images were ranged top left against those 
margins. It was necessary to fractionally reduce the size of each image to be 
slightly less than 11.38 inches high, rather than allow automatic fitting by Word. 

■ The document was saved and then output to an Apple LaserWriter II NTX but 
output to file , not actually printed. Requests to fix margins were not over-ruled. 
This created a PostScript file of the document, about 250 MB long. 

■ The PostScript file was then input to Frank Siegert's PStill program which 
converts to PDF to yield this document. 
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1. Uses of Library Programmes 

When writing and developing programmes it is sometimes desirable to know rather more 
about some library subroutine than is available in the specification. Occasionally it 
may be necessary to modify or re-write a subroutine. It is for these purposes that 
annotated programmes are provided; they shjuld be regarded as supplementary to the 
specifications. 

Binary tapes of subroutines, punched in a form suitable for use with Assembly, are 
available from the London Computer Centre. Users are strongly recommended to take these 
tested tapes rather than to punch for themselves the -outines they require. 

A letter P is printed in the Library Index (starting with Issue 6) to distinguish 
subroutines whose programmes have been issued. 

These documents are the copyright of Ferranti Ltd. , and may not be reproduced in 
whole or in part without permission. 



2. Layout 

Each programme is preceded, where necessary, by a description of the rrocess used. 
This description is omitted if the process is adequately described in the specification 
or if it is clear from the nature of the routine. 

Where the process requires a large amount of organisation, as in some print routines, 
the programme is preceded by a flow diagram. Occasionally a comprehensive flow diagram is 
provided and this, together with a print out of the programme tape, takes the place of 
the annotated programme sheets. 

There are usually two blocks of programme on each page and descriptive titles are 
printed in heavy type where appropriate. Programmes are printed on only one side of 
each sheet so that pages may be taken out and laid side by side. Each subroutine is 
precec'ed by its cue list and various other Assembly tags. The blocks of the routine 
itself are often arranged in the order in which they are obeyed, not the numerical order 
required on the tape. 



3. Notation 

The library programmes are intended to be read in conjunction with the corresponding 
specifications; whenever possible the notation introduced in the specification is 
used in the programme. 

Most of the conventions used in annotating library programmes are introduced in the 
Pegasus Programming Manual. Lists of the more important symbols and abbreviations are 
given below: 



B 2+ Block 2+ in the Main Store 

B.P Block and Position Address 

C(2) The content of Accumulator 2 

Er Tape Character Erase 

I.O. The Initial Orders 

l,s. Least significant (or less significant) 
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m. s. 
OVR 

P 

P 
P.P. 

PQ 

Q 
q 

r. o. c. 

Sp 
sp 

U 2.1 

W 2 

X 2 

a 



7 
</> 

o 

3 C or 3 C 
3 or 3 



M 



B 



{4.3, 20) 



Most significant (or more significant) 

The Overflow Indicator 

Accumulator 6 (P for Product) 

The content of Accumulator 6 

Preset Parameter 

The double-length Accumulator, 6 and 7 

The double-length number in PQ . 

Accumulator 7 (Q for Quotient) 

The content of Accumulator 7 

Round-off constant 

Tape Character Space 

Register i.l in the Computing Store. The U is often omitted. 

Block 2 in the magnetic tape buffer 

Accumulator 2. The X is usually omitted. 

The content of Accumulator 2 

The a-search is the section (or setting) of an input routine which 
ignores blank tape 

The /3-search is the section (or setting) of an input routine which 
searches for the start of a number or a letter shift 

The 7-search is the section (or setting) of an input routine which 
searches for a warning character or directive. 

Tape Character Figure Shift (= Blank Tape) 

Tape Character Letter Shift 

Tape Character Full Stop (= Decimal Point) 

The counter in Accumulator 3 

The modifier in Accumulator 3 

The block part of the modifier 

The position part of the modifier 

A modifier of 4.3 and a counter of 20. 
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4. Symbols on Programme Sheets 







R 3-0 1 
200 - 01 - 


















0+ 




f ftl \ 




/l. 5 4 00 \ 
I 1.0 6 10 J 




< 01 ) 


1.0 

1 

c 

2 






i+H 


72 

71 






1 2+ 7 




from 2+.6 


1.2 6 63 
1.4 1 00 






3 
« 

4 

S 

» 

6 






+0 






[5] 72 5 
1.6 60 






20 - 50 0. 
6 




<C\ 


1.7 5 00 
6 3 00 




VH> 


f > 

7 






| 2. 6 3 21 5 

2.6 6 10 5 
















» 






2.0 
1 
2 


1.7 5 66 

(j) 3 51 

(~ty 74 

(2.3) 1 40 




,to 2+.3 


[3]73 5 
24.r2l TO 








**[$ 







Main Store Block Number 

Oval indicates cue 01 

Brackets show that the order pair is overwritten 

Box round block number 

Box round single word address 

Arrow marks a loop stop 

Arrow indicates a jump to this order 



Solid line under the cue 
+ cue indicating an unconditional jump 



Vertical line marks a LINK 

Solid line under an unconditional jump 

Vertical line marks a pseudo order-pair 

Arrow indicates return from subroutine 
Letter in circle refers to flow chart 

Vertical line (may be a broken line) 
marks an order which is used as a number 



Line and arrow mark a jump 

Rings indicate that the N part of the order 
is a number, not an address 

Ring indicates that the address 2,3 is 
set in l c , not C(2.3) 



Solid line under an unconditional transfer of 
control to B2+.3 



Computer Department 
West Gorton, 
Manchester, 12. 



Perranti Ltd 



London Computer Centre, 
21, Portland Place, 
London, W. 1. 



April 1958 



Copyright Reserved 
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GENERAL PURPOSE NUMBER PRINT 



A self-preserving subroutine for printing scaled numbers in various styles. The 
style is controlled by two parameters. The first of these, known as the Indicator I, 
controls the digit layout; the second, k, is used as a scaling factor. Both are 
described in the Specification, but further details of I are given below to explain 
the operation of the programme. 



1. The Layout-Parameter (see Specification, Section 6) 

The binary digits of the layout-parameter are allocated as follows 

If J-< X. J-« 5 J (p 9 q 9 ) 



I = 



I 



*8 


*1 ° *3 U U U <7 


X 


S 


^00 *oi ^02 9oo 'oi 9 08 







ho hi hi ?io ?n lia 


P 20 P 21 P 22 <7 20 <7 al <7 22 


^30 hi ^32 ?30 ?31 ^32 



(P, q t ) 



(P* <?,) 



<P a t a ) 



2. Group Counters 

For the purpose of printing, the decimal digits of each number are treated in 
four groups. Each of these has a 3 bit counter pj and a 3 bit indicator q j. 

The four groups are taken in order. If pj ^ the routine will print (or 
suppress or omit) the next pj <1 ^ Pj ^ 7) decimal digits, as specified by various 
indicators. 

pj ~ indicates the end of the number and, in this case only, the octal number 
given by the group indicator q j gives the number of spaces which are to be punched 
after the number. When these final spaces have been punched, the LINK is obeyed. 



3. Group Indicators 

q j = Digits not printed will be suppressed, (i.e. replaced by Spaces) 

qj - 1 Digits not printed will be omitted altogether. 

qj 1 and q y 2 are used to indicate the action to be taken at the end of group j. 

<? 71 qj 2 = 00 Initiate 'special operation', see Section 6, then proceed to next 
group. 

qj^ qj 2 - 01 Punch two spaces, then obey LINK. 

Iji Iji = *° Punch one space, then look at next group indicator, 

q j x q j 2 - 11 Look at next group indicator immediately. 



4. Initial Printing 



*o = ° 
*o = 1 



= 
= 1 



gives CRLP before the number 
has no effect 

sign is punched immediately after the initial space. 

has no effect. This is used for unsigned numbers, or for numbers 

printed with delayed sign. 
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Immediately following the CRLF (if any) 7-fi spaces are punched before anything 
else. 

All numbers are punched ultimately as integers in the range < z < 2 38 . If 
m is the maximum number of decimal digits to be printed, where m 4 12, then 
X = 12 - m (0 4 I ■$ 7) i.e. J decimal digits are ignored completely. The 
integer to be punched is multiplied by 10* before entering the printing 
sequence so that the first X significant digits are ignored. 



5. Zero Suppression 

There are two indicators, cr and ft, set during the subroutine, which control the 
suppression, omission and punching of digits. 

cr = causes the digit to be punched. 

a = 1 causes the digit to be suppressed or omitted. 

ft - there is no R.H.Z. suppression or omission. 

jj. - 1 the remainder, r, is examined after punching each digit; if it is zero, 
puts a .= 1 . 

i z i^ - 00 Start with a = 1, ft ~ 0. L.H.Z. suppression or omission. Set cr = and 
print as soon as first significant digit reached. 

i 3 i li = 01 As above, but print Sp or - before first significant digit. 

i 3 x 4 = 10 Start with a = 0, ft = 0. Print all digits. 

iji,, ~ 11 Start with cr = 1, ft = 0. L.H.Z. suppression or omission terminated by 
the 'special operation' group (i 5 , i &l i y ). See below. 

6. Special Operations 

Special operations are initiated at the end of a group if the group Indicators 
9 j t and g- 2 are both zero. 

i c Only consulted if i, =0. 

t 3 i g = 00 Treat the last digit of the group as a significant digit, 
even if it is zero. 

i 3 » 5 = 01 No special treatment for last digit. 

i 6 i 7 - 00 Set cr ~ 0, ft = to stop further zero suppression (or omission). 

i 6 i 7 ~ 01 Punch decimal point. Set cr - 0, ft - 1 to stop L.H.Z. suppression 
but allow for "R.H.Z. suppression later. 

i i 7 ~ 10 Punch decimal point. Set a = 0, ft = to stop further zero 
suppression. 

i 6 i 7 = 11 If remainder r = 0, punch Sp. Set cr = 1, ft = 1 for immediate 
R.H.Z. suppression. 

If remainder r t 0, punch decimal point. Set cr - and ft = 1 to 
allow for R.H.Z. suppression later. 
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NOTATION 

3 / 

R.H.Z. 
L.H.Z. 

a 
s 

h 

z 

I 

8 
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five least significant binary digits of accumulator 3. 

sign bit of accumulator 5. 

right hand zeros 
left hand zeros 

Controls zero suppression, see section 5. 
Controls R.H.Z. suppression, see section 5, 

Layout parameter, see sections 1, 5, 6. 

Group counter, see section 2. 

Group indicator, see section 3. 

Style number. 

Scaling constant. 

Layout-parameter, see section 1. 

the number to be scaled and punched. 

|£| = \ks | where k is the scaling factor. 

Next decimal digit to be punched where 8 = [£. IQr/lO ] 

Remainder after forming 8; replaces £ to form next decimal digit. 
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PRINTING CYCLE 

Contents of Accumulators on entry to Printing Cycle 

1st Entry 



x l 


x m = *x ?i 


X 2 


i 3 + «7 11 < ?12 ).2- 13 + i s .2" 38 


*3 


i H . OOOOOOt 7 00000i 6 00. . . 


Space 
_ or - 




*4 


jr, 4 S becomes /x = after 
*m - ~ u * i; 4 c ~ A > 66 or der in 5+.4+ 


*5 


Sign bit = a; 5 m = qy, 5 C = #>j 


X 6 


|C| where \t\ .= |*«| 


X l 


Overwritten 



Subsequent Entries 



*1 


V " * j q j 


*2 


i 3 + (qj! <7j 2 ).2" 18 +i 5 .2" 38 


*3 


i 4 . 000000^000001^00. . .0 


"Space" 
-or -_ 




*4 


Sign bit = pi 4 B = 4 = 


*s 


Sign bit = a; 5 m = qy, 5 C = pj - If 


*< 


Remainder r 


*7 


Overwritten 
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FLOW CHART OF PRINTING CYCLE 
Entry 



> Form S, next digit for printing 



Test for last 
digit in Group j 



last 



Test cr 
(digit suppression) 



a = 



not last 



c 



Test q 




a = 1 



Is t 3 = i s 



<?;'! = Ij* = °? 



Omit S 



Yes (Print 

last digit 
even if zero) 



Test i. 



Set cr = 

to stop L.K.Z. 

suppression 



J 



»» = 1 



Print sign 
before digit 



•/ Test fj. \ 



fi= 1 



/* = 



Test r \ 
remainder J r ^ q 



r = 



Set cr = l 

to start R.H.Z. 

suppression 

on next digit 



Print digit 8 
or Space 



No R.H.Z. 
suppression 



_>/ Test for end 



More digits in group 



A of Group j 



J 



End of group j 



o- = 
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01 



02 



03 



R 





-0 3 


1 


- 


28 - 


0+ 


m 


72 


0.2+ 


60 


0+ 


& 


72 


0. 3+ 


60 









10+ 


6 


00 0. 



Title of cue list 



Normal cue 



Cue used only by R3 



Address -0.2 of Parameter List * 



LN 

4> NUMBER PRINT 



R 


5-0 3 


1 


- 01 - 


R 4 


6-0 2 


3 


- 01 - 


R 8 


3-0 3 


1 


- 01 - 


Rll 


0-2 


1 


- 06 - 



Call for address of Parameter List 



Call for cue 02 to R3 or R03 t 



Call for address of Parameter List 



Title of Optional Parameter List 



* Note that if the optional parameter list is not used, the Interlude will alter cue 
03 to be the address -0.2 of the programmer's parameter list. 

t When the Master Programme calls for R3 this tag will cause the order in 4+. 6+ of Rl 
to restore B0+ of R3 to U0. R3 uses Rl as a subroutine and this is the means of 
returning to R3. When the Master Programme does not call for R3, this tag will 
cause the dummy subroutine R03 to be accepted. Cue 02 to R03 is such that the 
order in 4+.6 of Rl will then restore B0+ of Rl to U0. 
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(0+) 



0.0 



1 00 

0.1 6 10 



( nl 1 






073 


kjzL-s " 


0.3 


1 


6 00 

7 52 




(Used only 
by R3) Qi 

0.5 


2+ [ 
10+ [ 
0.0 


D73 
2] 73 
3 00 
1 |70 7 




lo 


0.6 


1+ [ 

7+ [ 


[]72 
3] 72 
5J 72 
D 60 


0.7 


8+ [( 
0.3 




1 


' 


(8+) 

0.3 


1 ', 


i 00 
)|70 7 




lo ( 


0.4 


1 ' 

4 ; 


t 00 
! 05 




0.5 


0,7 e 


) 41 
i 63 




0.6 




6 < 
© c 


i 02 
41 




6 I 
(30) 2 


21 
40 




0.7 


(1 + ) 


1 


f 




1.0 


1.3 4 
16 2 
16 ■ 


! 40 
63 
10 
12 




1. 1 
1.2 




17 C 
32 3 


10 

01 




1.4 1 


67 










, » 






' 


16 7 


12 










1.3+ 1 
©4 


66 
52 










1.4 






1.6 4 
16 3 


63 
10 




l.S 




©4 
32 2 

4 2 
32 5 


52 
00 
05 
00 




i.a 

1.7 
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Mask 

ENTRY 

Store Accumulators 

Number to print, s to 6 

2s to 7 ra , s = style number 

Store 01 

Store U2 

Mask (= 1.000000100 0) to 3 



+ cue 03 to M. Read ~k s to 1 



+ cue 03 



I s to 4 



| -* s to 5 
\l s to 1 



T .000000t 7 000 to 3 

3j= - 11 to print - 

Jump if z < 

Negate z If z > (I.e. -\z\ to 6) 

3yr = 14 to print Sp 

t, = |**| (rounded). If £ = +1.0, 6 = -1.0 and OVR 

2 C - +30 for printing CR 



7 C = +19 



Jump if i fl = 1 (I.e. no CR LF) 
Print OR 
Print U" 
Print 4> 

1.000000i 7 000 

7 e = +18. Unconditional jump 

Print (7 - S) Spaces. 

* t . l 3 i % i 5 i ( t T etc. to 4 

Jump if j : =1 (sign delayed) 

Print Sp if z positive, - if * negative. 



V\ *s i t, i i etc - t0 4 



i 3 . 000000 to 2 



-1.0 to 5 

*-to 7+.0 (sheet 8) 
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from 1+.7+ 




(8+) 







i__. 




7 


5 01 


1.3 











7 00 


4 


7 53 




35 


7 05 




7 


5 01 


5 


^— «. 






(o) 


7 40 1 




(T) 


7 45 


6 


p>=^ 






C 25 ) 


7 52 




8+ 


[2] 72 


7 








2.0+ 


2 67 


2.0 


0" 


7 41 


^ »- 


32 


2 05 




7 


2 01 


l 








9+ 


[l]72 




5+ 


H]72 


2 








2.3 


60 



V*S *6 *7 et0 - t0 4 

Jump if i n = 1 

cr = in 5 S If i 3 = 1 and i ^ = 0; ct = 1 otherwise 

i, = 0; i 4 . 000000i 7 000 to 3 

V*e *7 e 

Jump if i. 



i r 000000 .000i 5 to 2 



t^.OOOOOO^OOOOOi^O to 3 



0.000000 i t i 7 -*~ 

0.000000 0-* — 

1.111111111111100... 



-*-to 4 



-*-to 4 



-X *~ to 4 



NEXT GROUP 

Next group markers (#.• qA to l m 

ex. 000000 to 5 

0j 1j> to 7 m 

-38 



^.2 



to 7 







(6+) 



a + Pj.2~ 3H 



Ijo t0 \ 



0.7+ 



to 5 



4+ Q~|72 7 



Read B4+ 
to Ul if pi 



.To 'LINK' 
' (sheet 11) 



cr + qr 



JO- 



+ V 2 " 



to 5 



V ';2 t0 \ 



Jump if i ^ 



tyi ?j 2 ) 



.2" 



+ V 2 " 



to 7 



»i + (9j t <? J2 >-2 



+ i„.2 



to 2 



->_to 5+. 3 to Printing Cycle 
(sheet 9) 
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from 8+,2+ 



(sheet 8) 



(5+) 2.3 



(7) 4 41 

6 7 00 

6 00 

0.0 4 66 



PRINTING CYCLE 

1st time, -2" 13 + (X + 1).2" 3B to 4; \t,\ to 7. 
Subsequently, y, + 2" 3e to 4: r to 7 (r Is remainder) 



Jump to form next digit, with 4 C = 1 



I I 

I I I 

I 1st ENTRY (With 4- 



2.5 + 



r-i--, 

L 6 -°_ 2 i J 

1 i 

[0"»i 

I ! 
MULTIPLY I 

BY 10^ I 

__Jr. ! 



(9+) 



2.6 



1.3 + 



0.0 4 67 
1.3+ 60 



(5+) 



32 4 05 

2.7 6 24 

1.6 5 67 

2.1 5 62 

1.7+ 2 60 

2.1 5 62 

0.2 7 60 5 

0.2 2 63 5 

2.0+ 3 62 



l 2 -^ 



16 3 10 
32 5 03 





1 

0.2+ 4 62 

0.2+ 6 61 

32 5 01 

0.2+ 60 




1 


2 







(6+) 



14) 7 40 



o+i" 



16 7 10 



fj. - (6+) 

X) 

If I = -1 (OVR) 
6 = 2" 9e , 7 = 



Move up gap bit 
before p 20 



>C* 


0.0 6 61 


0.0 

l 


(To) 6 40 


6 7 20 


2.6 60 



STOP If 0VR 



6, 



+10 



lOr to PQ 



J [ Jump X times on 1st entry 

Pass through on subsequent entries 



jj. to 4 

Divide by 10 11 ; S (next digit) to 7, r to 6 

Count digits in group 

Jump to print if v = (last digit of group) 

Jump if t 3 = i 5 = qjl = q j2 = 

Jump to print if a- = 

Jump if S = 1 to 6+.2 if q ■„ = 

Jump if i = l) to 6+.3 if q i0 = 1 



r J° 



Jump if i^ =0 



Print sign before digit 11 i li = 1 

a- - 0. Print this and later digits 
(i.e. stop L.H. zero suppression) 



Jump if fi - (no R.H. zero suppression) 

Jump if r J- (number not yet finished) 

Put o- = 1 if /i = 1 and r = (start R.H. zero suppression) 

Jump to print last digit 



Replace S by Sp (If q a = 0) 



Print S or Sp 



'0' 



5 41 
2.3 5 67 



Jump except after last digit of group 



-*-To 6+.4 after last digit of group (sheet 10) 



from 6+.3 



(sheet 9) 
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(6+) 



0.4 



T$ 



7 42 



0.5 60 2 



7 m = -1, 7, = Sp 

Jump modified by q 1 .^ qr.- 2 



q. q, = 00 



0.5 



5 00 
1.0 7 67 3 



i 6 i 1 = 00 , 



l 6 t ? = 01 



»« *7 = 10 



*6 »7 =" 



Set a = 



7 = -1, 7- = LP Unconditional Jump 



(9+) 1.0 
(6+) 0.7+ 



0.7+ 60 



4+ 072 7 



T 



Read 3+ to Ul 



(9+) i.o> 
1. 1 

(6+) 0.7 



32 4 00 
0.7 7 67 



16 7 

4+ 072 7 

(33 



H 1 'm ~ - 1 

10 -* — > Print • 



Set jj. = 1 



V = -I. 7 C = • 



Read 3+ to Ul 



0+) i.r 

(6+) 0. 7 



0.7 7 67 
16 7 10 

4+ 072 7 



u. i i a i * 

16 7 10 ■* — ' 



T 



Print 



■1. 7 C = • 



Read 3+ to Ul 



(9+) 1.1+ 

1.2 
1.3 

(9+) l.o+ 
1. 1 

(6+) 0. 7 



_Q 



1.0+ 6 61 

32 5 00 

7 41 

1.0+ 60 



32 4 00 

0.7 7 67 

16 7 10 

4+ [I] 72 7 



T 



Jump if r + 

If r = 0, put a - 1 



1 n = -1, 7 C = Erase 



Put /j, = 1 

7 = -1 7 = 

Print • or Sp 

Read 3+ to Ul 



/• if r t 
(sp if r = 



9, <?, = 01 



0.6 



0.7 



141 7 40 
16 7 10 
16 7 10 
4+ HJ72 7 



"C 



7 m = °. ?c = Sp 
Print 2 spaces 
Read 4+ to Ul 



<7, <7 9 = 10 



■JUMP TO 'LINK' (sheet 11) 



0.7 



16 7 10 
4+ 072 7 



T 



Print Sp 
Read 3+ to Ul 



*1 «? 



= 11 



-T7T?" 



0.7+ 



0.7+ 60 



4+ 072 7 



Read 3+ to Ul 



To 3+.0 to Next Group 
(sheet 8) 
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(4+) 



1.0 



1.3+ 7 61 



1.0+ 



l. 1 



(0+) 



1.2 






1.3 





i 
I 


C 14 *) 


5 40 


(7) 


7 40 1 





7 45 


1.3+ 


7 60 


16 


5 10 


1.2+ 


7 67 



1.3 + 



1.4 



1.5 



[7] 72 

0+ [0)72 

10+ [2] 72 

1.7 6 63 



Computing 
Store 
LINK 



1.5 + 

1.6 
0.1 



0.2 

(4+) 1.6 + 
1.7 



12) 6 52 



0.1 60 



1.7 6 11 
U2) 6 53 
1.6+ 60 



3 



[0]72 
2+ [T]72 
60 



J 



Sheet 11 of 14 



Jump if P t 

5 C = +14 (Sp) 

If q = 0, clear 7 
Jump if q = 

Print 9 Spaces 



Restore the Accumulators 

Restore U2 

Jump if LINK a go-order pair 



(4+) 



1.7 



(0+) 0.0+ 



0.1 



2+ [TJTC 
0.0+ 60 



0.1 6 10 




Restore Ul 



Overwritten 
by LINK 



Form A 60 in 1 . 7+ where A is the address 
set in the computing store LINK. 



Restore 6 



+ m» rw +n rt Read B0+ of R 3 1 +„ tin + 
+ cue 02 to R3. Qr BQ+ Qf R1 J to UO t 

Restore Ul 

Obey computing store LINK 



t See note on sheet 6. 
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■*-!«<- 



-*■(■*- 



^ = 



*fl 


, 




J S" w 


- * . 


^OO^Ol^Oi ^OO 'dl '02 


*1 U l 3 U *B l « *7 









*10 *11 *1« *10 'll *1* 


^20 * '»!' *22 '20 ^21 ^28 


*30*Si ^32 ^30 ?S1 ^82 



T 11+. 





Fl 




1.2 71 0. 


01 






5.2 6 10 


02 


-100000000000 




1.2 71 


03 






5.2 6 10 


04 


-100000000000 




1.2 5 71 


05 






6.1 00 


06 


-1000000 




1.3 76 4 


07 






6.5 00 


10 


-1.0 



Optional Parameter List 



-k.. 







12+ 


11 


1.0 77 7. 
5.4 20 


12 


-1.0 


13 


1.0 76 2 
6.0 20 


14 


-1.0 


15 


0.0 7 72 7 
3.0 20 


16 


,. 


17 


0.0 7 73 7 
2.0 20 


20 


-1.0 



-k. 



-k. 



-k. 



-k. 
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13+ 




0,0 

1 

2 




7. 7 7 00 0. 






0+ 00 0. 

2 00 0. 


from Input 




J560.0+ - 13+,! 


! 3 
1 

5 
6 


33 4 

34 4 


00 
03 




10 5 


70 4 










[1 
1 2 

-0T 


70 1 

00 

45 




7 


5.5 5 00 
(T) 1 52 












14+ 




1.0 
1 
2 
3 

4 


0.0 1 05 
0.2 1 03 










|1 
34 7 


71 5 
00 




34 7 01 
5.7 7 13 




0.1 1 00 
1.6 2 00 










|11 1 


71 2 


tQ Input 


5 


(-") 




6 




0+ - 00 




T 13+ 

J 560. 0+ - 13+. 

L 



R 0-0 2 
03 - 28 - 



INTERLUDE 



Collating mask 



Cue 02 for R03 



Plant LINK 

(511.0, 0) to 4 

Index address from 511.5 to 1 

Index word to 1 

Index word to 2 

Clear 2 if no parameter list 

Jump to obey LINK if no parameter list 

Reference Address to 5 

Parameter List address to B. P position of b-order 



Form new cue 03 which is the Parameter 
List address less 2. 
Store new cue 03 to Rl 



Set Transfer Address equal to 11+.0 so that the next 
subroutine overwrites the optional Parameter List 



Plant cue 02 for R03 in 11+. 1, since R03 will, if accepted, 
overwrite the unwanted optional Parameter List 



Restore LINK 

Overwritten by LINK for return to Input 



T 0+.2 

L 

L 



Cue list title to satisfy the call for R3 made by Rl. 



Indicate to Assembly that R03 has 2 cues 
Terminate cue list to R03 
Terminate subroutine R03 
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Ro 0-03 


T 3 + 












1 -38- 














0+ 073 


6 153 


o.o 


661 


0.7+060 


O 




o. 3+060 


33 505 


10 


640 


33 400 


137 


7000. 


0+ 073 


1 700 


6 


730 


0.7 767 


O 




0.3+060 


38 753 


3.6 


060 


1. 0+661 


O+OOOO. 





7 745 


14 


740 


33 500 


O 




ro+ 6000, 


0,7+760 


16 


710 


3 741 


O 


3000. 


LN 


7 5oi 


1 


541 


1.0+060 


liS 


no 


NUMBER 


PRINT 1 700 


3.3 


567 


33 405 


33 


400 


Ro 5-03 


3 753 


18 


743 


3.7 634 


34 


4<>3 


1 -01- 


35 705 


0.5 


0603 


1.6 567 





5704 


R4 6-oa 


7 5oi 





500 


3.1 563 


1. 


0701 


3 -oj- 


7401 


1.0 


7673 


1.7+360 


1 


300 


R8 3-03 


3 745 


14 


740 


3*1 563 


4 


345 


1 -01- 


35 753 


16 


710 


0.3 7605 


I.4+360 


Rll 0-30 


8+373 


16 


710 


0.3 3635 


5»5 


500 


I -06- 


3.0+367 


4+1737 


3.0+363 


5 


153 


I OOO 


1. 3+761 


I 


453 


T IX + 


0.0 


I05 


0.1 610 


14 540 


3.3 


4^3 




0.3 


I03 


1.7 611 


7401 


3 


506 


1*3 07I0 


I 


07*5 


la 653 


7 745 


33 


303 


5.3 6ioo. 


34 


700 


1.6+060 


1.3+760 


X 


453 


-100000000000 


34 


70I 


773 


16 510 


3.3+463 


1.3 07I0 


5-7 


713 


1 600 


1.3+767 


I 


341 


5.3 6100 


O.I 


IOO 


36 753 


773 


I 


453 


-100000000000 


1.6 


300 


3+173 


0+073 


3.5 


463 


1.3 57IO 


II 


1713 


10+373 


10+373 


35 


301 


6.1 0000 


«*5 


IOO 


o.o 300 


1.7 663 


34 


453 


-1000000 


+0 




1707 


13 653 


7 


445 


1.3 0764 


0+ 


-OO 


1+173 


o.i 060 


35 


403 


6.5 0000 







7+373 


O 073 


6+073 


-I.O 


T13 + 


8+073 


3 + 1 73 


3 + 173 




JS6o.o+-i3+.3 


0.3 060 


0.0+060 


I.O 


060 


t.o 0777 
5.4 oaoo* 


L 
Ro 


0-03 


19 74o 


16 310 


1 


74i 


-1.0 


03 


-38- 


i»3 4*53 


33 5°3 


33 


305 


1.0 0763 


T0+.3 


16 310 


0.3+463 


7 


30 1 


6.0 0300 


L 




16 713 


0.3+661 


9+173 


-I.O 


L 




17 010 


33 501 


5+373 


o.o 7737 






32 301 


0.3+060 


3.3 


060 


3.0 oaoo 






1,4 767 


1 44* 


1 


500 


-1.0 






16 713 


6 700 





0707 


0.0 7737 






1,3+166 


600 


1 


400 


3.0 0300 






1 453 


o.o 466 


4 


3°5 


-I.O 






1.6 463 


6 023 


II 


34i 








16 310 


1 15* 


0.7 


663 








3 4S« 


o.o 467 


6 


60a 








33 300 


1.3+060 


3 


34i 








4 305 


+100000000000 


6 


531 








33 5°° 




30 


340 
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PAGE LAYOUT SUBROUTINE 



This subroutine uses R 1 to print a number from XI. R 3 itself arranges the 
numbers in rows, columns and blocks with all necessary CR and IF characters. R 1 is 
concerned only with printing the number in XI in the manner directed by the layout 
parameter, the number of which, s, must be set in X7 on entry to R 3. 



The preset parameters L, B, N determine the layout: 



L 
B 

N 



No. of numbers in a line 
No. of lines in a block 
No. of numbers in a block 



1 If B f 0, 
J 5 = 0, 



then N = 
then N f 



01 



02 



R 


0-0 2 


3 


- 28 - 


0+ 


00 0. 












0+ 


00 0. 



LN 

<p BLOCK LAYOUT 



R 1 


1-0 3 


3 


- 06 - 


R 3 


7-0 2 


1 


- 01 - 


R 1 


5-0 1 


3 


- 02 - 


R 1 


6-0 2 


3 


- 02 - 


R 1 


7-0 3 


3 


- 02 - 



Title of Optional Parameter List 



Call for cue 02 to R 1 



Call for P. P. 01 



Call for P.P. 02 



Call for P. P. 03 
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Note: m is a marker set during R 3 

bi = indicates that a line is unfinished 

m positive (= 13 x 2" 38 ) Indicates start of new line 

m negative <= -1.0 + 13 x 2~ 3 ') indicates start of new block. 
At the end of a line or block, where m £ 0, the five least significant bits of m contain LP 
for use by the master programme if required. 



Jump if 7 C = +s (s is style number for R 1) 

|s| to 7 

Preserve Accumulators 

Set counters etc. in accumulators 

m + 1 to 1 (to test for m .= 0) 

Jump if LINK positive 

Negate LINK 

Preserve Accumulators 

Set counters etc. in accumulators 

Bet marker, m = -1.0 at start of new block 

Preserve Accumulators 

Set counters etc. in accumulators 







0+ 


Normal 


0.3 7 62 
7 7 02 


Entry 


0,0 

>- 

1 

2 


New 


[7] 73 
1+ [7| 72 


Line 




141 
0.7 60 










0.6 6 62 
6 6 02 










— 


1 
3 

4 

5 


New 


o El 73 

1+ El 72 


block 




32 1 00 
0.7 60 




a 




o El 73 

1+ 072 






2+ [0)72 


to 2+.0 


7 








U. O 


DU 




Optional P.P. 01; L = 4 



Optional P. P. 02; B = 5 



Optional P. P. 03; H = 



' becomes m 



becomes L 



• becomes B 



• beeemee # 



* p.p.ei; i 



* P.P.ea; 9 



* P.P.OSj n 
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* The counter B counts the number of lines in a block. Normally B will be reduced by 1 in B3+.1+ 

before printing the last number in the line. If the Master Programme requires to start a new line 
before the previous line was completed, it nay enter at 0+.1 (or at 0+.0 with -s in 7). On entry m 
will be zero, the a-order in 2+. 3 will not then Jump and the counter B will be reduced in 2-t-.3+. 



from 0+.7+ 



+0 



Jump if m - i.e. middle of line 

Print CRLF at end of line 

Jump if m negative i.e. end of block 

Jump to start new line * but not new block 



Print extra W 

Reset 3 and 11 at end of block 

Reset L at end of line 
Print </> 



► New block 



Overwritten 
by m 



Set w - 

Count numbers on line (L) 

Set m - 13 x 2" 3B at end of line (LF) 

Count lines in block (B) 

Set ix = -1.0 + 13 x 2" 3B at end of block 

Count numbers in block (If) 

Set 7i! =-1.0 + 13 x 2~ 3e 

Store counters and marker in Bl+ 

tn to o.o 

Restore Accumulators 

Number for printing to 6 

Write m to B0.1 (R 1 leaves m in XI) 
+ cue 02 to R 1 
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NUMBER PRINT (SHORT) 



This subroutine prints the contents of X7 as a signed or unsigned integer or 
fraction. It uses the Initial Orders number print routine. 





R 0-0 9 
4 - 28 - 


01 


0+072 
0.3 60 


02 


0+072 
0.0+ 60 


03 
04 
05 


1+072 
0.5 60 

1+072 
0.2 60 

0+072 
0.7 60 


06 
07 


0+072 
0.5 60 

1+072 
0.4 60 


08 


0+0 72 
0.2 60 


09 


1+072 

0.0+ 60 



LN 

I.O. NUMBER PUNCH 



R 3 


0-0 4 


4 


- 06 - 


R 


1 -0 2 


4 


- 02 - 


R 


3-0 1 


4 


- 02 - 


R 1 


2-0 4 


4 


- 02 - 


R 1 


5-0 3 


4 


- 02 - 



Title of Optional Parameter List 



Call for P.P. 02 



Call for P. P. 01 



Call for P. P. 04 



Call for P. P. 03 
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0+ 




0.0 


| 449 - 00 
073 


f rp V- 


{ VJ ) w 




1 


© 2 40 
0. 5+ 60 






( no *\ 


073 
1+072 


I. u » ) 




2 








C~K\~\— 


073 

© 2 40 

(677) 6 40 

0.6 60 






3 








4 










( Q6__)-- 


n FtI 73 








5 


















yu. 3+jo ^v 






0.0 4 00 
1+072 






6 








C (,= v 


73 
0.4 60 






7 














1 + 




0.0 


| 448 - 00 

073 


f no V 


^ UJ ) 




i 




(O) 6 40 
0.6+ 60 






^ 




073 
© 2 40 


( ()4 V. 


V. ui J-"- 




2 

v 




(6.i'-i) 6 40 

0.6+ 60 






3 








C 07 V- 


073 
0.6 60 


^ u/ J- 






4 


( 03 V- 


073 

© 2 40 

{fj\ 6 40 

0.0 4 00 


\. " c J 






5 

> 

6 










2+0372 
60 


J;o 2+.0 


1 















= (449.0) Modifier for printing fraction 
Store Accumulators 

+ P. P. 02 \ c 

J Jump 

Store Accumulators 



+ P. P. 01 



Store Accumulators 
2„ = n, 



6 =0.7 To print sign 



Store Accumulators 



6„ = 0.3+ To omit sign 
c 

4 = 449.0 (Modifier for printing fraction) 



Store Accumulators 



= (448.0) Modifier for printing integer 
Store Accumulators 



6 = 0.3. To print 3 digit unsigned number 



+ P. P. 04 



Store Accumulators 



. 2 = n 
c n 



6 =0.3+ To omit sign 



Store Accumulators 



+ P.P. 03 . 



Store Accumulators 



2 „ = «, 
c 3 



6 =0.7 To print sign 

4 = 448.0 (Modifier for printing integer) 
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| 2+ 


from 1+.7+ 


0.0 5 02 
@ 42 2 
(g>) 6 52 
0.3+ 1 63 




0.0 

I 




2 
3 




1 1 02 




1 l| 71 




(T) 5 53 
16 5 10 






4 


(?) 5 53 
16 5 10 


5 


0.7 1 00 
33 3 00 


e 


28 072 3 

60 6 




7 


77 072 4 
0.2 60 



Prepare to print Space or CR LP 

Dummy. Constant used to print Sp Sp or CR LP 

Form modifier for style of printing 

Jump if LINK set by 00 order 

Negate 1 and Store LINK 

if LINK set by 02 order 

Prepare to print CR LP 

Print CR or Sp 

Prepare to print LF or Space 

Print LF or Space 

Set Initial Orders LINK 

512.0 to 3 

Enter Initial Orders number 

print routine in B540 

LINK to enter 1.0. L- sequence in B525.2: 
this restores the Accumulators from B 
and obeys the Main LINK 





| 3+ 







11 


00 0. 


12 


00 0. 








12 


00 0. 







12 


00 0. 



OPTIONAL PARAMETER LIST 



P.P.01; n i .= 11 



P. P. 02; n = 12 

2 



P. P. 03; n ~ 12 

3 



P.P. 04; n 4 ~ 12 
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DOUBLE LENGTH NUMBER PRINT 
Flow Diagraa 



Sheet 1 of 5 



Preserve Accumulators 
Print CR LP 



Calculate round-off 

constant (r. o. c. ) 
Restore Accumulators 



q >0 




Set up + 



Set up - 



Print sign, 0; 
Add r.o.c. ; Test OVR 


' 


OVR set 


\ 


Put q ~ (1 - 2" 38 ) 





<D ( 



£ 



Test if n 1 & n 2 
are in range 



Yes 




Test n. 



n l 1- 



OVR clear 



Print • 
Print q 



02 



Preserve Accumulators 
Print Sp. 





LOOP STOP if 
n x > 12 or n 2 > 11 



Test n 



»! = 



STOP 




n 1 ? 



Test q 



} 



q < 



* STOP 



q > 




Test p 



P > 



p < 



Set up + 



Negate p and q 
Test OVR 



Clear OVR 
Add r. o. c. to q 



OVR set 




q > 



Set OVR 



Print p 



Restore 

Accumulators 
EXIT 
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R 








-0 2 




5 


- 


28 - 




0+ 





72 











60 




0+ 


@ 


72 





2+ 





60 
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LN 

D.L. PRINT 









0+ 


( m ^ 




073 
(~30) 3 40 


\ ox ) 


0.0 

i 

2 




16 3 10 
Q9) 3 42 




0.3+ 60 


( 0" \ 


073 


I °- ; 


3 

> 

4 
5 
6 




fTi) 3 40 
16 3 10 




10 10 3 
(l27) 4 45 




(127) 5 45 
(jD 4 44 


0---- 


f"Ti) 5 44 

0.6+ 4 63 




1+072 
0.0 60 



Write Accumulators to BO 

Print CR 

Set up LP and set 3a = 7 

Write Accumulators to B0 

Print LF or Sp. 

Dummy order, or Print <p after CR LP 

Remove 32 m.s. bits of 4 and 5 



1-2-Mj to 4 C 



ll-n 2 to 5 C 



LOOP STOP if b l > 12 
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from 0+.7+ 






1+ 


(7) :s^ 


0.0 5 63 
(T2) 4 44 


^y 


(D. -- 


QS 5 44 
0.3 5 60 


\zJ 


2 


2+[5 






72 




3 




©___... 


0.3 4 60 
7 00 


(T)._ ._ 


2 00 

0.4+ 7 63 


\ZJ 


5 


(7).. „__ 


fll) 3 40 
0.7 6 62 


\ZJ 


6 

7 


to 3+.7 


6 6 02 
3+0 72 




to 3+. 2 


3+0 72 
0.2 3 67 


















2+ 




0.0 


+0.1 




33 6 00 

0.0 6 21 




0. 3 

1 *" 

2 
3 

4 


(7 s ) J - 


/0.1+ 5 67 \ 

\0. 2 6 10 J 

072 

0.2 2 00 
(l27) 4 45 

0.7 4 60 


w- 


5 

» 

-,6 


©■--. 

to 4+. 2 


1+0 72 
0.4+ 60 


4+072 
0.2 7 63 














\ 0.6 t 


J 60 



Sheet 3 of 5 



TEST n 



and q 



LOOP STOP if n 2 > 11 
rc : to i c 
n 2 to 5 C 
Test n„ 



Jump if n 2 ? 

LOOP STOP if « 1 - n 2 = 

Clear 7 

Clear 2, rounding off constant = 

LOOP STOP if q < 

Set up - 

Jump is p > 

Negate £, if ^ < 

Jump if p < 

Set up + and jump 

FORM ROUND OFF CONSTANT 



Calculate round off constant (r.o.c. ) 

r. o. c. to 0.2 

Restore Accumulators 

r.o. c. to 2 

Remove 32 m. s. bits of 4 

Jump if n 1 = 

Jump if n l ¥ to print fraction 

Jump to 4+. 2+ if q < 0, to 4+. 7 if q > 
Set up - 
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3 + 


PREPARE TO PRINT INTEGER 

Test OVR 

Test q if OVR set, jump If q ^ 

Justify 

Clear OVR if p = -1.0 and q f on 

Add r.o. c. to q 

Justify 

Set OVR 








0.1 64 
0.3+ 7 62 


















0.0 






- 










6 23 
0.2 65 








1 


entry 


from 1+.7+ 






2 7 01 
6 23 






2 










0.4 60 








3 






32 02 


















7 2 00 
6 7 00 








4 


' 


q 10 ^ 
P to PQ 






5 


6 00 
6 23 








6 


4+[o]72 
0.0 60 


Jump to 4+. 0. to print p 
Negate q if p < 




from 1+. 6+ 


- 






7 7 02 
0.0 60 






7 




























4+ 


PREPARE TO PRINT 

Set 12-^ in 4 C 

Set 12 in l c 

Make C(3) a GO order pair 

Jump to print Integer 

Set up - 

Negate q (q < 0) 

Print sign 

Print 

Add r.o. c. to q 

Test OVR 

I If OVR set, (1 - 2" 38 ) to 7 

Jump to print Fraction 
Set up + 








(~V2) 4 44 
fl2) 1 40 






0.0 




to 5+.2 l 


32 3 01 

5+[0]72 




from 2+. 6+ 




(H) 3 40 
7 7 02 




2 








16 3 10 
16 10 






3 








4 


2 7 01 
0.6 64 








5 






(7) 7 42 
32 7 03 






6+072 
0.1 60 




to 6+.1 




6 




from 2+. 6+ 






Qo) 3 40 






7 












U. d U 


DU 
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R 5 



5+ 



from 4+. 1+ 



I 

i 0.0 
i 
i 
->i 



l 



from 4+.6+ 




0.1 7 62 
6+(o] 72 



+100000000000 



6+ 



O.o 



*t 2 7 00 

0.4 5 60 

(Tg) 6 40 

16 6 10 

(75) 6 40 

6 7 20 

16 6 10 

0.2 5 67 

72 

0.5 1 10 

+0 



+0 



(127) 5 45 
0.0 60 



Sheet 5 of 5 



PRINT INTEGER 

Print Sp for non-significant zeros 
after (12-?^) cycles 



(Remainder) x 10 to PQ 
Next Decimal Digit to 7 
Sets 4^, = -1 after (12-Kj) cycles 
Count 12 decimal digits 
Ensure printing of last integral digit, even if 



Jump if digit is zero 
Print sign before first digit 
Cancel sign printing 
Print digit of p 



= 10 11 x 2" 38 



PRINT FRACTION 

Restore q 
Jump if « = 

Print • 

Print fraction 

Restore Accumulators 
Plant LINK and EXIT 

Overwritten by LINK 

Unused 

Remove 32 m. s. bits of 5 



becomes { 
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PAST R 5 



FERRANTI LTD 

PEGASUS LIBRARY PROGRAMME 

Sheet 1 of 7 

FAST D.L. PRINT FOR 7168 WORD STORE 

This routine has the same specification as the original R 5 but it makes use of 
the full speeds of the Creed 3000 and Soroban punches. The extra space required to 
store such a routine is obtained by using some of the extra isolated blocks on the 
7168 word store and the routine in this form can therefore only be used on that store. 

The routine stores the accumulators, Ul and U2 in B1+, B2+ and B3+ respectively 
of R 5. This makes Ul and U2 available for use in the routine without altering the 
original specification of R 5. Four isolated store blocks are used for programme and 
two for rounding constants, 

Most mixed numbers printed by R 5 are less than 10,000. The new routine 
separates such numbers in order to prevent going many times round a time-consuming 
zero omission routine. 

Error Analysis 

Representing the integer to be printed by p - I. 2~ 3B we form 

/ 2™ \ 2 7 " 

I* = I.2' 3B ( — — + e J .2 -36 where e is the rounding error due to terminating — — 

Vl ° ' / 2 ™\ 10 

after 38 bits and has the approximate value +.214 if we take f — r- ) = 188894659315. 

\10 u / 

I* is then multiplied by 4 and 1.2" 38 added to the result, all succeeding digits 

I 
being dropped, i.e. we form — — + ie.1.2" 76 + (1 - S).2" 3a where S is the error 

10 11 

due to dropping digits (0 ^ $ < 1). 

During printing I + 4e7. 2" 76 .10 u + (1 - S)2" 3a ,10 11 is formed. 

We require < (4el. 2" 3e + 1 - S)2" 38 .10 11 < 1 if the process is not to fail. 

Since S < 1 and I is made positive before printing, we ensure that the error term is 
positive by making e > 0. 

The error will be greatest when I - 2 38 (after rounding) and we therefore require that 

(4e + 1 - S) < 2 38 .10" n . 

In the worst case S = and ie < 1. 

Thus <4e + 1 - S) < 2 and since 2 38 . 10" u 2 2.74, the required conditions are satisfied 
and the integer printed by this process will always be correct. 
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■" 






R 


0-0 2 




5 


- 28 - 




0+ 


H]72 


01 








0.0 


60 




0+ 


[0]72 


02 








0.2+ 


60 



LN 

PAST D, L, PRINT 











0+ 


Cqi~) 




1+ 

® 


\l\ 73 
3 40 




0.0 

l 




© 


1 42 




2 


16 


3 10 




0.3+ 


60 


C n" ^ 


1+ 


[7] 73 


K ° 2 J 


3 




© 


1 40 




4 
5 


2+ 


[D 73 




® 

@ 


4 45 

5 45 
4 44 






3+ 


[2] 73 




4+ 


[1]72 




c? 


0.6+ 


4 63 




16 


1 10 


to 4+.0 


1 


© 













Preserve Accumulators in B1+ 



Set up LP and set 1 j, = 7 
Print CR 



Preserve Accumulators in B1+ 



Preserve Ul in B2+ 



Remove 32 m.s. bits of 4 and 5 



12 - n 1 to 4 C 



Preserve U2 in B3+ 



LOOP STOP if n L > 12 



Print LP or Sp 



11 -n 2 to 5 C 
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PAST R 5 



from 0+.7+ 



1.0 



to 1018.4+ 



to 1018.4 



to 1018.1+ 



to 1018.0+ 



5+ 



- 1.0+ 64 



Sheet 3 of 7 



TEST n, , n„ and q 



Set 3 n - 896.0. 3 C = +10 

Store p in 1.0 
Store q in 1.1 

Set 7 m = (896.0 + n 2 ) 

Read round off constant from B( 102 1.7 + n 2 ) n 2 > 1 

Restore p 

Restore q 

Jump if nj t 

B1018 to U0 

Jump if q > and n 1 = 

Set up - otherwise 



LOOP STOP if n 2 > 11 



"i t0 4 c 



n 2 to 5 C 



Print figure shift if cue 01 



Jump if n 2 ^ 



STOP if n 1 = n 2 



Clear q for n 2 - 

Clear 1, round off constant = 

LOOP STOP if q < 

B1018 to Ul 

Jump if p ^ 

Negate p, if p < 

Negate q, if p < 

Set up - 

Jump if OVR clear, p i -2 3S 
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PAST R 5 



PREPARE TO PRINT 



from 5+.7+ 



if OVR set 
from 5+. 7+ 
otherwise 



from 5+. 5+ 



1.0 




1018 



1.3 7 62 

6 23 

1.1+ 65 

1 7 01 

6 23 

1.3+ 60 



32 02 
124 [T] 72 3 



1.0 



to 1021.0 1 3 





1020 


+188894659315 


1.0 7 20 


(V) 54 
(7) 6 41 


6 7 00 


111 [T]72 3 


0.0 60 


7 2 00 


6 7 00 


1.0 6 00 


(n) 1 40 


125 [0]72 3 


(U) 4 44 


1.1+ 65 


1.1 60 
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Jump If q = and OVR set to print -2 39 

Justify 

Clear OVR if p = -1.0 and q t 

Add r.o. c. to q 

Justify 

Reset OVR if p = -1.0 and q - 

Read B1020 to Ul and jump 

Negate q if q < and n^ - 

Print sign 

Add r.o. c. to q 

Jump if OVR clear 



May be 

overwritten by LINK 



Form (1 - 2" 38 ) in 7 

Print 

Jump to print fraction 




.36 



+ € .2 



10 



11 



(Pq)=2 2 x 2- 38 J x i-n 
v / 10 n 



I' = — (rounded up) to 7 [see note on Sheet 5] 



Read B1007 to Ul 

Jump to print integer 

q to 2 

p to 7 

Conversion constant to 6 

Set 11 in l c 

Read B1021 to U0 



11 - ni in 4 C 



Jump if OVR set. i.e. (pq) = ±2 S 
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Note: On entry to B1021 

If I < 10 11 

If 10 11 «S I < 2.10 11 
If 2.10 11 < I < 2 3e 



7 5 = and OVR clear 
7 S = 1 and OVR set 



7 S = and OVR set 







1021 


from 1020.3+ 


(T) 5 52 

(l4) 5 41 

6 00 

6 23 


0.0 

1 




2 


7.7 6 12 6 
1.3 6 61 




3 




7 6 00 

2.6 6 03 
1.0 6 62 

(T) 4 43 
fO.7 4 62^ 
^16 5 10 J 
(T) 4 41 
0.5+ 4 61 
(7) 1 40 

2.7 7 20 






4 








5 

/ »» 

6 






7 














1007 






9 5 10 4 
2.5 7 20 




1.0 






_ 


1 


1.2 1 67 

1.3 60 




- 




2 


(T) 4 43 
1.0 6 60 








16 3 10 

1.5+ 1 60 






3 










f *" 

4 

5 


16 6 10 

2.5 7 20 




1.4 1 67 

(T) 5 53 

16 6 10 

2 7 00 




S 


tj 6+.2+ 










7 








PRINT INTEGER 



} 



Set 5, 



Set x. 



n 2 .2* + Sp 




1 
-2 



if 

If 

if 



I < 10 11 



10 ' 



^ I < 2.10 



11 



I > 2.10 11 
Set x 6 = +2 if I > 2. 10 u 
Jump if m.s. digit present 

I - 10 4 
10 11 

Jump it I > 10" 



4 -Hj to 4 C 



Jump if n 1 ^4 



Print Spaces if 

n l > 4 and I < 10 4 



Set 1 = 4 

" =(£) 



to 7 



May be 

overwritten by LINK 

I< 10 4 



Print Space (4 -, = 7) or dumpy order (4>, = 0) 
101' 



Reduce counter in l c ; initially 



Jump to print sign if last digit 
Count zero omissions 

Jump to print Sp if non -sign if leant zero 
Print sign on finding non-zero digit 
Jump if last digit 

Print significant digits 

Reset 5 C = n 2 

Print last digit 

Restore F (fractional part) to 7 



Jump to EXIT if n 2 - 



l c = 11 if I > 10 4 
l c = 4 if I < 10 4 









6+ 


from 1007.7+ 


(l2\ 6 40 
16 6 10 




2.0 
1 




2.5 7 20 
16 6 10 




2. 1 5 67 
1+ 072 


from 1007.2+ 


2 




3 
4 




2+ [Tf 72 
0.5 1 10 


to LINK 


3+ [|]72 
0.5 60 




5 
6 

1 




+10 




+ 27487 




+10000000 



r Hi K B A « T I LTD 
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PRINT FRACTION 



Print • 

Print fraction 

Restore Accumulators 

Restore Ul 

Plant LINK 

Restore U2 

Jump to obey LINK 



10" 
10 



li 



PAST R 5 



Sheet 6 of 7 



+10 '.2 



7 n-38 
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PAST R 5 







1022 









+13743895347 


1 


+1374389535 


2 


+137438953 


3 


+13743895 


4 


+1374390 






5 


+137439 


6 


+13744 


7 


+1374 



ROUND OFF CONSTANTS 



+0.05 for n 2 = 1 



+0.005 for n z =2 



+0. 0005 for n 2 =3 



+0.00005 for n 2 = 4 



+0.000005 for n 2 =5 



+0.0000005 tor n 2 =6 



+0.00000005 for n 2 = 7 



+0.000000005 for n 2 = 8 



1023 



+137 



+14 



+1 



+0 



+0 



+0 



+0 



+0.0000000005 for n 2 = 9 



+0.00000000005 for n 2 = 10 



+0.000000000005 for n 2 = 11 



► Unused 
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PRINT FLOATING POINT NUMBERS 



01 



02 



03 



04 



05 



06 



R 





-0 6 


11 


- 


28 - 


0+ 


@ 


72 


0.2 



0+ 





60 

72 


0.4+ 





60 


0+ 


® 


72 


0.2 


7 


66 


0+ 





72 


0.4+ 
0+ 


7 



66 
00 0. 
















0+ 





00 0. 



LN 

F. P. PRINT MK 



R 


1 -0 1 


11 


- 02 - 


R 


1-0 1 


11 


- 02 - 


R 2 


7 -0 1 


11 


- 02 - 



Calls for P.P. 01 



T1+ 




Optional Parameter list 



Title 



P.P.01 



Issue 3 



FERR ANT I LTD 
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R 11 



NOTE: 



number of digits to be printed before decimal point, 
number of digits to be printed after decimal point, 
decimal exponent. 



Sheet 2 of 6 



T0+ 



GD_ 



QD 



on 



CK) 



0.0 



0+ 



r 0.0 6 10 
,0.0 60 



+0 



073 
QD 4 40 



16 4 10 
(TIP 4 42 
0.6 60 



@73 

<JA) 4 40 

16 4 10 



1+-0 73 
1 6 00 



0.1 6 05 

2+ GO 72 



Plant LINK 

Jump to plant and then to obey the LINK 

+ (2 x P. P. 01) 

Preserve accumulators 

Print CR 

Set up LP and 4p = 7 

Preserve accumulators 

Print Sp 

Preserve Ul 
Number to 6 
Argument A to 6 



T2+ 



to 3+ .5+ 



to 3+.2 

+ 



to 3-K5+ 



1.0 



2+ 



6 1 03 

7 3 00 



2 00 

3+0 72 



16 4 10 
0.5+ 6 60 



1.7 1 01 
0.2 1 62 



ClO) 1 41 
(T) 2 43 

0.0 6 20 
(T) 1 56 



1.4 1 63 
0.5+ 60 

+0 



BINARY/DECIMAL CONVERSION 



Exponent a + 2 n_1 to 1 

n to 3 

Zeroise decimal exponent 



Print Sp or LF 
Jump if argument zero 
Subtract 2 fi_1 to form a 
Jump if exponent ^ 



Convert negative binary exponent in 1 to 
decimal exponent in 2. 



+ P. P. 01 
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3+ 


0.0 

1 
from 2+. 3+ 


+. 9765625 
+.024 

CTo) i 43 

0.5 1 63 




2 




3 

4 




(7) 2 41 
0.1 6 22 




(?) 1 56 
0.2 60 




CIO) 1 41 

4+ H] 72 


from 2+.2+ 


5 


and 2+. 6+ 






1.2+ 1 60 
(7) 1 43 




i • 

6 

7 




1.7 6 20 
(?) 1 56 












4+ 




1.0 

■ 

1 
2 


(T) 2 41 
0.6 1 62 




1 1 02 

Cm) i 52 

(7) 6 51 1 
10 10 4 




3 


C~1Q) 4 40 
1.5 6 62 




4 

^ > 

5 




6 6 02 
(~TD 4 40 




5+ 072 
6 7 00 


to 6+.0 6 


6+[T]72 
1.0 60 




7 


+.8 



BINARY/DECIMAL CONVERSION 



- ffi 



2^1 



Convert positive binary exponent in 1 to 
decimal exponent in 2 



Final conversion of binary exponent in 1 to 
decimal exponent in 2 



Dummy order, or print after CR LP 
Set up + 

Negate 6, if negative 
Set up - 

Argument to 7 



10 
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0.0 



to 1+.2+ 



from 4+. 6+ 



1.0 



to 7-K3 



6+ 



0.1 4 10 

7+ [l] 72 

1.3 60 



ADD ROUND-OFF CONSTANT 



= +0.1 + 2 



■38 



d' - n t (floating) or d' = d (fixed) 
Store sign in 0.1 



Store argument In 0.2 
n 

2 

d> + n 

2 

Jump if number effectively zero 



Form round-off constant 



Jump unless round-off causes OVR 



Subtract 0. 1 from argument 
Jump if argument > 0.1 

Multiply by 10 

Reduce decimal exponent by 1 <-ci) 

("j - d) 

(d - n ) to 1 as exponent 

Jump if fixed point and d 4 n 

Clear 2 

Mark floating point 

Jump if argument not zero 

Write sign into 0. 1 if argument zero 



from 5+. 7+ 



from 6+. 7+ 



to &f. 



6+ 
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l.o 



7+ 



0.0 6 00 
1.2 3 63 



(T) 2 41 
1.2+ 60 



©141 
2 4 02 
2 5 03 2 
5 43 4 







1.6+ 5 63 
1.6+ 5 60 



QX) 7 40 

16 7 10 

1.5+ 5 67 

0.1 5 00 



8+ 072 
0.1 60 



0.0 





8+ 


448 - 


00 0. 


14 




6 7 


00 


16 5 


10 


9+0 


72 


( 10) 5 


40 


16 


10 4 


0.6 4 


62 


2 4 


00 


5 7 


20 


16 6 


10 


0.4+4 


67 


3 4 


00 


(127) 4 


45 


1.6+ 4 


60 


Cl2) 6 


40 



Sheet 5 of fi 



PRINT SPACES BEFORE NUMBER 

+0.1 to 6 

Jump If floating point output 
d + 1 

Increase exponent for floating point output 
- d' 



n - d r if d' > 



1 if d' < 



Jump if no spaces before the number 



> Print spaces in place of non-slgnlf leant zeros 



Sign to 5 



PRINT INTEGRAL PART 



Argument to 7 
Print sign 



Print zero before decimal point if necessary 
Jump if no digits before decimal point 
d' to 4 

Print integral part of number 



Jump if no digits after decimal point 
Set up decimal point 
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9+ 




1.0 


16 6 10 
1.5 2 62 


"■ 


1 




2 2 02 
2 4 03 




1.3 4 62 
3 2 00 




2 




16 10 
1.3 2 67 




3 








1.6+ 4 63 
1.6+ 4 60 




4 




' 


► 


5 7 20 
16 6 10 


5 




1.5 4 67 
0.0 4 00 


from 8+.7 


6 






7 


10+ {lj 72 

1.0+ 60 


















10+ 




1.0 


1.3 60 




1 7 00 






l 


(T)2 40 
92 © 72 4 


to 1.0s 2 


1.0 1 00 
0.7 3 63 




from 1.0s 


17 4 10 3 
17 4 10 3 


3 




4 


072 
1.6+ 6 63 




5 
6 




( 12) 6 52 

1.7 6 11 

C 12 ) 6 53 




7 


U+ [Uj 




72 
60 



PRINT FRACTIONAL PART 

Print decimal point 

Jump if no zeros after decimal point 

- d' 



n + d' 

2 



Max (n , - d') to 2 



Print zeros after decimal point 



Jump if no significant digits in number 



Print fractional digits 



(448.0, 14) to 4 



PRINT EXPONENT AND EXIT 

LINK for return from I.O. print routine 
Exponent to 7 



2 = 3, to print 3 digits in exponent 



Enter 1.0. print routine if X3 < 
(i.e. floating point), to print exponent 



Print two spaces after exponent, but not 
in fixed point printing 

Restore accumulators 

Jump if not computing store link 

Add jump address to 1.7+ 

Restore X6 
Restore 0+ to U0 
Restore Ul 

{Exit if computing store LINK 
Jump to 0+.0+ for ordinary LINK. 
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SIGNED £.S. D. PRINT FROM PENCE 

This subroutine prints the contents of X7 as £.s.d. in the form a.b.c where 

X7 - ± (240a + 126 + c).2~ 36 

and a,b,c are integers such that 

< a 

< b < 19 

< c < 11 



On entry X2 



■38 



n,2 , where n is the required number of places of pounds. 



K> 



CK> 



to 1+. 

■* 



from 2+.6+ 





R 0-0 4 




42 - 28 - 


01 


0+072 
0.0 60 


02 


0+[0|72 
0.2 OfO 


03 


0+ 00 0. 



04 




0+ 00 0. 



LN 

SIGNED £.S.D. PRINT 

0+ 



o.o 



m 



30) 6 40 

16 6 10 

(Jo) 6 42 

0.2+ 60 



w 



6 40 

16 e io 

1+072 
0.0 60 



1 6 00 
0.6 6 63 
(l2) 6 52 
0.6 fi 11 



> 0.7 6 10 




Entry and Exit 

Print CR 
Set up LP 

Set up Sp 
Print LP or Sp 



LINK to 6 

Jump if LINK is a go order pair 

Add Computing Store link into 

N address of 0.6 

Plant LINK in 0.7 

Jump to Computing Store link 

Restore Accumulators 



Overwritten by go 
order pair link 
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from 0+.6+ 



o.o 



1+ 



[7] 73 
0.6 4 00 
(u) 5 40 

10 10 6 
0.3+ 7 62 

7 7 02 
(2l) 5 42 
0.7 24 

6 1 00 
I 2+ 7| 71 

3+072 
0.3+ 60 



448 

2 



00 



+240 









3+ 






v 


M^) 7 40 
16 7 10 






0.0 




0.2 2 60 
0.0 2 67 














1 






L 




16 5 10 
16 10 




2 






0.7 60 




L 




3 




7 12 5 




4 




16 5 10 5 
0.7 1 00 








0.6 7 61 
0. 1 5 63 






5 








> 

6 






92 [0]72 4 
0.3+ 7 62 


■^ 




L 




2+[0]72 
P-n 


to 2- 







7 














■ J 



J J 



Sheet 2 of 3 



Form a 

Preserve Accumulators 

(448.0, 2) to 4 

Set up Sp In 5 

Print <p after CR LP 

Jump If X7 is positive 

Negate if negative 



Set up - in 5 and 7 in 5j, 



X7 = 



X6 = 126 + c 



Store 126 + c in B2+. 7 



(448.0, 2) 



Print a 

Print Spaces 

Jump if n = 0; i.e. print no Spaces 

Jump if n j- 1 to print Spaces 

Print - 

Print 

Jump to B2+ through Initial Orders LINK 

Restore a to original sign 

Print Sp if a is positive 

Plant LINK for Initial Orders 

Jump to PRINT if J a I ^ 

Jump if a = and number is negative 

Enter 540.3+ if x ^ (1.0. unsigned integer 

print) 
540.7 if x < (1.0. signed integer 

print) 

LINK for Initial Orders 
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from Initial Orders , ,„, 

— , ^ ^ H 10 5 40 

and 3+.7+ 




Sheet 3 of 3 



Print b.c 

+ 10 to 5 
c 

+ 12 to 3 
c 

| 12b + c | to 7 
ft to 7; c to 6 
c to 1 
Print • 



Print ft.c with suppression of most 
significant zeros 



Overwritten by 12ft + c 
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TEXT INPUT AND OUTPUT 



01 




LN 

TEXT IN/OUT 







0+ 


foil! i 


073 
(13) 55 


v U1 J * 

0.0 


1 


I 1+ ■ 1 1 70 7 

1 6 00 














rrn 1! 




2 


1 7 00 




(J7) 1 45 
0.3+1 60 




. 


4 

, * 

5 
6 


(3) 7 53 
17 7 10 

(3) 7 53 

0.4+1 67 
35 6 01 

0.2 7 60 








1 


1+072 
0.0 60 


















1+ 






{7 

0.1 1 


72 
10 




0.0 


Interlude v 


/ fm n "7" \ 


Entry 2 


[ viy u 






\ 1 


71 / 


2 


]0 1 


70 




ll+l 


;7i 


3 


0.6 1 


00 




10 1 


71 


4 



0.5 1 00 


tp Initial 5 


[ 59 [|J72 1 
3.1 60 


Orders 


6 


1+072 
0.7 60 


from Initial w 


0+ 171 "* 


Orders 




7 


3+ H 


I 1* 
]72 



R 52 



Sheet 1 of 3 



Method 



See page 2. 



Look up index 



Output Text 

Preserve Accumulators 

Shift text number, r, to 7 m 

Read text address to 1 

Text address to 6 

Text word to 1 

Text word to 7 

Character count, n, (in l.s. bits of text word) to 1 

LOOP STOP if no text (n = 0) 

Remove n from 7 

Output one character from text word in 7 

Shift next character into 5 1. s. bits of 7 

Count n characters 

Step 6,3 

Jump if not last word 



Overwritten 
by LINKS 



Search for first L during Input 

Restore Accumulators 

Plant LINK and EXIT 

Switch to Main Tape Reader 

Plant Initial Orders Link in BO. 

Plant Assembly Link in B1+. 1 
Set Link for first L in B0. 1 



Set l s = 475 

Enter Initial Orders at 534,1 to search for 
first L 



LINK for first L 
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1 


5 


5 


5 


a 


5 



5 


5 


*\ 


m 


»7 


a 6 


a 5 


a 4 


: a 3 


°2 


a i 


n 



No. Bits 
Contents 



n Is the number of characters to be printed from the word 
o; (i = 1,2,.... 7) is the i* character to be printed from the word 
m = if there are further words in the text 
= 1 if the word is the last word in the text 







2+ 


1.0 

1 


5.7 3 00 

4+ [3] 72 

5+ [|] 72 

0.1 1 00 




|o 1 


71 


from 5+.l+^ 


2 


17 4 


no 


y 










1.2+4 60 
31 4 42 4 




3 




1.2+4 60 

.Qai 1 — 

1.2+ 4 60 
(vty 4 41 




4 




5 




1.2+4 60 
(l5) 4 43 
1.7 4 61 

16 4 00 






e 

7 


















K 




2.0 


(15) 4 43 

1.7+4 60 

16 2 00 

(l5) 2 43 












, ^ 

1 

2 




2.1 2 60 
(j) 2 43 




3 


2.3 2 61 
(j) 4 41 








2.7+4 60 
(14) 4 43 






4 


to 5+.5+ 


4.5+4 60 

(Uj) 4 41 

2.6 4 61 

77 






5 

6 

7 
^ ► 




[l.2+ 60 




3 1 00 



Search for Letter Shift 

Set text address in 3,,, 



Restore Assembly Link to B0. 1 

Read character negatively, -c lf to 4^ 

Ignore <fi 

c 1 - 31 to 4 C 

Ignore Er 



c 1 - 30 



Ignore CR 
c 1 - 13 
Ignore LP 

c t - 27 

LOOP STOP if not k 



Read warning character from tape, c t , to 4 C 



Search for L, N or Z 

c 2 - 15 

Ignore Er 

Read next character from tape c 3 to 2 C 

c 3 - 15 

Ignore Er 

c 3 - 16 

LOOP STOP if not $ 

c 2 - 14 

Jump if warning character is N 



c 2 - 28 








Jump if L 








c 2 - 10 








LOOP STOP if not 


N, 


L, or 


Z 


77-stop for Z 








Return to search 


for 


next 


\ 


Text address to 


I 







PE 
PEGASUS 



R 52 




to 2+.2+ 



2.7 2 00 
(0)4 40 4 
17 4 OX 
4 6 00 



0, 



55 
3.2+ 5 67 



3.6+4 60 2 
(V)7 46 
7 1 00 
|0 0| 71 3 



4.0 



(2 



from 3+.5 



5 



5+ 



35 3 01 
3.1+1 62 
5.7 3 10 
1.2+0 60 



~0 : 



5 6 00 


V |7 53 
4.3 5 67 

6 7 06 
32 7 06 

3.6+ 60 



(T) 74 

10 01 70 

4.7 1 10 



r 
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Sheet 3 of 3 

Read and Store Text 

Write text address in Index 

Step index address 

Clear 7 

Set character count, n, In 5 

/ 3,4 before first character 
Jump to ^ 4>2 after ^ (last charac ter) 

6 to the last 7 bits of 2 ffl 

Last character to 4 C , clear 4^, 

Add next character from tape, c 3 , into 4^ 
Add last character from 4 C to 7 
Build up text word in 7 



Count a maximum of 7 characters 
Jump if last text word, to mark 
Set n = number of characters in word (Normally 7) 

Write text word to store 

End of Text 

Step text address by 0.1 
Jump if more words of text to be read 
Reset Transfer Address in 5.7 m 
Return to search for next text 
Character count = 7-n to X6 

Shift text down to correct position 

Set 7-« at end of X7 
Mark last text word 

Switch to 2nd Tape Reader 
Plant Initial Orders Link 

Overwritten by Link 



T 2+.1 

J 560.0+ - l+.l 

L 

Perranti Ltd., London Computer Centre, 21, Portland Place, LONDON, W. 1. 
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DOUBLE LENGTH INPUT 



A double-length fraction may be regarded as being of the form ±0.x, where x 
represents a string of 23 decimal digits. If there are fewer than 23 digits, the 
remainder are treated as zeros; if more than 23 digits, they are ignored. 

The sign of the fraction is stored and five is subtracted from its first decimal 
digit on input. The 23 digits are then read to form an integer, I, in the range 

~5. 10" < I < +5. 10" 

This integer may be converted to the required fraction, F, since 

F = 10" 23 x 2 7 * (5.10 22 .2~ 76 + I) 



= i + 10" 23 x 2 76 x I 



It would be possible to store 2 7 * x 10~ 23 and use the above relation to form 
the required fraction. If this were done double-length there would be a loss of 
accuracy and the equation has therefore been re-arranged as follows: 



F = | + /-7(1 



10~ 23 .2 7 *) 



~ + i-ij{4a-io-".2 76 )} 



The double length fraction 



A + 5.2" 3B = 4(1 



10' 



'".2 7 «) 



is stored and this last relation is used to form F. 



R 


0-0 3 


100 


- 28 - 


0+ 


©72 


0.0 


60 


o+ 


00 0. 












0+ 


00 0. 



LN 

D.L. INPUT A 
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(0+) 



CK> 



Read next no. 



Prom 4+.6+ 
(sheet 6) 



(1+) 



0.0 



[7J73 
5+ [D 73 



1+ [TJ72 
35 4 00 



6 00 
7 00 



2+ [2] 72 
1.1 60 



1.0 






16 2 01 
1.0+ 2 61 



1.1 



16 5 00 

lo) 5 43 

Qjl.2 5 63 

1.4 5 60 



0.4 5 67 



G 



1.3+ .0 4 00 



(1+) 



1.4 



16 2 00 
12) 2 43 



from 7+.1+ 



(sheet 4) 



(2+) 



23) 3 40 

2.1 2 61 

35 5 02 

(T) 7 42 



16 2 00 
15) 2 43 



2.0 



1.7 2 60 
(T) 2 41 

cz 



ENTHV 

Preserve Accumulators 
Preserve U2 in B5+ 

Sign marker set for + 
Clear 6 and 7 



SEARCH FOR SIGN 

Add character after CR to -13 
LOOP STOP if CR not followed by LP 



^-10 

LOOP STOP If decimal digit before sign 



11 



(0+) 



Mark 
minus 



(T) 5 43 
("A 0.4+ 5 63 



LF Sp Er <p 



CR 



1 



13) 2 42 
(T) 5 43 

1.1 5 63 
(l3) 5 43 

1.0 5 60 

6+ \T\12 



C t - 13 

LOOP STOP if • 
before sign 

-13 to 2 

c i -» 

Ignore if. Sp, Er, <j> 

c, - 30 



12 



Jump if CR 



-^.To 6+.0 to search 
for k L <p 
(sheet 7) 



Set counter for decimal digits In 3 
Jump if not • 
Mark fraction 
-5 to 7 

Not 
Ignore Er 



-^•To 2+. 1 to Read Number 
(sheet 3) 
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from 1+.5 + If Integer, or 2+.0+ if Fraction 

BEAD VISITS OF NOSffiEE 
(2+) 



2.1 



V 




3+ Q] 72 




(?) 2 41 




9+ [0]72 




6+ o ro 











c 3 -10 



10} 2 41 

2 7 01 

6 23 

1.0 3 67 

1.0 60 



Jump if not decimal digit - 

Restore decimal digit 
Add new digit to number 
Just Ify 
Count digits - 



from 9+. 7+ after Er 
(sheet 4) 

from 7+.7+ with OVB 
after ±1.0 (sheet 4) 



from 9+,6 (sheet 4) 



with X3<0 
End of 
fraction 



(3+) 




5.0 6 10 
1 7 10 




1.2 




3 


(?) 54 
5.0 6 01 




4 


1 7 01 
6 23 




S 


(?) 54 
0.O+ 65 


I 
1 
i 

2.3 3 62 

I 
i 
i 




6 


1.2 3 67 ' 







(W 



mpq) 




TO B+.6+ 

^ (sheet 4) 

if character not 
decimal digit. 



Jump if not decimal digit 

LOOP STOP if more than 23 
digits read in integer. If I 
fraction, ignore digits 
after 23 have been read. 



Jump if OVR set to 9+.0+ (sheet 4). Reaches LOOP STOP in 
7+.2+ unless ±1.0 



Jump to read next digit 



Count number of multiplications by 10 




Jump when I formed (see sheet 1) 



To 6+. 4+ to Convert Fraction 
(sheet 5) 
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(9+) 



Sp 



READ NON-DIGITAL CHARACTERS 



Not 
CR 



0.6+ 


(T) 2 43 


c„ - 14 
Ignore Er 






0.4+ 2 60 
0.2 2 67 






0.7 






















c„ - 30 

Jump If not CR 




^ 


(l5) 2 43 
0.0+ 2 61 






0.2 








1 

CR 








16 2 00 


C E 




0.3 

OS 

4 


(l3) 2 43 


e t - 13 

STOP If CR not foil 




0.4 2 61 
2.5+ 5 60 


owed by LP 











To 3+.0 
(sheet 3) 



T 
Fraction 

__Jf__ 

0.5 1.7 3 60 



*fcto 2+.5+ if Integer 
(sheet 6) 



0.5 + 
0.6 



< 23 

. * 

35 3 03 

1.2 60 



0.6+ 



(7+) 



P > 1.0 



— *> 



3+ (I] 72 
7+ [0]72 
0.2 60 



18) 2 41 
0.3 V 
<^A 0.2+ 2 61 



0.3 0.1 7 60 [ 

.::r:~ 

0.3+ 0.2 7 67 



±1. 



0.4 



16 2 00 

(id) 2 43 

0.4 2 60 

(T5) 2 41 

0.4 2 60 

32 02 

9+072 

1.0+ 60 



TeBt number of digits read 



> 23 



(-0.1, 23-c) to 3 



to 3+.2 to multiply 
->by 10 23 ~ c to form 
I - Number x 10 23 



(3+) 



1.7 



6+ Q] 72 
1.4+ 60 



T 



-Overflow 
from 3+,5+ (sheet 3) 
'or 9+.0 (sheet 6) with 2 C < 18 
but see note on sheet 6 



STOP on OVR or inadmissible character 



Jump if t preceded by 



to 6+.4+ to Convert Fraction 
(sheet 5) 



Jump to STOP if F>1.0 



(7+) 



0.1 



1 + H^ 2 
1.5 60 



to 1+.5 to Read Fraction 
(sheet 2) 



c 6 - 15 
Ignore Erase 

Ignore 

Set OVR, to indicate ±1.0 



,to 3+. 0+ (sheet 3) to test for terminating 
"character. Otherwise LOOP STOP on OVR 
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(4+) 



1.7 



(7+) 



0.0 



7.5 


1 11 4. 


95 


7 55 2 


54 


5 65 6. 


98 


1 22 2 



from 3+.7+ (sheet 3) 



(6+) 



(8+) 





y 


1.4 + 


l.o e 


> 10 




l.i ' 


r 10 


5 




) 70 




7+ ( 








1.1 1 


20 


6 








8+ [2] 72 




6 r 


r 00 


7 








t 


i 00 




> 


( 




1.0 1 


22 


2.0 




r| 70 




|4+ r 




1.1 1 


22 


1 








6 •" 


r 00 




€ 


00 


2 








6 C 


) 23 




1.0 1 


22 


3 


©' 


' 41 




6 C 


1 23 


i 


r Q. c 


55 




1.0 € 


04 


5 








1.1 7 


04 




6 C 


23 


6 








33 € 


01 




2+ [2 


]72 


7 








2.5+ C 


1 60 



J 
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most significant half of double length 
number for multiplication of fractions 



least significant half of double length 
number for multiplication of fractions 



CONVERT FRACTION 



Store double length number in 1.0 and 1.1 



Read B to 1 



> S.<?.2~ 3B to PQ 



B(p + 9.2 -38 ) to PQ 
A to 1 

[Aq + B(p + ?.2" 38 )].2" 3B to PQ 

(A + 5.2~ 39 ) (p + q.2~ 36 ) to PQ 
> \ (A + 5.2" 3B ) {p + q.2" 38 ), rounded up, to PQ 



(P + ?.2- 38 ) -j(A+ B.2' aB ) (P + q.2' 36 ) + 1 
= {p + q.2~ 3 *).2 7i .10~ a3 +-j to PQ 



To 2+.5+ to Store Number 
(sheet 6) 
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(2+) 



STORE NUMBER 

.Prom 8+. 7+ (sheet 5) if fraction 





i 


Prom 9+.4+ (sheet 4) if integer 


2.5 + 
2,6 


4+ [T|72 
0.0 65 


Jump if ±1.0 









2.6+ 


1.2+ 4 61 


2.7 


6 6 02 
1.1 60 



(4+) 



1.2+ 


5.7 


5 00 




1.5 


4 00 


3 











] 


00 4 


|o 


0|71 5 


4 








35 


5 01 


. 


1.3+ 


4 66 




5 








5.7 


5 10 




0+ 


072 


6 








0.1 


60 



(9+) 
±1.0 0.0 



Jump if + 
Negate 6 



(4+) 



1.0 



1.0 7 60 



I|F|>1.0 

To 9+.0+ 
but see foot of page 



) 7 40 4 



1.0+ 



32 6 00 

7 7 02 

6 23 

1.2+ 65 



7. = if F = +1.0 
=1 if F = -1.0 



-1.0 to 6 

Negate 7 

Justify 

Clear OVR if F=+1.0 



Transfer Address (T.A. ) to 5, 
Set 4j, = 6 i.e. counter of 2 



m 



t Store double- length number in two consecutive 
positions in the Main Store 



(T.A. +0.2) to 5.7- 



->.To 0+. 1 to Read Next Number 
(sheet 2) 



Note: The order in 9+. is intended to detect \f\ > 1.0 but this should previously 
have been detected in 3+. 5+ 
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from 0+. 7+ (sheet 2) 
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SEARCH FOR WARNING CHARACTER L 



Cj - 27 



LOOP STOP if not k 



C 5 -28 



(c s , c 5 - 28) 

LOOP STOP if k not followed by L <p 

Restore Accumulators 

Restore U2 

Plant LINK 

Does no harm, used in conversion 

Overwritten by LINK. Obey LINK 



London Computer Centre, 
21, Portland Place, 
LONDON, W. 1. 



© FERRANTI LID 1958 

got to be reproduced in whole or 

in part without the prior written 

permission of Hrranti Ltd. 



Issue 2 

23rd October, 1958. 

C.R.M. 



FERRANTI LTD 
PEGASUS LIBRARY PROGRAMME 



R 101 



FLOATING POINT INPUT 





R 0-0 2 




101 - 28 - 


01 


0+[3]72 
0.0 60 


02 


0+10)72 

0.0+ 60 



LN 

P.P. READ 



R 8 


2 -0 1 


101 


- 02 - 


R 8 


3 -0 1 


101 


- 02 - 


R 8 


3 -0 1 


101 


- 02 - 



CoT> 



o.o 



roT> 



from 



/9+. l+\ 

\ 5+. 1+J 



from 2+. 1+ 



from 1+.2 



to 1+.0 



to 1+.2+ 



0+ 



5 00 

| 0| 71 

2 00 

20 

1+[T]72 

16 1 00 

(n) 1 43 

1.0 1 62 

1 41 

0.2+ 1 62 



© 

0.2+ 
0141 
0.7 1 60 



tp 3+.0 
or 4+.0 



1.2+ 1 67 
32 1 00 

3+072 7 
0.0 60 



Sheet 1 of 6 



Calls for P. P. 01 



Search for sign 

Mark to read only one number 
Store LINK 
Clear 2 
Clear 6 and 7 



Character from tape, c , to 1 



Ignore <p, Sp, Er, LF 



c - 10 
l 

Jump if + 

Jump if not - 

Mark - 

To 3+.0 to read number 

or 4+. to read exponent (after Q) 
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from 0+.3+ 




Cy i 43 



to 5+. 



0.6+ 1 61 
34 7 00 



0.1+ 60 



after N 



L_». Qf) 1 41 
Qi| 0.7 1 61 

5+072 



Search for Warning Character 



c - 30 



Jump if not CR 



Character from tape, c , to 1 



c - 13 

2 

Jump If LF 

LOOP STOP for punching error. 



c - 27 



Jump if not X 



Warning character from tape, c , to 1 



Ignore Er 



Search for after warning character 



Ignore Er 



Interpret Warning Character 



■ After warning character L 



Over- 
■ written 
by LINK 



Read LINK 

Plant LINK 

77 stop after warning character Z 

Return to Number Read 

Jump to LOOP STOP if not Rafter warning character 

c 3 -28 

Jump if L 

c - 10 

3 

Jump if Z 
c - 17 

3 

Jump if not Q 
Set marker in 7 
Jump if Q 



c - 14 

3 



LOOP STOP if not L, Z, Q or N 
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Read Number 



Sheet 3 of 6 



from 2+.7+ 




(~2) 3 43 



1.7 3 60 
16) 3 43 



►Q* 1-5 3 61 



to 7+.0 



16 3 00 

(13) 3 43 

1.5 3 61 

7+[T]72 

1.0 60 



Clear 4 



Clear 0.0 



Used to store argument 



Used to store exponent 



Block exponent to 0.1 



Character from tape, c , to 3 

c to 7 

3 

c - 10 

3 

Jump if not a decimal digit 

+ 10 to 3 c 

Multiply by 10 and add next digit 

Jump if number does not 0VR 

Jump to ignore overflow digits after • 

Prepare to increase exponent if overflow 
digits before « 

Store number 

Reduce exponent by 1 for each digit after 



c - 15 

3 

Ignore Er 

c - 12 
3 

Jump if not decimal point 
Mark decimal point in 4 



S -14 



Jump if Sp 



c - 30 

3 



LOOP STOP if not CR LP. 



Character from tape, c , to 3 

c - 13 

Jump to LOOP STOP if not LF after CR 
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from 0+. 7+ 



o. 



4+ 



2 00 
16 3 00 

3 7 00 
(w) 3 43 
0.4+ 3 62 
(To) 3 40 

2 3 22 

7 2 00 
0.0+ 60 



U 



from 2+. 7+ 



3 43 
0.0+ 3 60 

.Q.H 1 - 

5+[Tj72 
1.6 3 60 
(7e) 3 43 
1.4 60 



1.0 



^ 



JZI 



3 40 3 

17 3 10 

17 3 01 

1.0 3 60 1 

1.0 1 61 

34 1 00 

35 1 03 
1.0 60 



to 0+. 1+ 
•^ 



->^C3 1-4 3 61 

16 3 00 

(T3) 3 43 

1.4 3 61 

1.7 1 62 

2 2 02 

0+[0]72 

0.1+ 60 
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Read Decimal Exponent 

Clear 2 

Character from tape, c , to 3 

5 

c to 7 

5 

c - 10 
5 

Jump if not a decimal digit 

+ 10 to 3 c 

Multiply by 10 and add next digit 

Copy exponent to 2 



c - 15 
Ignore Er 



c - 14 



Jump if Sp 



c - 30 



Print Name; Set Exponent 

Last character to 3 ; clear 3 m 
Print 4> on entry, then last character. 
Next character to 3 ffi , leave 3 c 
Jump If last 2 characters were <p 4> 
Jump except at first entry 

Set 7 In 1„ 



LOOP STOP if not CR or CR LF 
Character from tape, c , to 3 
c - 13 

6 

Jump if not LP after CR 
Jump if marker is positive 
Negate block exponent 



FER R AN T I LTD 
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Convert decimal exponent in 0.1 to binary exponent in X3 



from 6+.7+ 



1.0 



7+ 



0.0 6 00 
0.1 4 00 
8+072 
0.7 6 60 
(38) 3 40 
(38) 3 56 
0.4 4 60 
1.6 4 63 
(T) 4 43 

© 



3 41 
0.0 6 20 

1.2+ 60 




© 



4 41 
3 43 
0.1 6 20 
1.2+ 60 



o.O 



8+ 



+ .625 



+.8 



+0 



+0 



* 0.5+1 62 



to 9+. 3+ 



6 6 02 

3 56 

0.2 6 03 

3 56 

3 7 00 

9+[T]72 

1.3+ 6 60 





0.2 

© 



Argument to 6 
Exponent to 4 

Jump if argument = 

Set binary exponent 

NORMALIZE 

Jump if decimal exponent = 

Jump if decimal exponent < 

Subtract 1 from decimal exponent 

Add 4 to binary exponent 

Multiply argument by (10.2 ) 

Jump 

Add 1 to decimal exponent 

Subtract 3 from binary exponent 

Multiply argument by (2 .10" ) 

Jump 



10.2 



2 3 .10 _1 



+ P. P. 01 = -2 

+ 2(P.P.01) = -2" 

Jump if marker is positive 

Negate argument 

Re -normalize 

Add in round-off constant 

Re-normalize 

Binary exponent to 7 

Jump if argument is zero 



R 101 
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„n-i 
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Pack number, store and count 



from 8+. 7+ 



argument ^ 



1.0 



from 8+.7+ 



argument = 



to 0+. 1+ 

-K 



9+ 



±i 1.1 7 62 

7 7 02 

0.2 7 01 

0.3 6 05 

0.2 3 03 

3 6 01 

1.6 5 60 
6 1 00 

|0 0| 71 5 
35 5 01 

1.7 5 67 
|~b~~0| 70 

1.7 1 10 
■ / 0+072 \ 
60 J 



,0.1+ 



01 



R 





-0 


1 


101 


- 


06 


- 


_ . 










256 
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Jump If binary exponent > 



Negate binary exponent 



Subtract 2"" 1 from X7 



LOOP STOP if C(7) ^ 0, i.e. exponent overflows 



Clear last n digits of 6 



Add 2 to exponent 



Pack 

Jump if C(5) = 0, i.e. Cue 01 

Write number into Main Store 

Step 5 

m 

Count numbers 
Plant LINK 



Jump to read next number 



Optional Parameter List 



Title 



P.P. 01 = -2 (i.e. n = 9) 



Overwritten 
by LINK 
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R 102 



READ IN TABLES OF NUMBERS 



C°L> 



from I.Os 



Input 



GD- 



01 



02 



R 





-0 2 


102 


- 


28 - 


0+ 





72 


0.0 





60 


0+ 


@ 


72 


0.5 





60 



LN 

INPUT TABLES 



R 2-0 1 
102 - 02 - 



o.o 



0+ 



0.2 7 00 

[4] 72 7 
4.0 7 01 7 
0.6+ 60 



+0 



0+072 
0.7 7 67 



6 1 00 
0.5 1 10 



f 0.2 7 00 



^07 41 



35 7 03 
1 6 00 



35 7 01 
1+072 



Sheet 1 of 2 



Call for P. P. 01 



7„ = B.P - Address of 1st word of Index 



Block B to U4 
(B.P, t) to 7 

+P.P.01 = (B.P, 0) 



Count Tables 



LINK for Initial Orders, 
obeyed in 0.3 



Plant LINK and exit 



(B.P-1, 1) to 7 

Store LINK in 6 
Step address 



• Overwritten by LINK 
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1.0 



2 







1+ 





[|] 72 7 


5.7 


1 00 


35 


1 03 


33 


5 00 


1.2 


1 66 


5.7 


1 10 


4.0 


1 11 7 





@73 7 


0.3 


1 00 







10 


072 5. 


0.0 


60 


R 


-0 1 


102 


- 06 - 


. 


- 00 0. 








Sheet 2 of 2 



Block to contain Index to U4 



Bring Transfer Address to start of next 
available block and set X5 = (512.0, 0) 



Store Table Address in Index 
Write up Index into correct block 
Set LINK for Initial Orders 

Optional Stop before Reading Table 
Enter Initial Orders Input 
Optional parameter list 

P.P. 01 
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DOUBLE LENGTH INPUT 
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01 



R 





-0 1 


105 


- 


28 - 


0+ 





72 


0.0 





60 



LN 

DL INPUT 







0+ 


C ni "N > 


1+0 72 


C 01 ) " 


0.0 






l 


lb 7 UU 

C 10) 7 43 




OS 


0.1+ 7 63 




(7)7 43 




2 








0.7+ 7 63 




e 


0.3 7 60 




3 


(T) 7 43 




- 




4 


0.0+ 7 63 
CT3~) 7 43 




0.0+ 7 60 




5 






6 


0.5+ 60 




+0 









7 






I ». 


6 00 



Search for start of number 



1st character, c 1§ from tape to 7 



Cj - 10 



LOOP STOP If a decimal digit 



c l - 


12 










Jump 


if + 


or 


- 






LOOP 


STOP 


if 


decimal 


point 


c i " 


17 










Ignore Er, 


Sp 


. LP 


or 


<p bef 


c i ~ 


30 











Ignore CR before sign 

LOOP STOP for punching error 

Unused 



Clear 6 
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l.o 



(T)7 43 

1.1 7 60 

2+072 

0.0 60 



Read integral part of number to 6 



c - 10 



Read next character, c 



c - 10 

2 



Jump if not decimal digit 



Multiply by 10 and add 
new digit 



LOOP STOP if OVR 



c - 15 

2 



Ignore Er 



Overwritten by 
sign marker 



Overwritten by 
integral part 



Overwritten by 
LINK 



Overwritten by 
fractional part 



Overwritten by 
+0.1 



to 4+. 



3+.0 



2+ 



0.0 



+0 



+0.1 



Prepare to read fractional part 

Integral part to 1.1 
Plant LINK 
Clear 1.3 



c a -12 

Jump if decimal point 



+0.1 to 1 and 1.4 



Unused 
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3+ 


from 2+. 5+ ^ 


16 7 00 
J 10 ) 7 43 




0.0 




0.4+ 7 62 
("To") 7 41 




I 
2 
3 




1 7 20 
1.3 7 11 




1.4 1 21 
1 6 10 


^ 


0.0 60 




4 




(?) 7 43 








0.0 7 60 

( 15) 7 43 

4+[o]72 

0.0+ 60 


5 

6 




7 




+0 














4+ 


from 2+.3+ 




( 18) 7 43 
0.0+ 7 61 




0.0 


l 


16 7 00 
C 13) 7 43 


OS 

2 


0.2 7 61 
1.0 1 00 






0.5+ 1 60 
1.1 6 02 




3 
4 
S 




1.3 7 02 
6 23 




1.2 60 




1.1 6 00 






6 


1.3 7 00 

1.2 60 



Read fractional part of number to 1.3 

Read next character, c 3 
c 3 -10 

Jump if not decimal digit 

Multiply digit by (0.1)" 
Accumulate number in 1.3 



(0.1) n+1 to 1 



c 3 -15 
Ignore Er 



c 3 -30 



Unused 



Complete number and exit 

c - 30 

2 

LOOP STOP if not CR after last decimal digit 
Read character after CR, c 



c n - 13 



LOOP STOP if not LF after CR 

Sign marker to 1 

Jump if + 

Negate Integral part in 6 

Negate Fractional part in 7 

Justify 

Jump to obey LINK 

Restore Integral part to 6 

Restore Fractional part to 7 

Jump to obey LINK 
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SHORT NUMBER READ 



This subroutine reads a single length mixed number, fil, from tape and leaves it 
in the form 

p' = tf.Kf.2* 38 q' ~ l<f.2" 38 

where m is the number of figures after the decimal point. 



from 2+.0+ 



to 1+ 



from 2+.7 



C°E> 



to 1+.4+ 
-< 





R 0-0 3 




116 - 28 - 


01 


0+072 
0.4 60 


02 


0+ 00 0. 



03 




0+ 00 0. 



LN 

SHORT NUMBER READ 



, r 



0.0 



to 2+.7+ < i. 

-« 



from 1+. 5+ 



0+ 



0.1 5 67 
0.6+ 7 63 



1.4 5 67 

0.6+ 7 63 



5 6 12 1 
0.3+ 7 61 




6 00 
1 00 



1.4+ 60 



(Jo\ 5 41 



Jump if not LP 

Jump if in a or /3- search 

Jump if not Sp 

Jump if in a or /3-search 

Negate 6 if tl negative 

Jump if m t 

Put 7 = 1.2" 38 f or m = 

Jump to LOOP STOP if OVR set 

Overwritten by LINK 



Plant 


LINK in 


0.4 




-1.0 


to 7 






Clear 


6 






Clear 


1 






Jump 


to read 


character 


from tape 


c x to 


5 
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Note that q'( = 2A0 r ) will not overflow into p' on the order 1 54 in 
1+. unless OVR has already been set by the order 2 5 50 in 1+.2+. 

p ! - 2p and q' = 2q. Clear sign of 7 
Form AMCP.2 -38 in 6 







1+ 




1.0 


C^)o 54 






6 5 01 


V 


1 


(T)6 50 






5 6 01 






7 5 00 


2 


(T)5 50 




3 


5 7 01 
1.4+ 1 62 


4 


f 




from 0+,l ^ 


2.3+ 5 67 
16 5 00 










from ' ^ 


4 




0+ 7 




r ^ 




to 0+. 


7+ 


5 


(p)\ 5 43 
0.7+ 5 63 


















1.7+ 5 60 










<- 




2.0+ 5 67 






35 1 00 












7 














1 > 


32 7 06 


















2+ 








2.2+ 7 62 
0.0 5 67 




to 0+. 





1 


2.0 

^ — -> 




•^k— 
















2.2 7 63 












^ 


2.1+ 7 61 










L 


^ > 

2 


M/)7 40 

Ti 1 (V\ 








w' 








1.4+ 60 












._ _ i. 


2.5 5 67 
















32 1 00 








4 


1.4+ 7 63 








^ > 






U4J 5 43 




C^ 


2.5+ 5 61 






16 5 00 




6 

t£> 0+.1 + 


(T3) 5 43 
n 1+ 5 fin 




from C 


+. 


3+ 




-J 


jr\ 


2.7+ 60 





Form ltf".2~ 38 in 7 

Jump unless digit read in a-search 

Jump if not Erase 

Read character from tape, c , to 5 

c^lO 

Jump if decimal digit 

Jump if + 

Jump if not - 

Set lp - 1 to mark negative number 

Change sign digit of 7 for + or - 



Jump unless sign read in middle of number 

Jump if not • 

Jump if in a or /^-search 

LOOP STOP if two • in one number 



-38 



to 7 after • 



Retain lp ( = 1 if - read, otherwise 0) 

Jump to read next character 

Jump if not 4> 

-1.0 to 1 to mark a-search 

Jump unless <p read in middle of number 

Cj - 30 

LOOP STOP on inadmissible character °>>+ ? b 

Read character from tape, c 2 , to 5 A ^***^ 

c 2 - 13 

Jump if CR followed by LF 

LOOP STOP 

Ferranti Ltd., London Computer Centre, 21, Portland Place, LONDON, W.l. 
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GENERAL SINGLE-LENGTH NUMBER READ 

This subroutine reads a single-length number, ff, from tape:- 

p' = ^. ltf". 2~ as and q r = m.2" 38 where m is the number of figures punched 
after the decimal point. 



R 121 



Sheet 1 of 3 



01 



02 



03 



R 





-0 


3 


121 


- 


28 


- 


0+ 


m 


72 




0.5+ 





60 




0+ 





00 


0. 




















0+ 





00 


0. 



LN 

NUMBER READ A 









0+ 


from 3+.6+ 




6 7 12 3 




0.0 






1.0 7 10 






1.1 4 10 




1 


072 
1.0 6 00 




2 


1.1 7 00 






0.4+ 1 63 
(jJ6) 1 50 




3 


<. 


60 1 




4 


1.7 1 10 








1.7 60 




5 


CUP) 




073 


\JiLJ 








1+072 




6 


1.3 5 00 




7. 5 10 5 




7 






7 R •*** 











Exit and Entry 

Negate 7 if 3^, = 1 

jrarf".2""'to i.o 



m.2* 3e to 1.1 



Restore Accumulators 
t!.ld n .2~ 3e to 6 



m.2" 38 to 7 



Jump if LINK is go order pair 

EXIT by Computing Store Link 

Plant LINK 
Jump to LINK 
Preserve Accumulators 

Set 0.4 in S^ 
• Clear 1, 2, 3 and 4 
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Search for First Character of Number 







1+ 


1.0 


4+072 
5 00 


' 






16 7 00 




1 


(To) 7 43 




2 

J 


1.4 7 62 
1.2+ 5 61 




@7 41 




^- » 


1.6+ 60 




1.6 7 60 




4 


(T) 7 43 




5 


0.0 7 61 
35 3 00 




• *■ 


32 2 00 




r *" 


2+[0]72 




3+ [T] 72 


to 2+.3+ 


7 


0.3+ 60 
















Ul 




w, 


(T) 7 43 




0,0 




/■ 


0.2 7 61 




35 1 00 




1 

. >. 

3 


1.6+ 60 




(~Z\ 7 43 






J Qt7«5 












1 1 


7 Rfl 


, 


3 

4 


C^)l 43 
fT) 5 40 
1.1 7 60 




S 




(l4) 7 43 




- 


0.7+ 7 61 




16 7 00 







(13) 7 43 




- 1 0+ 7 60 




1 
[ *<2 


* 0.7+0 60 

, 



Set /3 - search Indicator 

Read character from tape, c 1 , to 7 

c 1 - 10 

Jump If not a decimal digit 

LOOP STOP on decimal digit in a - search 

Decimal Digit in $ - search 

Jump to MAIN INPUT with C 1 In 7 

Jump if + 



c l-» 




Jump if not - 




Set 3 m = 1 if - 


read 



Set -1.0 in 2 if + or - read 



TO MAIN INPUT 



Cj - 12 

Jump if not # 

Set ^ = 1 ifo read 

Jump to MAIN INPUT 

c 1 - 15 (Cj > 13) 

Return to j3 - search if LP or Sp 

Ignore Er in a or /3 - search 

c x - 16 

Set a - search indicator 

Jump to a - search on 4> 

c t - 30 

Jump on inadmissible character 

Read character from tape, c 2 , to 7 

c z - 13 

Return to /3 - search on CR LP 

LOOP STOP on inadmissible characters 
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R 121 



Main Input 



from 1+. 7+ 



o.o 







2 40 
0.3+ 60 



.1.0 



to 0+.0 



&' 



+10 



Sheet 3 of 3 



c } to 5 




Multiply number by 10 




Add new digit c 3 




Justify 




LOOP STOP if 6 ? 




Add 1 to 4 for each digit 


after :• 


Clear 2 if last character not Sp or Er 


Read character from tape, 


c 3 , to 5 


c 3 - 10 




Jump if decimal digit 




C 3 -14 





Jump if not Space 
Ignore Spaces after sign 
Jump to EXIT on 2nd internal Sp. 
Set 1 in 2 on 1st internal Sp. 
Jump to read next character 



c 3 - 12 




Jump if not > 


LOOP STOP if 


2nd :• 


Set Lj, = 1 if > read 


Jump to read 


next character 


c 3 - 15 




Ignore Er 




c 3 -30 





Jump to LOOP STOP on inadmissible character 

Read character from tape, c„, to 5 

c 4 - 13 

LOOP STOP on inadmissible character 

Jump to EXIT 
Overwritten by LINK 



Perranti 
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SIGNED £.S.D. READ TO PENCE 



This subroutine reads from tape a sterling sum a.b.c and leaves it in XI as 
{240a + 12i> + c). The sum M must not exceed £1, 145,324,612. 5.3d in modulus, 
i.e. I 240a + ub + c | < 2 3B 



01 



02 



03 



R 





-0 3 


142 


- 


28 - 


0+ 





72 


0.5+ 





60 


0+ 





00 0. 
















0+ 





00 0. 



LN 

SIGNED £.S.D. READ 









0+ 


from 5+. 6+ 




6 7 12 4 




0.0 


1.7 7 10 

072 




l 


0.3+ 6 63 








6 1 00 




2 


(u) 1 52 




0.4 1 11 




3 








1.7 1 00 

1.7 6 10 








4 






5 


6 62 




1.7 60 




o \jjn 


C 01 "> 


^ U1 J 






1+[T] 72 




6 


2 00 
3 00 




7 


4 00 



Exit and Entry 

Negate X7 if number is negative 

Store number in 1.7 

Restore Accumulators 

Jump if go order pair LINK 



Add Computing Store link into 
N address of 0.4+ 



Number to 1 

Plant LINK in 1.7 

Jump if Computing Store Link 

Jump to LINK (go order pair) 

Preserve Accumulators 



Clear X2, X3 and X4 
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Search for first significant character 





Fl 




1.0 


2+072 
5 00 




V 






16 7 00 
(To) 7 43 

1.4+ 7 62 
1.2+ 5 61 






1 




> 


2 
3 

4 

' 

5 

6 

i >■ 

7 




(lO) 7 41 
(?) 2 40 
1.6+ 60 




1.6+ 7 60 
(?) 7 43 
0.0 7 61 


ip 3+.5 


35 4 00 

3+072 

4+[T|72 

0.5 60 










m 




> 

0.0 

1 


(?) 7 43 

0.2 7 61 

32 3 00 

1.6+ 60 




(*?) 7 43 
1.0+ 7 63 
1.1 7 60 

(?) 7 43 

(?) 5 40 

1.1 7 60 






2 








3 
4 










6 


(h) 7 43 
0.7+ 7 61 






6 




16 7 00 
(13) 7 43 
1.0+ 7 60 
0.7+ 60 




7 



Set j3- search indicator 

Read character from tape, c , to 7 

c i - 10 

Jump if not a decimal digit 

LOOP STOP on decimal digit in a-search 



+ 1 to 2 C 

Jump if + 

Jump if not - 

Marker in 4 for negative number 

To read number 



c i - 12 



Jump if not • 

Set X3 = -1.0 after first point 

c i - 15 

Jump to f3- search if LP or Sp 

Ignore Er 

c - 16 

Set a- search indicator 

Jump to cl~ search if 



c i - 30 

Jump to LOOP STOP for punching error 



Read character from tape, c , to 7 
c„ - 13 

Jump to /3- search after CR LF 
LOOP STOP for punching error 
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Read the digits of the Number 







3+ 






(T) 2 40 








1 


(To) 5 41 
(To) 6 40 






6 7 20 










5 7 01 






2 








Qi 


6 23 




0.3 6 61 










0.5 2 61 










1.7 7 10 








4 










*- 


7 00 


from 1+.7+ 


16 5 00 
(To) 5 43 


and5+.7+ 






5 










5 63 






6 


fT) 5 43 
1 4+ 5 61 








7 


3 1 00 












4+ 






2 00 








1.0 










1 


32 3 00 




(20) 6 40 










0.1+ 1 60 










35 3 00 








2 


(n) 6 40 

1.7 7 01 








3 












0.1+ 1 63 






4 




1.4 60 




(~2) 5 43 














1.6 5 61 












5 




v 




0.5 2 60 






5+[0J72 






jo 5+.0 


0.0 60 




7 


(' - ) 












\ / 



+ 1 to 2 after decimal digit 



+ 10 to 6„ 



N' = x ,N 



(x = radix) 



Add new digit to M 



LOOP STOP if N > 2 38 

Jump after decimal digit 

After • store 20a or (240a + 12b) 

Clear X7 

Read character from tape, c , to 5 

c - 10 

3 

Jump if decimal digit 



c 3 -12 



Jump if not • 
Point indicator to 1 



Clear 


• 2 


after point 


Set X3 = 


= -1.0 


+ 20 


to 


6 c 


Jump 


if 


first • 


(0.1, 


0) 


to 3 


+ 12 


to 


6 c 


20a + b 


to 7 


Jump 


if 


second • 


LOOP STOP if third : 


c - 
3 


14 





Jump if not Sp 

Ignore Sp before decimal digit 



Number store f Overwritten by LINK 
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from 4+. 6+ 



0.0 



to 0+.0 

-« 



£2 

2 



to 3+. 5 

-< 



6 



5+ 



o 



0.3+ 5 60 3 

5 43 
0.7 5 60 
fis) 5 43 
0.2 5 60 3 
0.2+ 60 



16 5 00 

Qz) 5 43 
0.2+ 5 61 
1.7 7 01 
0.6 64 
r~^*j 0.5+ 60 



0+0 72 
0.0 60 



3+072 
0.5 60 
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Read Terminating Character 

Jump if Sp 

c - 15 

3 

Jump if Er to ignore Er 



c - 30 

3 



LOOP STOP after CR in a or b 
LOOP STOP on other punching errors 



Read character from tape, c , to 5 



c - 13 

4 



Jump If CR without LP or Sp in number 
| 240c + 126 + c | to 7 
Jump if 0VR clear 
LOOP STOP on overflow 

To B 0+ to complete formation of number 



Return to read rest of number after 
reading Er 
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square; root 

Method 

An iterative Newton-Raphson process is used to find 



y = Vx, where x = p + 2" 



3 8„ 



^ = y + d where d 
1+1 n n n 



~ 2 Vn^ n 



with y = i « + i . 

2 2 

In fact, y Q = - p + — , and the process terminates when d > 
to be the answer. There is special treatment of the case p = 



R 200 



Sheet 1 of 1 



y n is then taken 
when the 



root is also (1 



2" 3e ). 



(1 - 2' 



■ a e i 



01 



R 





-0 1 


200 


- 


28 - 


0+ 


m 


72 


0.0 





60 



LN 

SQUARE ROOT 




LOOP STOP if p < 



Overwritten by q 



p to 5 

q to 0.0 

p to 7 

\i> to 7 

v = (| *> + j) or y* + i = <y„ + d„) to 0.7 

Test OVR and exit if p = (1 - 2" 39 ) 

q to 7 

x/y n to 7 

(*/y„ - y„) = 2^ to 7 

dn to 7 

Test d n 

vfto 6 

Plant LINK and exit 

Overwritten 
by LINK 



Initially j f Overwritten by y n 
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EXPONENTIAL 



1 + ax + bx 2 + ex 3 



,2X 



Method: Let y = ■■ ■■■ , be an approximation toe in the range 

1 - ax + bx 2 - ex 3 

-k-£x£ k. 

Expanding log y, it can be shown that 

x - | log y - ax + (c - ab + | a 3 )* 3 + {a 2 c + ab 2 - be - a 3 b + | a 5 } x 5 + . .., 

and, equating coefficients: 

a = 1; b = |-; c = fe coefficient of x 7 = j^ . 

If we set a = 1 + a, 6 = f + /3, c = j^+7, where a, fi and y are small, 
then the logarithmic error is 

S = | log y - x = a* + [y - (/3 + -| a + a/3) ■ + (a + a 2 + -*■ a 3 )] * 3 + ... 

Now we can make JS| as small as possible in the range -k •& x <$■ k by choosing a, /3 
and y so that S is a multiple of the Chebychef f polynomial 

T T (0) = 64i9 7 - 1126 5 + 560 3 - 70, 

where 8 = x/k. We find, in fact, 



a 


- 


1 - 


1 *« 
14400 ' 






b 


= 


2 
5 ~ 


300* 


1000 


17 / 




200,000 


c 


= 


1 
15 


300 


■ -J- ft* 
2250 


3 


200, 000 



k 6 - ... 

k 6 - 



These approximations are satisfactory only when k is considerably less than 1, but 
the following method allows the range to be extended. 



Sheet 2 of 3 



R 220 



Let X be the operand and define x = X/16, so that k - -fa, and |S| < 0.65 x 2~ 4lt 
approximately. Now 

e X = g 16X _ (y€ -8*)8 = y B e -Ut 

where 16|S| < 0.65 x 2~ no and can be neglected, so that we can write 

Y = e x = e 16X = y* 3 

within the required accuracy. If we substitute X/16 for x in the original definition 
of y we find 

1 + a (X/16) ■+ 6 (X/16) 2 + c(X/16) 3 
y ~ 1 - a (X/16) + &(X/16) 2 - c(X/16) 3 ' 

To preserve significant figures of X when evaluating y, form 

aX + cX(X/16) 2 



7} = 2 3 (y - 1) 

which gives y - 1 + 2~ 3 % 



1 - a(X/16) + 6 (X/16) 2 - c(X/16) : 



If we define V r so that y z = 1 + 2~ (3-r > 7 7 r then 



,3 

and Y = y 2 =1+77 



In actual operation i1 r lies outside the range ±1 so £ r = -j?^ is used instead. 



This gives C = — 



2~ 2 aX + 2 -1 °cX 3 



5 - 2 _5 aX + 2- ? &X 2 - 2~ 13 CX 3 ' 

^l = ^0 + ? ^0 ' 

4» 3 = L 2 + J ^>2 ■ 

so that {e x = 1 y 23 = i + Ha' 
Since k ~ Z" 1 , the last term in a, b and c is negligible and may be omitted. 
Then a = 1, 

l _ _2 1_ 9-e _ 1 2 -16 

5 300 1000 

c = J 1_ 2 -8 _ * 2~ 16 

15 300 2250 

Putting a = 1 and substituting z = 2~ 2 X, the formula for £ may be written 

z + 2~ 4 cz 3 



I = : — 

A + 2~ 3 <-^ + 2~ : 



2 6z 2 - 2~ 4 cs 3 ) 
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GD- 



01 



R 





-0 1 


220 


- 


28 - 


1+ 





72 


0.0 





60 
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LN 

EXP 





Fl 


s 


7 6 03 




0.0 


..©!."... 




33 6 01 




6 7 26 


7 6 00 


2 


(T) 7 51 


6 7 22 


3 


6 7 00 


4 


"HS'l'li"' 




6 7 22 


5 


^T) 6 51 




6 6 22 


0.7 6 01 


6 

7 


0.7 1 10 


+.25 




^ / 




| 1+ 


0.0 


/Q) 6 51 \ 


, 


I 0.0 6 10 J 


( 6 6*21 \ 




0.1 6 10 J 




0.7 6 21 


2 


0.0 6 21 


0.0 6 11 




0.1 6 00 


0.6 6 21 


4 
3 


0.0 7 00 


Of (0[ 72 
0.0 60 






12 6 31 3. 


6 


11 5 41 5 




4 21 0. 


7 


13 5 00 



► -| + 2" 3 (-z + 5e 2 - Cz 3 )— >6 



£ ->7 



} 



£^8 



>■ *=2 



C->6 



I + i C 3 = i exp P-+6 



Plant LINK and exit 



► Overwritten by LINK 



► l£ = z->6 and 0.0 



z z — >6 and 0. 1 

Cz 2 —*6 
Cz 3 — >6 
z + Cz 3 — »-6 

5 S 2 -+6 

z + Cz 3 -V7 



B = 2~ 2 b (^2~ 2 .f) 



C = 2T*c (^2~\&) 
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LOGARITHM, WIDE RANGE 

This subroutine has been superseded by R 224 for most applications, but it is 
still used by certain other subroutines. 

The subroutine evaluates 



P< = ^-log a (fq) 



-U6 



in the range 2 " ^ (fq) < 1. 

The method used is the same as that described for the programme of R 224, but only 
four terms of equation 11 are used in R 221. 



01 



R 





-0 1 


221 


- 


28 - 


1+ 


® 


72 


0.0 





60 



LN 
LOG 









1+ 




(T) 5 42 
0. 0+ 6 63 


Cm \ 


\ U1 ) 


0.0 

1 




® 


5 56 






0.5 


6 14 




0.6 


6 21 




2 


0.7 


6 11 




0.5 


6 00 




3 





7 00 






0.7 


6 25 


to 0+.5 


4 


0+ 


§72 




5 




[+97184015999 ] 




6 


+47161592838 




7 


(+16674141054) 



-2 to 5 

LOOP STOP if number negative 

y to 6, a - 2 to 5 

(y - K/2) to 0.5 

(3-2 /2)y to 6 

(3 - 2/2) y + (%/2 - 1) to 0.7 

(y - 14/2) to 6 

Clear 7 

z to 7 



= V4/2" 



= 3 - 2/2 



= X/2 - 1 



<■ Overwritten by (y - 14/2) 



„ Overwritten by 
(3 - 2/2)y + (%/2- 1) 
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from 1+.4+ 



o.o 



0+ 



-95265423098 



+362106 



+16333229 



+925542347 



+94323185383 






= -\ log g 2 = -0. 34657359027 



T 7 = +0.00000131733 



T s = +0.00005941994 



T 3 = +0.00336710344 



T t = +0.34314574944 



z to 0.5 
z z to 0.6 



-g^(a - 2)log g 2 to 7 



- Overwritten by z 



• Overwritten by z l 



Overwritten by 
►-^(a-2)log e 2 









2+ 




0.7 


7 10 


1.0 






e 


1.0+ 


6 61 




1.6+ 


65 


i 








0.1 


6 00 




0.6 


6 21 


2 








0.2 


6 01 




0.6 


6 21 


3 








0.3 


6 01 




0.6 


6 21 


4 








0.4 


6 01 




0.5 


6 21 


5 








0.0 


6 03 


6 


© 


6 51 




0.7 


6 03 


i * 


A. 7 


1 10 \ 


7 


(,7 


60 j 
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-$%(a - 2)log g 2 to 0.7 

LOOP STOP if number too small 

Dummy order (Due to rounding OVR is never set) 



Sum series in z l 



Multiply by z 



i(log g y + 2 log g 2) to 6 
32 lo S e (A7> to 6 



Plant LINK 

Jump to obey LINK 
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FLOATING-POINT LOGARITHM 



01 



R 


-0 1 


223 


- 28 - 


0+ 


Of) 72 


0.2 


60 



m 








FP LOG 






R 


5 


-0 


1 


221 


- 


01 


- 



Given j> = A, 2 

(Standard floating-point form) 

This subroutine evaluates 



i>' - log e i> 



log. A + a log, 2 



Tag calling for R 221 



R 223 



Sheet 1 of 1 







0+ 


0.0 


127 7 77 7 
127 7 00 


1 


+256 


f nt '"I i 


7 00 
6 4 00 


ULJ 2 > 


3 


0.0 6 05 

6 4 03 


4 


1 1+ 5|71 
0.6 1 00 . 


S 


+0 


ft 


1 « 


0+ 072 
0.7 60 






7 


0. 1 4 03 
1+ Q]72 



Collating Mask for Argument A 

Constant added to exponent 

and round- off constant for argument 

Clear 7 

4.2° to 4 

A to 6 

a + 256 to 4 

Plant the LINK 

Set link for return from R 221 



Cue to R 221; ^ log e A to 6 



Link for return from R 221 
a to 4 



London Computer Centre, 
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1+ 


<s> 


7 40 


6 


7 20 


1.7 


4 22 


® 


5 40 


6 


5 56 


0.1 


6 01 


"0' 


5 56 


0.0 


6 05 


0.1 


5 01 


5 


6 01 


(" 


\ 


\> 


) 












+190530848196 



} 



+32 



2"" log e A to TQ 

2~ 3s [log A + a log. 2] to PQ 



5 c = +38 
Normalize 
Round argument 
Re-normalize 



Pack result in 6 in 
standard floating point form. 



Overwritten by LINK 



© 



= log. 2 
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LOGARITHM, VARIABLE RANffi 

This subroutine evaluates 



p' = — log fi (pq) 



(pq) is first normalized to give an argument y and an exponent o. Only the most 
significant 39 bits of y are used in the calculation. The required logarithm may be 
written 



log-y + a log e 2 where \k < y < %. 



We have*: 



V" t r 
log e (l-2t cos 6 + t 2 ) = -22_ — cos r 

If we put z - cos 8, then 

log e (l + t 2 -2tz) = -2^ p T r (z) (-14*4 1) 

where the T r are Chebycheff polynomials. 

Changing the sign of z and using T r (-z) = (-l) r T r (z), we have 



log e (1 + t 2 + 2te) - -2^ (-l) r £■ T r (* 



r=i 



Taking the difference of (2) and (3) 
"(1 + t*> + 2tz 



log. 



If we now put 



_2t_ 
1+t 2 



(1 + t 2 ) - 2tz 
• = A, then 



4 L-27+rv +1 w 



r=i 



l0g e{l-rw) = 4 Z_"5FPr T ir*i<*> 

> * +■= ft 



where k 2 < 1 and t = t(1 - /l-A 2 ) 
If j! is any constant 

loge i— = log^ + 4^ — t^,) 

It is required to choose k and A such that this series will give log e y. 
We require 

y - A.- — which gives 

1 — kz 

Z k y + A 



(1) 



(2) 



...(3) 



(4) 



(5) 



(6) 



...(7) 



...(8) 



where 14 4 y < & and -1 4 z 4 1 
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For y = % * - -l 
y = % a = 1 
Solving these equations, we obtain 



% -A = -k(% + A) 

k -i = ka + A) 



/T- i 



A - Yt^/2 and k = 

Substituting in (6) gives t = 

Substituting in (8) gives z = {3 _ 2 ^ )y + ^ _ 1} 

Substituting in (7) 



vT+ 1 
3 + 2/2 - 2/4 + 3V2 
y - i4i.?i. 



3 ^"t zr+l 
log c y = -jloM + *2_2m T «r*iW 

r=0 



R 224 
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...(9) 
...(10) 



...dl) 



The first five terras in equation (11) are used by the subroutine to evaluate 



log e y. The truncation error introduced by taking only five terms is about 2 
o log.2 is added to log.y to give the required log. (£<?). 



-37 



* "Tables of Integrals and Other Mathematical Data" - H.B. Dwight (Macmillan, 
New York. 1947). 



01 



PERRANT I LTD 
PEGASUS LIBRARY PROGRAMME 



R 224 



R 


-0 1 


224 


- 28 - 


i+ 


(0J72 


0.1 


60 



LN 

LOG. MK. 2 



R 2 


7-0 1 


224 


- 06 - 


R 2 


0-0 1 


224 


- 02 - 


R 2 


6-0 1 


224 


- 02 : 
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Title of Optional Parameter List 



Calls for P. P. 01 









,1+ 




0.0 


+97184015999 




<h 


C\ J. 

MM 5 42 


^ oi y-- 






CS 


0.1+6 63 




2 


(Too) 5 56 
0.0 6-14 




0.6 6 21 






0.7 6 11 




0.0 6 00 






0.7 6 26 




0+072 


to, 0+.6 


8 


2+072 




6 




+47161592838 




7 


i \ 

1 +16674141054 

i 1 



A - %v^ ■ Overwritten by (y - &/2) 

-2 to 5 

LOOP STOP If number negative 

y to 6, (a - 2) to 5 

(y - y 4 v^) to 0.0 

(3 - 2v*2)y to 6 

(3 - 2V5)y + (%,/jj- 1) to 0.7 

(y - #i/2) to 6 

z to 7 



= 3 - 2 vH 



= n/5-i 



Overwritten by 

(3 - 2i/2") y + (%/2 - 1) 
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0+ 




0.0 


-95265423098 




1 


+8415 




2 


+343173 




3 


+16347428 




4 


+925538404 




5 , 


+94323185678 


from 1+.5+ 


t 


/o.6 7 10 A 




6 






V 7 7 21 J 




/0.7 6 10 A 




1 


\ 0.0 5 20 J 



l lo ec 2 



•0. 34657359027 



+0. 00000003061 



T ? = +0.00000124845 



T s = +0.00005947159 



T 3 = +0, 00336708909 



T 1 = +0.34314575051 



z to 0.6 



z i to 6 



R 224 



Sheet 4 of 4 



z l to 0.7 

-Ua-2) log e 2 to PQ 



Overwritten by z 



' Overwritten by z i 







2+ 


1.0 


/fl2^ 55 o\ 
I 1.7 1 10 J 


1 


1.0 7 10 

1.1+6 61 


2 


1.6+0 65 
0.4 5 02 


I > 

3 


7.5+6 01 5 
0.7 6 21 


4 


1.3 5 66 
0.5 6 01 


5 


0.6 6 20 
0.0 6 03 


6 


(^)6 51 




1.0 6 03 


7 


7 5 00 OA 

\ ° / 



.Overwritten by 
2" M (a-2) log c 2 



+ P.P. 01 Shifts down n-1 places 

Plant LINK in 1.7 

-2" n (a-2) log g 2 to 1.0 

LOOP STOP if number too small 

Dummy Order (Due to rounding, OVR is never set) 

4 to 5 A 



Sum series in z 



Multiply by z 

log e y + 2 log e 2 to 6 

+ P.P.01. Divide C(6) by 2" 

2"" los e (fiq) to 6 

= Optional P. P. 01 (rt = 5) 



Overwritten 
by LINK 
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FLOATING-POINT EXPONENTIAL 

The method is explained in Section 2 of the Library Specification. 
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R 


-0 1 


225 


- 28 - 


0+ 


072 


1.0 


60 



LN 

FP EXP MK 4 



R 


5 -0 1 


225 


- 02 - 


R 


5-0 1 


225 


- 02 - 


R 


6 -0 1 


225 


- 02 - 


R 


7 -0 1 


225 


- 06 - 


R 1 


4 -0 1 


220 


- 01 - 



y Calls for P. P. 01 



Title of optional parameter-list 



Call for cue 01 to R 220 







0+ 


( ni ^ t. 


^1.7 1 10 \ 

y 6 i oo I 


1.0 


1 


1.5 6 05 

1+ [0]72 


2 


6 1 06 
1.6 1 01 


to 1+.0 3 


(38) 1 41 

0.0 1 62 

20 

0.2+ 60 


to 1+.2+ 4 




5 


+0 


e 
from 2+. 5+ 


+0 

x - 


or 2+.7+ „ 

7 


(- 


>6 

/ 



Plant LINK in 1. 7 
x = A. 2 a to 1 



A to 6 



Overwritten by 



- I (-Ioi72) [= - (6 - 2)] 



a + 2"" 1 to l c 
a to l c 
a + 38 to l c 
Jump if a > - 38 

If a < - 38. clear PQ 
and jump to find exp (0) 



+ 2 x P.P.01 = -2 71 .2" 38 



+ P. P. 01 = -2 n " l .2* 38 



Optional P.P.01 = -2"" 1 .2" 38 L Overwritten 
where n = 9 [by LINK 
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R 225 







1+ 


from 0+. 3+ 


(75) 1 44 
0.5 1 63 




0.0 




f 

1 




(25) 1 52 






0.7 6 20 






(jf) 55 1 


from 0+.4+ 


2 

»- 

3 


1.0 6 12 
0.6 7 21 




4 


0.5 1 00 




-H) 












2+ [5| 72 




5 










0.0 60 






6 




+190530846196 






7 


+198282496599 










2+ 




0.7 6 63 








0. 




1 


1.0 1 02 




® 1 41 
© 1 56 




1.6 6 03 




2 


© 1 56 




1.6 1 03 




3 


0.7 1 63 




i 






1.5 6 05 




4 


1 6 01 




1.5 1 05 


t£ LINK 


5 


1.7 1 60 




e 

6 


0.6 60 









6 00 
\1 o an 


to LINK 7 
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(37 - a) to l c 
Jump If a > 37 

(37 - a) to 1^ 

1 



x. 



log e 2 



(fixed-point form) to PQ 



-I = - Integral part of */log e 2 to 1.0 

P log e 2 to 6 

Set link in 1 for return from R 220 



+ cue 01 to R 220; Forms exp (P log g 2) 



Link 



= log.2 



2 log e 2 



Jump if floating-point underflow 
b 



(**h) 



+ 2 to 1 



Normalize 

Round argument, B, in 6 

Re -normalize 

b + 2"" 1 to 1 

Jump if floating-point underflow 

Pack result in 6 in standard floating-point 
form, B,2° 

Jump to obey LINK if in range 
LOOP STOP if floating-point OVR 

Clear 6 if floating-point OVR 
Jump to obey LINK 
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SIN/COS 

This subroutine evaluates sin ttx or costtx for -1.0 •£ x < 1.0. 
Method 

The method is based on the Chebycheff expansion for 



R 240 
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cos -^y in the range -1 ■£ y < 1. 

4* 

00 

cosfy ~- Jo(f)+2 2 J 2 n(fK^> 

The truncation error in taking the first seven terms of this series is about 

7.5 x 10~ 13 (< 2- 40 ). Re-arranging the first seven terms as a power series in y, we have 



cos 



77 ""^T - 



2r 



In this series c, < -1 and it must therefore be re-arranged to read 



77 

cos -y 



2 fc r> 



zr 



r=o 



-y 



where 6 1 = c, + 1, 6 r = c r for r ^ 1. 

To reduce rounding errors in evaluating the series, it has been re-arranged to read 



where 



7T 

cos -y 



l r> 



2T 



r=o 



- y + a. 



a - — 6 , a r = 26 r for r ^ 0. 



Overflow 

If 

leaving 
but 

so that 



|x| > - 2x will overflow, 



y = 2x ± 2 and 0VR set 

COS = - COS(£ ± 77) 

77 
COS TTX = - COS ^ <2x ± 2) 

77 

= - cos -y . 



Thus, if OVS is set in forming y, the series for cos — y is still used but the result is 
negated and OVR cleared at the end of the routine. 



To find sin vx 



sin 6 = cos(e ~^\ = cos [9 + -^j 
so that sin 77JC = cos n(x J = cos t/x + -J. 
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The cue for sin rrx therefore subtracts - from x and then uses the cosine series as 
before, but with y - 2 x - 1. 

Note that (x - ~) will overflow if -1 ^ x < -\, leaving (x + -j) in X6 and OVR set. 
The fact that OVR is set does not matter because 

(x + § ) > \ 

and overflow will occur when x + -| is doubled to form y; thus the fule given in the 
previous paragraph will still apply. 

Special Cases 

If y = -1.0, y 2 will overflow and become negative (y 2 = -1.0) and the series will 
give the wrong result. For y = -1.0, cos Try- 0. The routine therefore tests y 2 and 
sets the result equal to zero if y 2 is negative. 

IT —38 

If y = the series will give the value cos — y = 1 - 2 



01 



02 



R 





-0 2 


240 


- 


28 - 


0+ 


§ 


72 


0.0 





60 


0+ 





72 


0.0+ 





60 



LN 
SIN/COS 
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GD 



o.o 



33 6 03 

(l) 6 50 

1+Q] 72 

6 6 21 

6 5 00 

1.1 6 00 



1.1 1 10 

0.4 1 00 

I 5 6 21 

1.0 6 01 1 



0.4 i 66 
6 51 



0. 



5 6 03 
1.7 6 01 



1.0 5 62 
6 00 



l.o 



1.1 64 
6 6 02 

+246961 



-13838561 



+505357534 



-11469815955 



+139456286127 



-128478236105 



+183251937962 



x - - to 6 

2 

y - 2x or <2x - 1) to 6 



y . to 6 

y 2 to 5 



a 6 to 6 



Plant LINK in 1.1 



1 ? =2 



,?r 



► Form ^ a^ In 6 






- y + a in 6 



If y 2 - -1 (OVR when y = -1) clear 6 



-N 



EXIT if OVR clear 
Negate 6 if OVR set 

a, = +0.00000089843 



= a 5 = -0.00005034439 



= a 4 = +0.00183847999 



a 3 = -0.04172694736 



= a 2 - +0.50733901344 



= a, = -0.46740110012 



a Q = +0.66666666666 



► Overwritten by LINK 
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We have :- 

2 00 

1 - t — 



= 1 + 2 £] t r cos r (t 2 < 1) (1) 



1 - 2t cos 6 + t 2 

r-i 

For cos 6 write z, then cos rO = T r (z) where r r are the Chebycheff polynomials, and 
we get, for -1 < z < 1, 



1 - t z 



(1 + tz) - 2tz 



= 1 + 2 Y^ * r T r<*> < 2 > 



r"i 



Change the sign of z and use T r (-z) = (-l) r T r (2). Then 

2 °° 

Now take the difference of equations (2) and (3):- 

.2 



4(1 - t')tz , X -1 J .2f+i „ , s 

(1 ftV - 4tV = 4 L * ViW 

r-o 

9 00 



r«o 



Change the sign of t : ■ 



(1 + t 2 )z 



d-«v + «y = S'-'^w 



Now integrate with respect to t from tc t: ■ 
t 



1 + t z 



where t = tan—. If we put k = tana, then:- 



where 



(4) 



(5) 



d* = > (-l) r I . , T (2) (6) 

Z_^ 2r + 1 2r+i 



(i - *v +f jlV 2 ™ 

\i-tV * 

a. 2t 

Substitute t = tan « on the left-hand side, so that tan a = ,2 and 

X * V 

d(tan a) = 2(1 + t 2 ) 
At (1 . t 2 )2 * 

t 

1/2 1 

— / d(tan a) = — arctan (z tan a) (7) 

J 1+2 tan a 

o 



Thus arctan (2 tan a) = 2 2_ <-D 2 r" ~ T T 2r»i (ar) < 8) 



<^ .art 1 

arctan (kz) = 2 ^ (-D r 2r+1 T gm (2) < 9 > 

t = i (/l + k z - 1). 
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To find arctan of a given number, reduce it to the range 4x4 1, the raid-point 
of which, in terms of angle, is 77/8. 

Put arctan x = — + 9 so that - — ■£ 6 < — - 

and put arctan ky - 8, choosing k so that -1 ■$ y < 1 

i. e. k = tan — = </2 - 1. 

o 

7f 

Now, arctan ky = - — + arctan x 

o 

= - arctan (•/T - 1) + arctan x 

x - (v/T- 1) I 



= arctan - 



1 + x (^2- 1) 



y'^T' iV*^". 1 !) = u^r-~" + "-'»'* ~) ,tnd - 1 * y<1 (10) 



* - b/T- n V 

- 1) + (3 - 2/2) x J' 



The quantity y can thus be found from x by means of a multiplication, an addition, a 
subtraction and a division and we then use:- 



arctan 

» 



~ t 2rn 



where t = tan ~ = A + 2</2 - (i/2~+ 1) 

lb 



We can now write: - 

2r+i 



i-arctan* = {+£ <-D r % 2r + . 1} T 2r+1 (y) 



(12) 



r=o 



N.B. (3 - 2t/2) = 0.172 and (v/2~- 1) = 0.414 so that there is no danger 

of exceeding capacity during the formation of y from x (0 •£ x < 1). 

-12 
The truncation error in taking the first six terms of the series is about 1.2 x 10 

(= jx 2 _3B ). 



1. "Tables of Integrals and Other Mathematical Data" - H.B. Dwight (Macmillan, 
New York, 1947J. 
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p' = — arctan (p) 

or p' = — arctan f — J 
where -14 « j> ' 4 14 



R 241 
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R 0-0 4 




241 - 28 - 


01 


1+GD72 
0.0 60 


02 


l+[o]72 
0.1+ 60 


03 


0+ 00 0. 




04 




0+ 00 0. 



LN 
ARCTAN 









A 


Cni *\ 




6 4> 00 


I 01 J 


0.0 






32 7 00 
0.2+ 60 






flF\ 


I 


6 4 00 


\®D 






7 4 06 




2 






^ * 


33 2 02 




*■ — 


0.4 6 63 




3 






^ » 


6 6 02 




0.5 7 63 




4 


7 7 02 




L — > 


6 7 03 




5 


3+ [I] 72 
(T) 2 53 


$o 3+;2 


6 


1.2 7 62 








7 6 01 


j;o 3+.1 


1 


11 fi en 











p to 4 [Entry for — arctan (p)] 

- y = -1 to 7 

£ to 4 [Entry for — arctan (*-j\ ■ 
Sign 4 = sign (pfq) 

- 14 to 2 

- |*>| = -u to 6 

' - \q\ - -v to 7 

u-^ to 7. If u = 1, v = 0, 7 = -1 St 0VR set 



a = | to 2 



Jump if u ^-v (except u=l, v=0) with 

-v to 6 

Jump if u = 1, v = 



fx6 = -« 
1 X7 = u-v 
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A 
B 

C 



tan 



Z- - 



(2 - /2) 



= k tan 2j | = k(3 - 2/2) 
= k(l - tan|-) = k(2 ~ /2) 







3+ 


from 1+.7+ 


33 2 03 
32 4 06 


1.0 


from 1+.7+ 


7 7 02 
1.3+0 65 




l 
— — ♦ 

2 

c 

3 


from 1+.6+ 


7 3 00 
1.2+ 6 60 




6 3 26 

7 3 00 






4 


1.6 3 01 

1.7 7 21 




5 


2+ GO 72 
0.1 60 




f 
6 

7 




+161019749891 




+47161592863 














2+ 




0.0 


r ^ 

+161019749976 

0.0 6 01 
6 3 26 




1 


2 


0.0 7 10 

7 7 21 


3 


O+Q 72 
6 3 00 


4 


1.1 6 00 

0.5 5 00 






\ 3 6 21 
1.0 6 01 5 




5 

6 

7 


to 0+.0 


0.5 5 66 

0.0 6 21* 

2 6 04 

11 1 1ft 











Mark for interchange by putting a 
and reversing the sign of 4 



■5 in 2 



"{s = !:S = S^*«» 



v - u to 7 

Jump if « = 1, ii =0orti =0, v - f 

LOOP STOP if u = V = 

(f - 1) to 7 where g = * or J , < £ « l 

<f - l to 3 

^ - 1 +J to 3 
5 (£ - 1) to 6 



sei 



A - +0.58578643762 



B = +0.17157287534 



+0.58578643793 



B (<f-l) + C to 6 

g - 1 + A 



A, B and C 
as above 



Overwritten 

by *■ 



g - tan tt/8 



5(^-1) + C k (tan tt/8)(1 + g tan tt/8) 
Plant s in 0.0 

y* to 6 



to 7 



y2 to 3 

o 13 to 6 

Set counter of 6 i.e. 5^, = 2 



Porm|r^ a 2r+1 y 2r+i in 6 



r=o 



Plant LINK 
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0+ 


from 2+. 7+ 




1.1 4 62 




f * 

1.0 






6 6 02 




/ \ 








l 


+42077 




2 


-450949 




3 


+3461449 




4 


-26138321 




5 
6 


+213371737 




-2072724634 




7 


+36242177026 



Jump to LINK if 4 is positive 
i.e. p' = (I At) arctan (p) 
p r = (1/tt) arctan (p/q) 



a 13 = + 0.00000015307 



► Overwritten by LINK 



a u = - 0.00000164054 



a, = + 0.00001259267 



a 7 = - 0.00009509065 



a 5 = + 0.00077624185 



a 3 = -0.00754052829 



a x =+0.13184827194 
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ARC SIN /ARC COS 

The subroutine evaluates, using R 241 and R 200 as subroutines, 

y - - arc sin x - - 4 y 4 — 

TT 2 2 



Sheet 1 of 2 



or y = — arc cos x 

J 7T 



<? y < 1 



where — arc sin x = — arc tanf ■ 

TT 7T I /, 2 



V*-« 



and — arc cos * 

77 



arc tan f - 

? 7T I 



Y^ 



Initially p = x and the result, y, is left in X6. 



01 



02 



03 



04 



R 





-0 4 


242 


- 


28 - 


0+ 





72 


0.1+ 





60 


0+ 





72 


0.1 





60 


0+ 





00 0. 
















0+ 





00 0. 



arc sin 



arc cos 



LN 

ARC SIN/ARC COS 



R 


5 -0 1 


200 


- 01 - 


R 1 


1 -0 2 


241 


- 01 - 



Call for cue 01 to R 200 



Call for cue 02 to R 241 
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o.o 



+0 



1+ jf) 72 
1.2+ 60 



1+072 
6 7 00 



l.o 



1+ 



+0 



~0< 



1+072 
1.6 60 



1 



Sheet 2 of 2 



3 (8.2, 0)- 



Modifier 0.2 for 00 order In 1.0 
Modifier 0. 1+ for 60 order in 1. 3 



or (8.2, 0) to 3 
Plant LINK 
-x to 4 

(1 - x 2 ) to PQ 

Jump if * = or -1.0 
Set LINK 



Cue to R 200 



; /l - x 2 to 6 



LINK obeyed in 0.7 of R 200 



/T~T 



x to 7 



x to 6 

Set LINK from 1.5 or 1.7 



Cue to R 241; — arc sin % to 6 



* = or -1.0 to 6 

— arc sin x - or — to 6 

77 2 

Jump to 1.5+ (arc cos) or 1.7 (arc sin) 

38. 



+ (1 - 2 
EXIT 



) to 6 and set OVR 



LINK for arc cos 



— arc cos x to 6 

77 



Jump if * = -1.0 and clear OVR 



LINK 
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COSINE AND SINE 




This subroutine evaluates 

p ' = COS TTp 

and q' = sin np 

where -1 ■£ p 4 1 

and wp = (7) + £) tt/2 



V 



-2, -1. 0, 1, 2 arid -% ^ £ < % 



Let 



tan f if .f 



(^r 



/'(f 2 ) 



where />(£ 2 ) = a + a 2 £ 2 + a, £» + a & £ 6 



Then cos f = c = 
and sin £ - s = 



1 - t 1 

1 + t £ 

2t 



.2 ' 



1 + t k 

cos pu and sin ^w are then evaluated using the following relations: 





T? = -2, 


1""" 

2 f 
! 


TJ = -1 


— 1 — 
1 
| 


Tj = 


1 

t 
1 


7J= 1 


COS pTT 


-c 


1 

1 

_|_ _ 


S 


1 
I 


c 


1 
1 


-S 


sin Ptt 


-s 


1 
1 


-c 


1 
1 


s 


1 
-4 


c 


OVR 


Set 



1 

1 
1 
1 


Set 
0.1 


1 
1 

--!■- 

1 
. J 


Clear 



1 

1 
1 


Clear 
0.1 
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01 



02 



R 


0-0 2 


243 


- 28 - 


0+ 


00 0. 












0+ 


00 0. 



LN 

COS + SIN 









0+ 


ENTRY 




33 7 


00 




3.0 










33 7 


22 




32 6 


01 




1 


G> 


52 






3.3+ 6 


62 




2 






35 1 


01 




32 6 


06 




3 










33 e 


03 




5.3 6 


10 




4 










6 6 


21 






6 E 


, 00 




5 










4.0 4 


02 






5.4 £ 


00 




6 










5 e 


i .21 




5.0 e 


01 4 




1 








~ 


3.6+ 4 


66 







a - order partial cue 



b - order partial cue 



& to 7 

p + Vi to 6; OVR set if p 5> + % 

p - % to 6: OVR set if p < - >A 

2p ~ 3 /2 modulo 2 

Jump if*<-#or/>»+ttor-i4«/>< l k 

Add 0.1 to l m 

Change sign bit of 6 

£ to 6 

£ to 5.3 

t 2 to 6 

€ 2 to 5 

Set ip = 5 



a 6 to 6 



. Form P(i 2 ) = a + a^ 2 + a^ 4 + a£ 6 in 6 
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, LiL 


4.0 


1 5.3 6 21 
(2) 5.51 




32 5 01 




5 6 26 




7 4 02 


2 
3 


7 51 
4 7 21 
33 6 03 


4 


6 4 26 

7 5 00 


5 
6 


4 5 21 
4.6+ 6 61 
(7) 6 42 


^ -*" 


32 6 03 


7 


5 7 02 

6 6 10 1 









2+ 


5.0 


5 
5.2 


5 



10 1 
64 


1 


6 

7 

® 


6 

7 
1 


02 

02 
52 


2 










60 1 


3 


(" ) 


4 


-358550 





5 


-144697 


583 


6 


-958183 


6739 


7 


+215888 


603285 



Sheet 3 of 3 



€ P^ 2 ) to 6 

'A (S 2 ) to 5 

% (£ 2 ) - 1.0 to 5 

- t to 7 
t to 4 

- M to 7 

- &(t 2 + 1) to 6 



} 



- s 



2t 



/ 2t 2 \ 

" W^j t0 



(t 2 + 1) 
6 



to 7 and 5 



If C(6) = replace it by -2" 33 to avoid OVR 



to 6 



s to 7 

c to 6 or 7 



s to 5 or 6 



Negate s and c if OVR set 

(i.e. if original angle is > 94 or < -&) 



Set Computing Store Link modifier 
Obey LINK 

Overwritten by £ 



= a 6 = -0.00001304397 



a 4 = -0. 00052640674 



= o 2 = -0.03485851898 



a = +0.78539816344 
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FLOATING-POINT ARCTAN 

An accurate series cannot be obtained for arctan (plq) over the full range of 
- 004 p/q ^ 00, it is possible, however, to evaluate arctan for some reduced range 
of p/q and then transform the result to the required value. 

R 251 uses the following rules to reduce its operand to the range 

- tan 77/12 * x 4 tan 77/12. 

The original range is 

- 00 ^ p/q ^ oo, 

(i) If p > 0, arctan {plq) - - arctan (- p/q). 

If p > the routine negates p and sets a marker to indicate that the 
arctan then obtained must be negated, p 4 throughout. 

(ii) If q < and p 4 0, arctan (p/q) = - 77 - arctan (pi- q). 

The routine makes q ^ and sets a marker to show that the resulting 
arctan must be transformed. 

The range is then - <x>4 p/q <$ 0. 

(iii) If p 4 0, q > 0, arctan (p/q) - - n/2 - arctan (q/p). 

If \p\ > M. the routine interchanges p and q and sets a marker to 
indicate that the above transformation must be made to the result. The 
range is then - 1.0 ^ x 4 0. 

(iv) If x(- tan 6) is less than - tan (7T/12) the routine forms a new variable 
X using the formula 

„ ■ x + tan 77/6 
X - tan (9 + 77/6) 



1 - x tan 77/6 
x + 1/1/3 



1 - x//3 

and sets a marker to indicate that it must later subtract 77/6 from the 
tan" 1 X thus obtained to give 8. 

The range of the operand is then 

- tan 77/12 ^ X 4 tan tt/12. 

Six terms of a Chebycheff expansion are used to evaluate tan -1 X in this range 
and the appropriate transformations are then made to form tan" 1 ^/?). 
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R 


-0 4 




251 


- 28 - 




0+ 


072 


01 








0.2 


60 




0+ 


©72 


02 








0.2+ 


60 




0+ 


00 0. 


03 














04 








0+ 


00 0. 



LN 

PP ARCTAN 



R 





-0 1 


251 


- 


02 - 


R 


1 


-0 1 


251 


- 


02 - 


R 


1 


-0 1 


251 


- 


02 - 


R 4 


1 


-0 1 


251 


- 


02 - 


R 4 


1 


-0 1 


251 


- 


02 - 


R 4 


2 


-0 1 


251 


- 


02 - 


R 5 


7 


-0 1 


251 


- 


06 - 



Tags calling for parameters 



Title of optional parameter list 
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The operands are in packed floating-point form, 
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Coo- 



to 3+.6+ 



t o 3+.6 



P = A.2 a 



0.0 



2 



0+ 



/3.7 7 77 7\ 

1 7. 7 7 77 6/ 

0.0 7 02 

7 5 00 



0.1 7 05 
1+ [1)72 



7 5 06 
0.0 6 10 



0.1 6 05 
0.0 6 16 



0.0 5 04 
4 40 



© 



3+ ED 72 
1.1 6 63 



1.0 



35 4 06 

6 6 02 
1.3 7 62 

33 4 06 

32 02 

7 7 02 
0.6+ 5 63 

6 3 00 

7 3 01 
1.5+ 5 61 
0.6+ 3 62 

5 5 02 

32 02 

32 4 06 

3 6 04 

0.6 60 



q = 5.2 £ 



+ P. P. 01 = -l (f.p. ) [-Overwritten by A.2 a 



P.P.) V 



+ 2 x P.P.01 = -2 n .2" 3a 



+ 1 to 7 [Entry for arctan (j))] 
q - B.2 b to 5 [Entry for arctan (p/q)] 
B to 7 

b to 5 

p - 4.2* to 0.0 

A to 6 

a to 0.0 

(a - b) to 5 

Set marker; 4^ = 0, 4 C = +2 

Jump if p < 



4 m = 0. 1 as marker to indicate p > 

-M 

Jump if <7 ^ 

Set marker for <7 < 0; set 0VR for q < 

+ [5| 

Jump if a < b ( i.e. !*>| < k|) 

- Ml to 3 

\B\ - \A\ to 3 

Jump if a > b (i.e. \p\ > \q\) 

Jump if \B\ 5- Ml and a = 6 (i.e. |/>| < |q|) 

(6 - a) to 5 

Set OVR if interchange 

+ MI 
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3+ 




0.0 
1 
2 
3 
4 
5 


-6660229215 




+9621296443 




-12495809411 




+17499210016 




-29165451385 




+87496355273 


from 1+.7+ 


/ 3 7 fid \ 


from 1+.3 
or 1+.5 


6 

7 


1O.6+ 7 
\6 + 


60 J 

H 1 

72 J 











6+ 




7 


6 26 




1.0 










7 


6 00 







7 00 




1 


© 


5 56 






6 


6 01 




2 


© 


5 41 






0.6 


6 10 




3 


0.7 


5 10 




5 


5 02 




4 


© 


5 43 




1.6 

© 


5 62 
5 50 




f 
5 

6 




© 


6 51 5 






6 


5 00 




2+ 


[I 72 


to 2+.2 


7 


1.2 


fiO 











- U\ 

LOOP STOP If p = q = 
4/2 



Sheet 4 of 6 



Overwritten by 
(4/ff) 



Overwritten by 

(|o - &l + 2) 



y A/2B to PQ 

Normalize 

A/B 

<-|a - 6| + 2) to 5 C 

(A/B) to 0.6 

<-|a - b\ +2) to 0.7 

|a - b\ - 2 

|a - b\ - 39 

Jump if I a - b\ ^39 

5 m - |a - t>| - 39 

Form fixed point * between -1.0 and 

* to 5 
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2+ 




1.0 
1 


+73653313182 




-79350416784 


from 6+.7+ 


1.0 6 01 
1.7 6 62 




2 




3 
4 
5 
6 




1.1 5 21 
33 6 01 




(T) 5 51 
1.1 5 03 




6 5 26 

7 5 00 




32 02 
(7) 4 43 




7+ [l] 72 




7 




■ 







to 4+.4+ 

-4 

to 4+.4 

_« 



7+ 



1.2 64 



1.0 



0.6 


5 10 / 


® 


6 40 


0.7 


6 10 


5 


5 21 


1.0 


6 10 


1.1 


5 02 


0.0 


6 00 


1.0 


6 21 


7.6+ 


6 01 5 


1.4 


5 66 


0.6 


6 21 


0.7 


5 00 


4+ 


§72 


0.4+ 


4 67 


0.4 


60 



^3 = tan 



12 



1 
2/3 

x + tan (tt/12) 

Jump If x < - tan (tt/12) 



x + (1/V5) 
» Form a - 7=- in 5 



- (xf/i) 



Set OVR (as marker) 



4 C = +1 



Jump if x ^ - tan 77/12 
Replace x by X in 0.6 

Exponent of X - 2 to 0.7 

x 2 (or X 2 ) to 6 

to 1.0 
Set 5^, = 3 
First coefficient of series to 6 



Sum Chebycheff series to form 



y 1 

77 



arctan 



x (or — arctan X J 



Exponent to 5 



Jump if x ^ - tan (tt/12) 
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from 7+.7+ 



from 7+. 7 



0.0 





4+ 


21 2 52 5. 


42 5 25 3 


+0 


+0 


+215888603272 


0.0 6 03 





5+ [Q72 


0.7 4 62 


6 6 02 


33 6 03 


© 452 


1.1 4 62 







5+ 




6 6 


02 


1.0 






L_^ 


32 6 


06 


1.6+ 6 


60 


i 






2 


0.3 6 


20 


@ 5 


56 




4 6 


12 4 


3 


© » 


56 


4 


0.2 6 


03 


© * 


56 




0.2 5 


03 


0.1 6 


05 


5 








5 6 


01 




6 6 


00 5 


6 






^ ►■ 


1.7 1 


10 




\l 


\ 


7 


I -25 


6 

/ 



1/6 



+ 2 x P. P. 01 = - 2". 2 



■ae 



1 n"38 



+ P.P.01 = - 2" -1 .2 



= 77/4 



— arctan x - — arctan X - — 



Jump if \p\ > \q\ 

1 1 



1 *• 

— arctan L - 

Tt q 



Shift marker 
Jump if q > 



■=r arctan x 

2 77 p 



— arctan %; 

7T q 



- 1.0 - — arctan ± 

77 p 



Jump if arctan (p/q) = 

Multiply by tt/4 to form -j arctan (p/q) 

Normalize 

If p < 0, 1p = 0, No effect 

If p > 0, ip = 2, arctan (p/q) = -arctan (-£/?) 

Re -normalize 

Round 

Re-normalize 



Pack result in 6 in standard floating-point 
form 



Set result zero if floating-point underflow 
Plant LINK 

Optional P. P. 01 = -2*" 1 .2* 38 1 Overwritten 
where n - 9 f by LINK 
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GAUSSIAN QUADRATURE 

This subroutine evaluates the integral 

rX+h 

f(x) dx 

where \h\ < 1. This integral may be expressed f in the form 

X+h 

n 

fix) dx = 2VV 

where the B r are the weights of the Gaussian quadrature formula and 

4> r = i{f(X + ht r ) + f(X - ht r )j 

where the t r are the roots of the Legendre polynomials, P n (x). 

In this subroutine, the roots and weights are stored as positive integers to 
preserve greater accuracy. 

t "Numerical Analysis" (Chapter VII) - Zdenek Kojial (Chapman and Sail, London, 2955J. 



01 



R 


-0 1 


300 


- 28 - 


0+ 


1372 


0.1 


60 


LN 




GAUSS 


. QUAD. 


R 1 


5-0 2 


300 


- 02 - 



Call for P. P. 02 
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o.o 



r 



? 



ED 

*j [?]73 

6 3 00 

2 00 

(o) 4 40 

1+ 072 

2+ [7] 72 



(7) 6 50 

0.5+ 64 
0.5 60 



0.0 6 10 
072 
0,0 6 00 
0.7 1 10 
0.7 60 



0.0 



f 



1+ 



1.1 6 21 4 
6 2 01 
35 4 01 
0.2+ 4 66 
2+ |7]72 4 
0.4 4 67 
2 3 21 
0+ [o]72 



1.0 3 21 4 
0.6 1 00 



+ 



1+ [0j72 
0.7 60 



2+ [1)72 4 
0.0 60 



} 



rX+h 
Overwritten by / f(x) dx 

Preserve Accumulators in BO 

A to 3 

Clear 2 

Becomes (n"^ 4 40. Set by Interlude 



Read first block of roots and weights 

n fX+h 

2h! Hf^ - J f(x) dx to 6 



LOOP STOP on OVR 



Restore Accumulators 

rX+h 

I f(x) dx to 6 
J X~h 

Plant LINK and exit 



Vr = \Pr{f< J + ht r ) + f(X " Ai r>} to 6 
I (^ r ) to 2 
Step modifier in 4 

Read next block of roots and weights 

Count n' points 

n 
h 2 ffjpr to 6 
roi r r 

At r (rounded) to 6 

Set LINK for return from Auxiliary 

+P.P. 02 = Cue to Auxiliary 

LINK for return from Auxiliary 
Read roots and weights 
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T 6+.0 



from Assembly 










7+ 




33 


4 00 




1.0 








34 


4 03 












1° 


5 


70 4 


1 










1 





70 1 












1 


2 00 


2 


© 


2 45 


c 


1.3 


2 60 


3 








33 


4 00 




35 


4 03 


4 


© 


1 52 




4 


1 05 


5 








1 


2 00 












1° 


o 


70 1 


6 








1 


3 00 




0.7 


7 00 


1 








3 


[|]72 7 



1 
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INTERLUDE 



Clear 5 

Set LINK for return to Assembly in 0.5 

Jump if required list not yet found 

Jump if last list 

(7+. 0, 0) to 2 

Plant new T.A. in 5.7 m 

Replace LINK in 1 

Overwritten by LINK for return to Assembly 



► Overwritten by new T.A. 



(7+.0, 0) 



i- 511.0 to 4^ 

Index Address to 1 
Index word to 1 
Index word to 2 

LOOP STOP if no parameter list 

Collating mask for last 
12 bits of modifier 

► Form parameter list address in l m and 2 m 
■■ 1st word of parameter list to 1 and 3 



7+.0 to 7 



m 



Read B10+ to U4 



to 10+. 1 
< 
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INTERLUDE (continued) 



R 300 



2.0 



r- * 

t 



I L 







4.1 1 62 
1 45 

1 3 04 

3 1 00 

35 2 01 

35 2 01 

34 3 01 

072 2 

1.0 6 00 2 

4.0 6 10 3 

2.6 2 66 

|l]72 2 
2.7+ 3 66 

0+ g]73 3 

1+ [4] 72 3 

2.4 3 67 





3.0 






9+ 




1+ 

3 


E 


73 3 

72 7 




34 


3 01 






4.0 


3 01 




3 


2 00 


from 10+. 3 


2 


1 


3 00 




? 




3 


© 


1 45 






1 


3 01 




|0f 


2| 70 




4 


© 


3 52 




3 


1 01 




5 










|0+ 


2| 


71 




2 


1 00 




6 










4.0 


7 


00 




h 


«l 


71 7 


to 6+.3 


7 


0.3 




or 6+.4 
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Jump if a 6, 8, 10, 12 or 16 point formula 

+1 to 1 if n odd; to 1 if n even 

n' = <n + 1) to 3 if n odd; n' = n to 3 if n 

even 

n' to 1 

Form address of 3rd parameter in 2 

(8, n r ) to 3 

Read first block containing roots and weights 



Copy roots and weights from 
parameter list into B2+ onwards 



Write last block of roots and weights 
Read B10+ to U4 

(2+.0 + n') to 2^ 

n' to 3 

1 to 1 if n' odd; to 1 if n' even 

n' to 3 

Read order-pair in 0+.2 to 1 

Add n' into If position of 6-order 

Restore adjusted order-pair in 0+. 2 

New T.A. to 1 

0+.0 to 7 ffl 

Plant new T.A. in 6+. 6 



from 8+. 



.to 9+.2+ 



to 6+.4 
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4.0 



10+ 



0+ - 00 0. 


5.7 2 00 

35 2 03 

10 01 70; 2 

1 3 03 
3.2+ 3 60 

4.0 2 00 

34 2 01 

34 2 01 

0.4 5 60 
4 [1] 72 7 

3 1 01 

1 7 00 

Uo) 6 40 

1.0 60 



1.0 




11+ 



.0. 



6 24 

1.7 5 00 

16 5 10 

5 53 

1.1 5 61 

1.3+ 7 60 

16 7 10 

16 6 10 

5.7 2 10 

1.6 1 00 

77 





59 (3]72 1 
3.1 60 



39 00 
3.4+ 6 67 6 
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INTERLUDE (continued) 

= (0+.0. 0) 

Transfer Address to 2 

Reduce T.A. by 0.1 

Read last word of list to 1 

Test list of roots and weights 
and jump if required 

Set T.A. for next list as 2+.0 

Jump if list rejected, unless last list 
Read B11+ to Ul 
n' to 1 
n' to 7 
+10 to 6,. 



Arrange to print n ' 
Printing constant to 5 



Print CRLF * 300 Sp + n 



Plant T.A. in 5.7 m 
LINK to 1 

STOP 



Self modified LINK for return 
to Initial Orders input 



Constant for printing if no 
parameter list supplied 
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LISTS OF ROOTS AND WEIGHTS 

J560. 0+ - 6+. 1 

+256315268353 

+47093317874 

+181751852204 
+99165386105 

+65591142427 

+128619202965 

+6 

J560.0+ - 6+.1 

+263962465814 

+27825488178 

+218986013844 

+61127633295 

+144457248869 

+86231026213 

+50422130590 

+99693759258 

+8 

J560.0+ - 6+.1 

+267705388118 

+18326479577 

+237786807609 

+41080874044 

+186754680092 

+60222000768 

+119130818817 

+74015472214 

+40922266697 

+81233080341 

+10 



J560. 0+ - 6+. 1 
+269809332680 

+12967457725 

+248521859063 

+29395258100 

+211629235633 

+44001995897 

+161440729985 

+55846237017 

+101108752553 

+64181939731 

+34423897211 

+68485018474 

+12 

J560.0+ - 6+.1 

+271964458139 

+7463611212 

+259642805294 

+17112118360 

+237942893098 

+26156972519 

+207643982665 

+34257750763 

+169840528812 

+41120632293 

+125898693188 

+46497389997 

+77406594626 

+50193644529 

+26116839838 

+52075787271 

+16 

J560. 0+ - 6+.0 
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LINEAR INTERPOLATION 
Method 

See Section 5 of the Library Specification, 



(jry 



C°l> 



to 1+.0 



01 



02 



03 



R 





-0 3 


320 


- 


28 - 


0+ 


@ 


72 


0.2 





60 


0+ 


§ 


72 


0.1 





60 


0+ 


§ 


72 


0.4 





60 



LN 

LINEAR INTERPOLATION 



R 2 


6-0 2 


320 


- 06 - 


R 


-0 1 


320 


- 02 - 


R 1 


4-0 2 


320 


- 02 - 



0.0 



0+ 



+0 



0.0 7 00 
0.3 60 



OD 7 52 
0.0 7 01 



jT]72 7 
1.0 7 00 7 



072 7 

1.7 7 10 



6 7 00 
6 00 



1.2 7 03 
6 23 



1+[0|72 
0.0 60 



Title of optional parameter list 



Call for P.P. 01 



Call for P. P. 02 



+ P. P. 01 



Address B.P. of index 



(B.P, 0) to 7 



Table number, q, to 7 



m 



(B.P + q) = address of q n word in index to 7 : 
Block containing index to Ul 
Address of table to 7 
1st block of table to Ul 
Address of table to 1.7 

x to PQ 

(x - x Q ) to PQ 
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1+ 




0.0 

x<x 


1.3 6 24 

7 5 00 

0.4 7 63 

1.1 7 03 




1 

* >X m 




0.4 7 62 
07 41 




2 


X — X 


0.5 6 60 
0.5 7 61 

+0 




3 








4 




(T) 5 41 
C~25) 5 52 


5 


« 


1.7 5 01 
1.3 6 26 




7 


2+ §72 
0.0 60 


















2+ 






/O.O 7 12 \ 

U— -< 

/0. 1 7 10 \ 
y 7 00 J 




0.0 
1 


2 


[T| 72 5 

1.0 6 00 5 






0.4 5 66 
[FJ72 5 




3 




0.0 6 22 
1.0 5 00 5 




4 


S 
6 


0.1 5 22 
0.6 1 10 

/ 2 - oo oA 




0.4 60 


7 






Sheet 2 of 2 



x - x 



r - 



-°- to 5 and 7; (x - x r ) to 6 



(r - m) 



(r - m + i) 



, Jump to P.P. 02 If * is out 
of range 



+ P.P. 02 



(May be loopstop or cue to 
extrapolation routine) 



(A + r + It) to 5 m . [Address of 



S* = ?LJlJx to 7 



f<*r) 



Sx to 0.0 



Overwritten by - 8* 



overwritten by Sx 



Sx to 0. 1 

Clear 7 J 

Table containing f(x r ) to Ul 

f(x r ) to 6 

Bring down next block of table if necessary 

(i - Sx) f(x r ) to 6 

/(* r+i > to 5 

(i - Sx) /(x r ) + Sx /<x r+i ) = /(*) to 6 

Plant LINK 

Optional P.P. 01 = (2.0, 0) 

Overwritten by LINK 

Optional P. P. 02 LOOP STOP if x out of range 
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POLYNOMIAL INTERPOLATION 

Method: See Section 6 of the Library Specification. 



CK> 



QE> 



CK> 



01 



02 



03 



R 





-0 3 


321 


- 


28 - 


0+ 


S 


72 


0.1 





60 


0+ 


@ 


72 


0.2+ 





60 


0+ 


13 


72 


0.4 





60 



LN 

POLYNOMIAL INTERPOLATION 



R 


0-0 1 


321 


- 02 - 


R 5 


2-0 2 


321 


- 02 - 



Call for P. P. 01 



Call for P. P. 02 



o.o 



0+ 



+0 



25) 7 52 
0.0 7 01 
0.3 60 



0.0 7 00 
0727 

1.0 7 00 7 
073 
[2] 72 7 
5+ [I] 72 
7 2 00 
6 7 00 
6 00 
6+[6l72 

2.2 7 03 



Entry 



+ P. P. 01 .= Address B.P of index 



Table number, q, to 7 ffl 

7 = B.P + a (Address of oth word in index) 



7 = B.P 
m 



Read address of table (A, n) from the Index 

Preserve accumulators 

Read first block of table to U2 

(A. «> 

x to 7 



X - X. 
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5+ 






6 23 




1.0 


2.3 6 24 




7 3 00 






* 


1.3 7 62 










2 

* > 


+0 




2.1 7 03 




3 


1.2 7 62 
1.5+ 6 60 




4 








1.2 7 60 






6 




^ > 


6 5 00 




2 6 00 




6 


@ 6 45 




2.0 6 04 




7 










0.1 60 














6+ 






40 0. 






0.0 















0.0 2 01 

0.2+ 6 63 






1 

' 








2.0 6 13 






2 


9 n r no 










6 1 00 




3 


(jT) 6 43 




<OS 


0.4 6 62 




4 


(T)6 41 
0.5 6 63 




S 


7 00 




6 


(T)0 55 
3 6 04 




l+®72 


to 1+.0 


1 


0.0 RO 













Check that x la within range 





R 


0-0 2 




321 


- 06 - 


01 


2 



- 00 0. 


02 


1.2 


60 






x - x to PQ 



x~x. 



to 3, x-x to 6 



Jump if r > 
+ P.P. 02 



_ / LOOP STOP or 

i cue to extrapolation routine 



r-m 

Jump if r N 
r - m + 1 
Jump if x - x r 
Jump if r = m - 1 
x-x to 5 
(A. n) 
n 
N-n 



Form n', where n' = number of points to be 
used for interpolation 



(A o . n) to 2 

Jump if N < n 

n to 2.0 

n' = min (n, N) to 6 

n' 

n' - 9 

LOOP STOP if n ' > 8 



LOOP STOP if n ' 4 1 



-(n -2) to 6 t if „ ' odd X7 * 
r - i(n'-2) to 6 



Optional Parameter List 

(Punched after Block 6+) 



Title 



Address B.P of index 



LOOP STOP in 1.2 if given x is out 
of range 



Issue 2 
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Form s 





f 


1+ 


from 6+. 7+ 


2+[l] 72 
0.3+ 7 60 




0. 

l 

2 

" 

3 




5 7 00 
2.3 7 03 




5 7 01 
0.3+ 7 62 




(7) 6 43 
1.0 6 63 






^ 


4 

S 




2. 1 1 13 
2.1 7 00 




6 7 03 
0.6+ 7 62 




6 




2.1 6 00 
6 3 03 




7 




35 6 20 
7 2 01 














2+ 


1 




2.3 4 00 
1 7 00 


1.0 




1 


072 2 
0.0 6 00 2 


■" 


2 




2.0 6 10 7 

1.3+ 2 66 




3 




072 2 
1.1+ 7 66 

1 4 20 








4 




6 




3+0 72 
35 2 00 


to 3+.0 


0.0 6 60 
(T) 4 51 

(T) 5 51 

i 44. n cn 




6 

7 











if Si < 
m-n' if S[ > m-n' 
s t otherwise 

r - |(n'-l) + S 

1 if n is odd and x-x r ^ - h 
otherwise. 



Jump if n ' is even 

x-x r 

x -x -h 

2(x-x r )-h 

Jump if x-x r > \h 

«! = r - i(n'-l) 

Jump if s < 0, leaving s = 0. 

(m-n') to 2.1 

(m-n') to 7 

(m-n 's ) to 7 

Jump if Sj ^ m-n' 

s = m-n ' 

r-s 



(A g , n) to 2 



Read block of table to U0 



Transfer f +i to 2.i, i = 0, 1.....7 



n'-l 
(n'-l)A 



2 ra = l 



Jump if (rc'-l)A is single length 
\(x-x r ) 
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Inner Loop: Neville' s Iteration Process 



from 2+.6 



o.o 



4+ 



1.0 



2.0 7 00 6 

2.0 7 10 6 
1.0 6 66 
0.0 2 11 
0.4+ 60 



6 



0, 



Sheet 4 of 4 



Clear shift marker 



h(r-s) to 7 


x-x to 5 


x-x to 0. 1 

s 


(n'-l-t) to 3 


x-x to 5 


(t + l)h to 0.2 


<l\ t0 6 


f(t) - f<*} 



0.0 is used to store 
shift marker 



0.1 is used to store 
x-x. 



0.2 is used to store 
(t + l)h 



S+l+l 



f(t+D 

s+i 

Jump if 0VR clear. 



Halve all the values of / in U2 



Add 1 to the shift marker 



/<* +1 > to 2.i 
s+t 

Add 1 to i 



^s+i + i 1 t0 5 



Count on i 
Count on t 
Restore accumulators 

jf' +1 > toe 

Shift mrker to 7 

Answer, /(x) , to 6 

Plant LINK and obey it in 2.0 

Ferranti Ltd., London Computer Centre, 21, Portland Place, LONDON. W. 1. 

Issue 2 Copyright Reserved D.M. 24th September, 1957 
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TWO WAY LINEAR INTERPOLATION 
Method 

See Section 5 of the Library Specification. 



Q2> 



Gl> 



from 1+.6+ 



01 



02 



03 



R 





-0 3 


322 


- 


28 - 


0+ 


m 


72 


0.1+ 





60 


0+ 


m 


72 


0.1 





60 


0+ 


m 


72 


0.3+ 





60 
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TWO WAY LINEAR INTERPOLATION 



R 


-0 1 


322 


- 02 - 


R 1 


2-02. 


322 


- 02 - 



Call for P. P. 01 



Call for P. P. 02 



0+ 




El 72 5^ 
1.0 5 00 51 

V o_|]j72 sy 

1+072 > 




6 7 03 1 



Sheet 1 of 3 



+ P. P. 01 



k to 5 m 
B.P + k 
Index to Ul 
A to 5 m 
Table to U2 

Plant LINK in 0.0 
x to 0. 1 
y to 0.2 
Put 6 in lp 
x or y to 7 



Overwritten by LINK 

Overwritten by x 
■ Overwritten by y 

Overwritten by my 
' Overwritten by f(x r y) 

Overwritten by (x-x r )/h x 

Overwritten by (y-y s )/hy 



x~x or y~y to 7 J 



■ Overwritten by loop marker 
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1+ 


1.0 


6 23 

2.0 6 24 1 




1 


7.5+ 7 10 1 
l.G 7 62 




2 




+0 




1.4 7 03 1 

1.2 '. 62 




3 




4 




©7 41 

1.5+ 6 60 








1.2 7 60 
2.0 6 26 1 




6 






tp 0+.6+ 8 


7.7+ 7 10 1 
0.6+ 1 66 




to 2+;2+ 7 


2+ [I] 72 
1.2+0 fif1 








". I 



(1) x-x to PQ (2) y-y to PQ 

(1) r to 7; x-x r to 6 (2) a to 7; y-y s 
(1) r to 0.3 (2) s to 0.4 

Jump if r, s > 

+ P.P. 02 > cue to extrapolation routine 



to 6 



(1) r-m,, 



(2) a~m^ 



(1) r-m x + 1 

Jump if * = x r or y - y s 



(2) s-m y + 1 



(1) (x-x r )/h x to 0.5 (2) (y~y s )/h y to 0.6 



Jump if first time round loop 







3+ 






0.7 10 
0.4 6 10 




1,0 


from 2+.7+ 


[|]72 5 
0.6 6 02 




1 

2 




2.0 6 20 5 
2.0 6 01 5 




1.4 5 66 
0072 5 




3 




2.0 1 00 5 

0.6 1 22 




4 

s 

6 




0.3 5 01 
35 5 03 




0.7 1 00 
1.0 1 61 






7 


6 5 00 
2+[2]72 



Indicate second time round loop 
f(x r y) to 0.4 

- (y-y s )/h y 

{i - (y-y s )Ay}f(*i y s ) 



fi*i y) 



r+i,s ■ m 



Jump if first time round loop 
/(* r +i y) t0 5 
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.7+ 








2+ 


from 3+ 




0.5 


6 02 






2.0 








0.4 


6 20 




0.4 


6 01 






1 








7+ 


2 


0.5 


5 22 




0.0 


60 


from 1+, 


0.3 


7 00 


to 1.; 


!+ 






34 


5 01 






1.3 


2.3 


7 20 




0.4 


7 01 






4 


35 


7 20 








7 


5 01 






S 


2.3 
35 


7 00 
7 20 






6 


0.3 


7 10 




3+ 


[I] 72 


tp 3+.1 




7 


1.1 













DU 



R 
322 


0-0 2 
- 06 - 


2 



- 00 0. 


1.2 


60 






- (x-x r )/h x to 6 

{l - (x-x r )/h^f(x r y) 

f(* y) to 6 

Jump to obey the LINK 

r to 7 



^oo t0 5 » 



nL r to 7 

m y r + s to 7 C 



m y r + s to 7 m 
A rs to 5 m 



"ly tO 7 C 



m y t0 7 m 



m y to 0.3 m 



Optional Parameter List 

(Punched after Block 3+) 



Title 



P.P. 01 



P. P. 02 



Sheet 3 of 3 



Perranti Ltd., London Computer Centre. 21, Portland Place, LONDON, W. 1. 
Issue 1 Copyright Reserved D. M. 31st July, 1957. 
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Sheet 1 of 5 



TVO WAY POLYNOMIAL INTERPOLATION 
Method 

See Section 6 of the Library Specification. 



01 



02 



03 



R 





-0 3 


323 


- 


28 - 


0+ 


m 


72 


0.2+ 





60 


0+ 


m 


72 


0.2 





60 


0+ 





72 


0.4+ 





60 



LN 

2 WAY POL. INT 



R 


0-0 2 


323 


- 02 - 


R 


1 -0 1 


323 


- 02 - 



Call for P.P. 02 



Call for P.P. 01 






Ql> 



0.0 



*dD 5 52 



0+ 



+0 



+0 



5 00 



0.1 5 01 
Q]72 5 

1.0 5 00 5 
073 



[3] 72 5 
1+072 



0.1 6 10 
5 3 00 



CUD 3 45 
3.0 3 03 



+ P. P. 02 Extrapolation Cue 



+ P. P. 01 Address of Index 



Table number k, to 5 B 

5 = B.P + k 
ft 

Read Index to Ul 
Address of Table to 5 
Write Accumulators to B0 
Read 1st Block of Table to U3 



m. 



x to 0. 1 

(A, n) to 3 

n to 3 
c 

n-N 
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1+ 




3 3 06 3 


1.0 






3.0 3 01 




3 1 00 


1 


@ 1 43 


c* 


1.2 1 62 


2 


(?) 1 41 


Q* 


1.3 1 63 


3 






1.2 2 00 




0.2 7 10 


4 






7.3+ 7 00 2 




6 00 


5 






2+ QO 72 




2.6 7 03 2 


e 






6 23 




3.0 6 24 2 


7 






7.5+ 7 10 2 



Sheet 2 of 5 



Porn n' = number of points to be used for 
interpolation 

n - N if n < N; if n Js N 

n r = min (fi, n) to 3 

n' to 1 

n' - 9 

STOP if n' > 8 

n' - 2 

STOP if n' < 2 

6 to 2n 



y to 0.2 
(1) x to 7 



(2) y to 7 



(1) w 



(2) y-y, 



(l) (*-x Q )/h x i2 ) (y-y )/h 



(1) r x to 0.3 



(2) r to 0.4 









2+ 


tp 0+i0 


0.0 7 63 
2.4 7 03 2 




2.0 


to 0+.0 


0.0 7 62 
(T) 7 41 




1 




2.3 6 60 
0.0 7 60 


tp o+.o 


2 








7.7+ 6 10 2 
7 00 




3 

4 
5 

6 

7 




1 6 00 
(7) 55 
6+[T]72 
7.5+ 6 04 2 




7.7+ 4 00 2 
7.3+ 10 2 


tp 6+.0 


7+[|]72 











Check that x is within range 

Extrapolate if r < 

r - m 

Extrapolate if r > m 

r - m + 1 

Jump if (x-x r ) or (y-^) = 

Extrapolate if r - m 



(1) *-* to 0.5 (2) y-y to 0.6 
r r 



n' - 2 

j (n'-2) to 6; to 7 if n' even 



r - } <n'-2) 
(1) *-* to 4 



(2) y~y r to 4 



Set 5 initially zero 



PEGASUS LIBRARY PROGRAMME 



Sheet 3 of 5 



Form s, where 



where 



from 2+. 7+ 



s 
s 
s 



if s < 

m - n'Mf s x > m - n' 
Sj otherwise. 

* L = r - |(n' - 1) + S 

1 if n is odd and x - x r i- ~ h; 



otherwise. 



l.o 



6+ 



1.2+7 60 
3.0 4 03 2 
7.7+ 4 01 2 
1.2+ 4 62 
(7) 6 43 
1.6+ 6 63 
3 7 00 
2.4 7 04 2 

6 7 03 
1.5+ 7 62 

7 6 01 
7.3+6 10 2 
7.5+ 6 13 2 
3.0 3 20 2 
3.0 4 00 2 
7.7+ 7 00 2 



Jump if n' even 
(1) x~x r -h x 
(1) 2(x~x r )-h x 



(2) y-y r -hy 
(2) 2(y-y r )-h y 



Jump if x~x r > i h x ; y-y r > ~ hy 



a' = r -i (n' - X) 



Jump if s' < 0, leaving s zero 



n' 



m - n' - $' 
Jump if m > n' + s' 
m - n' 
1 ) s x to 0. 1 



(2) s y to 0.2 



1) r x - s x to 0.3 (2) r v -$ v to 0.4 



y *y 



1) n'h x 

1) *-jc r to 7 



(2) n r h y 

(2) h y 

(2) y-y r to 7 







|, + 






2.1+ 6 60 
(T) 4 51 

7.7+ 4 10 2 




2.0 
1 




2 




7.5+ 4 22 2 
7.5+ 7 10 2 
2.7 2 66 
3.3 2 00 
0.1 2 20 
0.2 7 01 
35 7 20 
7 5 01 




3 

4 
5 
6 


to 4+.0 


4+ OH 72 
2.0 60 








1+[T]72 
1.4+ 60 




7 



Jump if n'h is single length 

1) \h x (2) |h y 

i) 1 (*-* r ) (2) i- <y-y r > 

1) \ to 0.5 (2) L to 0.6 

1) x^c r + (r-s)k (2) y-y r + (r-s)h 

1) x-x s to 0.3 (2) y-y s to 0.4 
Jump to form y-y s 



my to 2 



iHySjj 



+ s„ 



m y s x ' °y 



(A, 



s x"y 



- 1.0) to 5„ 



Return to B1+ to form s y , y-y s etc. 



Outer Loop 



Sheet 4 of 5 







3+ 


from S+.7 


4+1 
0.7 6 


72' 
00 


1.0 






1.6 3 67 
1.4 3 61 




l 




2 
3 




3.0 1 10 

(7j 72 

3.0 6 00 

1.4 1 10 




T 5+ [3] 72 A 

Lq_ 3 _ 4 _°J 

3.0 1 10 6 

2.6+ 60 




4 
5 






k 










1 K4. n 








6 
t 


[ 5+ 1 1 o 
0.1 5 00 




0.0 5 01 
2.2 6 66 












1 






4+ 


from 7+.6+ 




35 2 20 
0.0 7 10 






2.0 
1 




3 6 00 
35 3 01 




0.7 6 10 
0.1 5 10 




2 
3 




1 072 5 
1.0 7 00 5 






4 




3.0 7 10 2 

2.5+ 5 66 




5 




1 072 5 
35 2 01 








2.3+6 67 
0.7 2 00 






6 






tP 8+.0 7 


8+[T] 72 

1 n 9 err 













01 



02 



R 

323 


0-0 1 
- 06 - 


2 




- 00 0. 


0.0 


60 






Stage Modifier to 6 

Count n r 

Jump to Interpolate in x 

Answer to 3.0 

Restore Accumulators 

Answer to 6 

Plant LINK 

> Overwritten by LINK 
(0, 1) to 3 

f< x s+n *> t0 3 - n 

f(* s+ i yy to 5+. i 

Add m y to form address of f(x s+i+1 y g ) 

Increase Stage Modifier 



*»,, to 0. 0„ 
y m 



n r to 6. 



(0.1, n') to 3 

(Stage Modifier, n') to 0.7 

Table Address to 0. 1 



Transfer n' consecutive values from Table 
to U3 



n r to 2, 



n' ~ 1 to 2, 



Optional Parameter List 

(Punched after Block 4+) 

Title 



(2.0, 0) Address of the index 



LOOP STOP 
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Inner loop: Neville' s Iteration Process 



Sheet 5 of 5 



1.0 



8+ 



,127j 2 45 
4 00 



0.2 10 
9+f2l72 



2 1 00 
0.3 5 00 3 



0.5 4 01 3 
3.0 6 02 1 



3.1 6 01 1 
5 6 20 



4 6 25 
3.0 7 01 1 



2.2+ 64 
6 00 



3.0 7 00 6 
(T)7 51 



2.0 



7 
<3 



3.0 7 10 6 

1.7 6 66 

35 6 00 

0.2 6 11 

1.3+ 60 



3.0 7 10 1 

35 1 01 

0.5 5 03 3 

1.3+ 1 67 

1.2 2 67 

7 1 00 

0.2 7 00 

3+|T|72 

(V) 1 50 7 



1.0 64 

2.7+ 60 



(0, n> - 1) to 2 



Make Shift Marker zero 



n' - 1 - t 



y - y or*-* to 5 



(t + l)h 

-'♦■ (t> 
s+i 

(t) (t) 

Z — 2, 

s+i-n s+i 

(v l+1 (tl -* s +i (t)) <* - w/<* + 1)h 

s+i 
Jump unl ess OVR set 



Halve all values of z in U3 



Add 1 to Shift Marker 



(t+D 



s+i 



*-y S H+i ° 



r x - x 



s+t+i 



Ferranti Ltd., London Computer 
^site 1 Copyright Reserved 



Count on i 
Count on t 
Interpolate z to 1 
Shift Marker to 7 

Shift up z if necessary 

Stop if | z | > 1.0 

Centre, 21, Portland Place, LONDON, W. 1. 

D.M. 19th August, 1957. 
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This subroutine estimates the value of y - f(x), for given x, by interpolating 
in a table of values of f(x) given at values of x which are not necessarily 
equidistant. 

The method used is Neville's iterative process which is described in Section 6 
of the Pegasus Library Specification of R 327. 



01 



02 



R 


0-0 2 


327 


- 28 - 


0+ 


©72 


0.1 


60 


0+ 


[0J 72 


0.3 


60 



LN 

POLYNOMIAL INTERPOLATION (UNEQUAL INTERVALS) MK 3 



Call for P. P. 01 



R 


0-0 1 


327 


- 02 - 


R 4 


7-0 2 


327 


- 02 - 



Call for P. P. 02 











0+ 




0.0 




+0 


Cm") 


^_ 


35 


7 22 


< Hi I 






1 










0.0 


7 01 









Q] 72 7 




2 










1.0 


7 00 7 


Cr)""} 







[7] 73 


\. U" } 






3 













[3] 72 7 






3.0 


5 00 




4 










® 


4 40 




35 


4 22 




5 










7 


2 00 




35 


5 22 




6 










1+ 


CD 72 




7+ 


072 


tp 7+.0 


7 


0.0 




(!sheet 4) 






ou 



from 7+.0+ 



or 7+,2 
(sheet 4) 



PERRANTI LTD 
PEGASUS LIBRARY PROGRAMME 

Entry 



+ P. P. 01 = (B.P, Q) 



R 327 



i.o 



1+ 
/ 3+ [2j72 \ 



I 3.1 1 00 4 J 
TflTh 7 45 \ 



5 



;t27) 7 45 

I 1 7 04 

4+ [3] 72 

1.3+ 7 63 

7 1 03 

1.0 1 10 

1.1 5 10 
(9) 1 43 

1.5 1 62 

_0_ I « 

1.6 1 63 

® 1 53 

<2?) 1 52 



Sheet 2 of 5 



q to 7 m 

Address of <?th entry in index to 7 m 

7 m = A q (address of gth table); 

7 C = "g (no. of points for interpolation) 

Store Accumulators 

1st block of table to U3 

m to 5 C (m = no. of values tabulated) 



A q +3 



2m = A q + 3 



= **n> 



7 m = 4- + m + 3 = X(y.) 



Form n r , where n' = number of points to be 
used for interpolation 



Nto 1 £ 
n to 7, 



(N ~ n) to 7 C 

Jump if N < n 

Replace N in l c by n 

Store n ' in 1.0 

Store ik 1.1 

«' -9 

LOOP STOP if n' ^ 9 

n' - 2 

LOOP STOP if n' < 2 

n' + 1 to l c 

*m ~ 2 " r (say) 



Overwritten by n' 



Overwritten by m 
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3+ 



2.0 



6 51 

|0]72 2 

2 4 00 

0.0 7 00 4 

7 51 

6 7 04 

2.5+ 7 63 

2.4+ 4 66 

[0]72 4 
2.1+ 5 67 

3.7 7 61 

2 4 03 

3.7 4 60 

1 4 03 
3.0 4 62 

4 00 



Sheet 3 of 5 



Check that x is within range 

x/2 to 6 

Block containing x Q to U0 

4m = Mx ) 

Jtj to 7 

( Xi /2) to 7 

(x/2) - (x { /2) to 7 

Jump if * < *:j 

Step modifier in 4 

Read new block of x where necessary 

Count on m 

Jump if x > * m . t (out of range) 

i = ,4(*) - A(x Q ) to 4 

Jump if x < x Q (out of range) 

i - r to 4 

Jump if t ^ r 

Set s = if i < r 



from 2+. 4 + 
(sheet 4! 



to 2+.0 



* 



sheet 4' 



3.0 



4+ 



1.0 5 03 
(0) 5 41 1 
3.3 5 62 

1.1 4 00 
1.0 4 03 

!5) 4 52 

4 2 01 

4 3 01 

1.0 1 00 

2+ [T]72 

35 4 02 

1 5 00 

Qj]72 3 

1.0 60 



^ 



I 



+0 



Form s 

m - i - n to 5 C 
m-t-n+rto5 c 
Jump if i - r 4 m - n 



Set s=m-n if i - r > m - n 

*m = ^<y s > 

«' to 1 

Set marker for transfer; ^ - -0. 1 

n' to 5 

Block containing y s (or x s ) to U0 



+ P.P.02 



Loop stop or cue to extrapolation 
routine if out of range 
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Sheet 4 of 5 



from 4+. 6+ 



(sheet 3) 



& 



4+. 5+ 



< sheet 3 i 



l.o 



to 5+.0 



(sheet 5) 



2+ 



0.0 7 00 3 

1.5 4 60 

2.0 7 10 5 

1.2+ 3 66 

15172 3 

35 5 01 

1.0 5 67 

1.7 4 66 

2 3 00 

3.5+ 60 



(7) 7 51 
6 7 03 
3.0 7 10 5 

1. 1+ 60 



5+ [0]72 
0.0 1 67 



x s+i (>s+i> t0 7 
Jump if x 

y s+i to 2.i 

Step modifier in 3 

Next- block containing y s (or x s ) to U0 

Step 5 m 

Count on n' 

Jump and set shift unit second time 

A(x s ) to 3 

Jump to read y s 

l/2(x s+i - x) to 3. i 



Jump and set (n' - 1) in 1 











7+ 


from 0+.7+ 




7 


3 00 


(sheet 2) 


0.0 






1.0 
2 

7 


1 63 
3 00 

2 00 


to 1+.0 


1 


(sheet 2) 




1.0 


4 66 




2 











R 


0-0 2 






327 


- 06 - 






2 


- 00 0. 




0.3 








0.4 









3.7 


60 








\ = A q + m + 3 = ^(y ) 

Jump if not inverse interpolation 



Interchange addresses of x and y for 



inverse interpolation 



4^ = 0.1 for inverse interpolation 



Title of Optional Parameter List 



Optional P. P. 01 



Optional P. P. 02 



Address of Index 
B.P = 2.0 



LOOP STOP if 
out of range 



FERRANTI LTD 

PEGASUS LIBRARY PROGRAMME 

Inner Loop: Neville' s Iteration Process 



R 327 



Sheet 5 of 5 



from 2+.7+ 



o.o 



10 



2.0 6 10 5 
0.6 5 66 



1.0 



6+ 



0.0 4 11 
0.2 60 



Set initial shift marker 
in 0.0 



Set 3^, = 



Overwritten by 
shift marker 



Jump if 0VR clear 



Set 5, 



m 



Divide each value of y hy 2 



Add 1 to shift marker in 0.0^, 

y& 1} to2.i 

Step y modifier 

Step x modifier 

Count on i 

Step t 

Count on t 

Restore Accumulators 

y?'» to 6 

Shift marker to 7 

Scale up y^ n "^ if necessary 

Jump if LINK is a go order-pair 

Plant LINK negatively otherwise 

Overwritten by LINK 
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DETERMINATION OP THE ZERO OP A FUNCTION 
Method 

See Section 4 of the Library Specification. 



01 



R 





-0 1 


340 


- 


28 - 


3+ 





72 


0.0 





60 



LN 

ZERO OF F(X) 



R 4 


3-0 2 


340 


- 06 - 


R 4 


1 -0 1 


340 


- 02 - 


R 3 


3-0 2 


340 


- 02 - 



Title of Optional Parameter List 



Call for P. P. 01 



Call for P. P. 02 









3+ 






4+QJ72 
1 073 


/■ nl 1 


^ 01 J 


0.0 

1 

2 




1.4 6 10 

1.5 7 10 




0.7 1 00 

4+{T|73 




3 


^to Auxiliary 


+0 






from Auxiliary 


4+ [7] 72 
1.6 6 10 


(1st time) 


4 
5 




1.2 1 00 
1.5 6 00 


from 2+ 5+ 


0.2+ 60 








e 

7 









3+072 
0.4 60 



Preserve Accumulators in Bl 



x a to 1.4 



»5 to 1.5 

Set LINK 1 for Auxiliary 

Preserve Ul 

+ P. P. 02 = cue to Auxiliary 

Restore Ul 



F a to 1.6 
Set LINK 2 



Xfr to 6 

Enter Auxiliary 



LINK 1 for Auxiliary 
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from 1+.7 



to. 



2+.1 



l.o 



4+ 



2+[o)72 
0.1 60 



+0 

2+072 
0.0 60 




Start next cycle 



+ P. P. 01 * c 



LINK 2 for Auxiliary 



Optional P. P. 01 = e 

Optional P. P. 0.2 
= cue to Auxiliary 



Overwritten by x 



*6 











2+ 


Auxiliary 


4+ 
1.7 


072 
6 10 


(2nd Time) 


0.0 


Next cycle 


1.7 
1.4 


6 00 

6 20 


from 4+.0+ 


1 

2 
3 

4 
5 




1.6 
1.5 


5 02 
5 22 




1.7 
5 


5 01 

6 25 




1.3 
7 


7 10 
6 00 


tp 3+.6 


0.6 

3+ 


1 00 

072 




e 

7 




1+ 
0.0 


[0]72 
60 




+0 



Restore Ul 
F b to 1.7 
F b to 6 

~F a to5 
F b- F a 



*c = l*fb ~ Va>M ~ F a> t0 7 



x c to 1.3 and 6 



Set LINK 3 for Auxiliary 
Enter Auxiliary 

LINK 3 for Auxiliary 
Not used 



PE 
PEGASUS 



RRANTI LTD 
LIBRARY PROGRAMME 



R 340 



Auxiliary 

(3rd time) 




to 4+.0 



o.o 



0+ 



+0 

1.7 7 20 
1.6 5 11 
1.6 6 25 

0.4+ 60 



0.1 6 62 

1.4 6 00 

1.5 6 10 
1.7 7 10 



1.3 6 00 

1.4 6 10 



1.6 5 10 
1.5 6 03 



1+072 
0.5 60 



Sheet 3 of 3 



Restore Ul 


Exit if F c 


F c to 5 


F a to1 



= 



C(6) < if F a and F c have opposite signs 



Restore Accumulators 



x c to 6 



Plant and obey Main LINK 



Next cycle if OVR 



Form € - \x c - x b \ 



Next cycle if -\x c - x^\ > e 
Otherwise EXIT. 



EXIT 



, Overwritten by 
Main LINK 



Not used 



F a +F C to 1.6 



F + F 



► F Q and F have same sign 



F> h = fo 
a (or 

x c to 6 



a c 



F' = F„ 



x c ~ x b 



. F and F c have opposite 
signs 



as above 
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LINEAR INVERSE INTERPOLATION 

This subroutine estimates the value of the argument x corresponding to a given 
value y = f(x) of a function f(x), tabulated at equal intervals of x. 

The subroutine searches through the table, starting at y = f(x ), until it finds 
two adjacent entries such that 

y r < y « y r+1 
or y r > y ^ y r+1 

If there are two or more values of r for which - one of these conditions is true, 
the smaller value of r will be selected. The required value of x is given by 



>r+i - y 



X « Xf j. i 



yr+i - y t 



where h is the tabular interval (xj +1 - x^). 



01 



02 



03 



R 





-0 3 


341 


- 


28 - 


0+ 





72 


0.1+ 





60 


0+ 


m 


72 


0.1 





60 


0+ 


El 


72 


0.3+ 





60 



LN 

LIN. INVERSE INT. 



R 


-0 1 


341 


- 02 - 


R 3 


6-0 2 


341 


- 02 - 


R 1 


6-0 2 


341 


- 06 - 



Call for P. P. 01 



Call for P. P. 02 



Title of Optional Parameter List 
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o.o 



(02> 



(wy 



0+ 



+0 



7 00 
{25) 7 52 

0.0 7 01 

\T\l2 7 

1.0 7 00 7 

0J73 

072 7 

4 00 

1.4 2 00 

6 2 03 

1.3 3 00 

35 7 01 

1+ ^72 

0. 1 7 66 



from 2+. 7+ 



l.o 



0. 1 



"v. 2. 



to 3+. 



1 + 



0.0 6 00 
1.1 1 10 



0.5 7 05 
0.2+ 2 63 



16} 4 46 
0.2 65 



1.2 2 00 
1.1 7 01 



3+ [5]72 
0.0 7 66 



7.7 7 70 0, 


2 - 00 0, 




0.6 60 
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ENTRY 

+P. P. 01 = Address B. P of index 

Clear 7 

Table number, q, to 7 ffl 

(B.P + q) - address of q word In index to 7 m 

Block containing index to Ul 

04, n) to 7 

Store the Accumulators 

1st block of table to Ul 

Clear 4 

y a = f(x Q ) to 2 



y to 2 



h to 3 



7 m = A + 1 



7 m = A + 2. 



Unconditional jump 



INITIAL SETTINGS 

x to 6 

Plant LINK and Exit 

Clear 7 C 

If y < y leave in 4 

If y £ y set 2" 3 '* in 4 

Correct 4 if y - y has overflowed 

:k to 2 

m to 7 C 

l m = A +3. Unconditional jump 



Collating Mask 



Optional P. P. 01; Index address =2.0 



Optional P. P. 02; £°? *£ P J f t{£?> ° Ut 
of range of table 







H 




^ 


35 7 03 




s * 




0.0 


7 1 04 






1.0 5 02 7 






I 


0.2+ 7 66 






' 








4 5 02 






2 


1.0 5 01 7 






? 






1.0 3 20 7 




3 


5 6 25 




0.0 7 12 




4 


1-f 072 




1.5 1 05 




5 


1 3 20 






35 2 22 




6 


35 6 25 




0.0 7 11 


,to 1+.0 7 


n ("vl in 


\~ 






U |_ 





PERRANTI LTD 
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INTERPOLATE 



R 341 



from 1+. 4+ 



o.o 



3+ 



7 1 00 
0.3 4 16 
0.3 7 66 
1.7 4 00 

072 7 
-0 



>| 1.0 6 13 7 

1.0 5 00 7 

0.5 65 

0.7 5 62 

0.1 1 67 

35 1 01 



+0 



2+ 072 
0.0 60 



Sheet 3 of 3 

A + 4 + r to 7 ffl [-=ilOc ) + r] 

r + 1 to l m 

y - y r to 5 unless y r at end of block 

Jump unless at end of block, A(x r+ .) to 7. 

(y - y r ) to 5 

(y r+ i - y r ) to 5 

(y r+1 - y).A to fl2 

_ ov+1 - y) 



(yr+i - >v) 



•h to 0.0 



Clear l c ; l ffl = {r + 1) 



(r + 


1)A,2 


" 13 to 


PQ 


r r+i* 


2 -13 


to PQ 




*r+i 


to 7 






x to 


0.0 






Restore Accumulators 



© 



SEARCH FOR y r 

(i + 3, m) to 1 

Change 0.3+ to 1.0 5 02 7 If y > y 

Unit modify. 7 m =,4 + 3 + 1+1 

y{ - y to 4 at end of block 

Next block of table to Ul 



^t+i y I if y > y b-order Is 

±w +1 -y)J 1 -° 5027 

Jump if 0VR set 

Jump when y lies between y r and y r+1 

Count entries in table (m) 

Form U + 0,4, 0) in 1 if f(x) outside range 

of table 

+P P 02 ] May be lo0p stop 0r cue 
J to extrapolation routine 

y r found 

Enter interpolation sequence 
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FLOATING-POINT ARITHMETIC 

The subroutine Is self-preserving and carries out the operations of addition, 
subtraction, multiplication and division, on the packed floating-point numbers 
x = A.2 a and y - B.I , leaving the result in X6. 





R 0-0 2 




610 - 28 - 




0+ 00 0. 


01 














02 






0+ 00 0. 



ln 

floating-point arithmetic 



R 


5 


-0 1 


610 


- 


02 - 


R 


6 


-0 1 


610 


- 


02 - 


R 


6 


-0 1 


610 


- 


02 - 


R 1 


7 


-0 2 


610 


- 


02 - 


R 2 





-0 3 


610 


- 


02 - 



Calls for P. P. 01 



Call for P. P. 02 



Call for P. P. 03 
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On entry 1 = for x and - 
= 1 for + and + 







0+ 






32 1 01 

6 4 00 

3.6 6 05 

6 4 03 


C i ^ 


^ ' 


3.0 


r ■■ j. \ 


{ X + J— 


1 




2 


7 5 00 
3.6 7 05 




3 


7 5 03 
4.2+ 1 63 1 










4 


3.5 5 01 
3.7 60 1 






5 
8 




+0 




+0 

7 6 20 
4.2 60 






f™ 


7 














1+ 




* 


7 7 01 
7 6 26 
7 6 00 

(T)5 44 
5 4 01 

5.1+ 1 62 
7 7 02 
4 5 04 






4.0 
1 


tn 2+ 1+ 




2 








+ 3 










4.7 5 62 
5 4 03 




4 
5 
6 




7 6 06 
6 7 06 




7 6 06 
5 5 02 

^0^5 43 0. 

5.0+ "> 62 


to 2+.0+ 


7 












1 



Make 1 negative to indicate + or - 



A to 6 



a + 2*" 1 to 4 



B to 7 

b + 2 71 " 1 to 5 
Jump if + or - 
b to 5 



• UNPACK OPERANDS 



x or * 



+ P.P.01 = -2 n - 1 .2" 38 



+ (P.P.01) x 2 = -2 n .2* 3S 



AB to PQ 



- + 



25 to 7 

A/2B to 6 

Set 0VR if B - 

1 - b to 5 

a+b+2 n ~ 1 If x; a-i>+l+2 n_1 if + 

Jump if x or + 

Change sign of B 

a-b to 5 

Jump if a ^ b; in this case d = a-b > O 



b + 2 n_1 to 4 



Interchange A and B \ If a < b 

d = b-a > to 5 

+ P. P. 02; rf-(37-w) to 5 

Jump if d > (37-n) 



>+h- 
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* k = min (d, 37-n) ^ 









2+ 




(25) 5 50 


S. 

from 1+.7+ 


(V)7 51 5. 
7 6 01 




from 1+. 2+ 


5.7 6 60 
(j36) 4 56 




2 






3.5 6 03 




3 


(jT)4 56 




4 


5.6+ 4 63 

3.6 6 05 

4 6 01 




3.6 4 05 






5.7 4 60 




32 02 




e 






»• 

' *■ 

i 


6 00 




(25) 1 52 






60 1 
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, NORMALIZE AND 
PACK RESULT 



d-(37-n) to 5 m 

+ P. P. 03 Shift argument down k* places 

Add arguments 

Jump if argument = 

Normalize 

Round argument 

Re-normalize 

Jump if underflow 

Pack result in 6 

Jump if no overflow 

Set OVR if overflow 
Set result = 

EXIT 



■ + &- 







3+ 




R 0-0 3 




610 - 06 - 


01 


-256 




28 00 


02 














03 






28 00 



Optional Parameter- list 



Title 



= -2 n - 1 .2' 38 



(37-n) 



(37-n) 



, n = 9 
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FLOATING-POINT SQUARE ROOT 

Method 

This subroutine evaluates F'(6) = /f(6) where F(6) = x = A.2 a Is the standard 
floating-point number held in X6. 

If we put A' = A and a' = a for a even 

and A' = 2A and a' - a-1 for a odd, 

then Vr = /J 7 ^^'. 

An iterative Newton -Raphson process is used to find y = i/A f . 

y n+ i = y n + d n where d n = iif n -y n ) 

= W +■ i 
o ir 2 



with 



i > 



The process terminates when d n ^ 0. y n . 2 2 is then taken to be the answer and is 
left in standard floating-point form in X6. There is no need -to test for floating- 
point overflow or underflow since the answer is always nearer to one than the operand. 



01 



R 





-0 1 


611 


- 


28 - 


0+ 


m 


72 


0.0 





60 



LN 

P.P. SQ. ROOT 



R 1 


5 -0 1 


611 


- 02 - 


R 1 


6-0 1 


611 


- 02 - 


R 1 


6-0 1 


611 


- 02 - 



Calls for P. P. 01 
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0+ 


f m "N ^ 


6 5 00 
6 7 00 


0.0 


1 


(?) 7 45 




l+[l]72 


2 


0.2 6 63 
1.6 6 05 




3 

4 


6 5 03 
0.5 7 60 




(7) 5 43 




>■ 


6 6 01 




1.5 5 03 


e 


(7) 5 51 




6 


1.5 1 10 
1.5 6 60 




7 




6 1 00 
6 7 00 














1+ 




1.0 


©7 51 








1.7 7 11 






1 


1.7 1 26 
1.7 7 03 






2 


~0~7~5~f~ 








1.0+7 63 






1.7 6 00 




3 


(76) 5 56 




4 


1.6 6 05 

5 6 01 




+0 

v_ 1 




5 


e 


+0 


1 


♦0, 



2+ 



R 0-0 1 
611 - 06 - 



-256 
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x to 5 
x to 7 
If a even, to 7; if o odd 2* 3B to 7 

LOOP STOP if x < 

A to 6 

a + 2 (B " 1) . to 5 

Jump if a even 

(a - 1) ■+ 2 (n_1 > to 5 

24 to 6 

a' + 2.2*""^ to 5 (a r = a or a - I) 

(a'/2) + 2 (n_1 > to 5 

Plant LINK 

EXIT if A = 

4 r to 1 . (A* = A or 24) 

l r to 7 



?4' to 7 

Initially y ={i4'+l) to 1.7; then y n+1 =y n +d„ 



2d rt to 7 



(s-*) ; .- 

d^, to 7 
Jump if c£„ < 
vT' to 6 
Normalize 



Pack in P.P. form: /x to. 6 



+ P.P. 01 



j(rt-l) j2 -38 1 



Overwritten by LINK 



+ 2 x P. P. 01 



-2 M .2- 3e 



Overwritten by y n ; finally = /I* 



Optional Parameter List 



P.P.01 = -2"" 1 (n = 9) 
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SHORTER FLOATING POINT ARITHMETIC 



This is a self-preserving subroutine which performs the operations of addition, 
subtraction and multiplication on packed floating point numbers and can also be adapted 
to do division instead of multiplication. The floating point number x = A. 2 a is held 
in a single word with the least significant n bits representing the non-negative integer 
a + 2 n ~ 1 , and the most significant 39-w bits representing the fraction A. 



01 



02 



03 



04 



R 





-0 


4 


612 


- 


28 


- 


0+ 


Ill 


72 




hE 


72 
70 


— 


|2 + 


5| 


4.4 


1 


10 




0+ 





00 


0. 




















0+ 





00 


0. 



LN 

SHORTER P.P. ARITHMETIC 



R 1 


6 -0 1 


612 


- 02 - 


R 1 


7-0 1 


612 


- 02 - 


R I 


7 -0 1 


612 


- 02 - 


R 2 


7 -0 1 


612 


- 06 - 



Calls for P.P. 01 



Title of Optional Parameter List 
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0+ 






32 1 01 

6 4 00 

5.7 6 05 

6 4 03 


t , "\ 


( \ ■■ ) w 

4.0 
( r \ > 


\ * - ) * 

1 


2 


7 5 00 
5.7 7 05 




3 
4 


7 5 03 
4.6 1 63 


4 to 2+.0 


/5.6 4 01 \ 

I 6 720 J 






^ from 2+.4 


5 4 01 
5.3 60 






5 






{ , 




7 7 12 1 
4 5 04 
(25) 5 50 
5.2 5 62 








6 

7 




















1+ 






5.0 
1 


5 5 02 
4 41 5 
(j)) 6 51 5 
5.2+ 60 










(jj) 7 51 5 
7 6 01 








2 












(40) 4 56 
5.6 6 03 




3 


4 


5.7 6 05 4 
4 6 01 4 


5 


(25) 1 52 
64 1 


6 


+0 


7 


+0 



Make 1 negative to indicate + or - 

A to 6 

a + 2"* 1 to 4 

> UNPACK OPERANDS 

B to 7 

b + 2"" 1 to 5 J 

Jump if + or - 

Overwritten by orders jumping to 2+. 
. „ . „ when routine adapted for division 

a + b + 2"" 1 if x; a - 6 + 2 + 2""Hf-f 

Negate 7 if - 
a - 6 to 5 
- b to 5 m 
Jump it a ^ b 



b - a to 5 m 

a + 2"" 1 + b - a = b + 2 

Shift i down b - a places 



n-i 



to 4, 



Shift 5 down a - b places 

Add Arguments 

Normalize 

Round argument 

Pack result In 6 if exponent > - 2"" 1 

Clear 6 if underflow (exponent < - 2 n ~ l ) 

Exit provided OVR clear 

(otherwise 77 stop in 5.6) 



+ P. P. 01 



. 2 n-i_ 2 -3 8 



77 stop here if 
OVR is set 



+ (P.P. 01) x 2 



-2 n .2" 3e 
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2+ 


from 0+.4+ / — *. 


4.0 7 60 


4.0 






5.6 5 01 




© 5 44 


1 


6 51 




7 6 26 


2 






7 6 00 




4.5 7 00 1 


3 











0+[4)72 


to 0+.5 4 


4.4 7 10 






2+ []] 72 


5 






4.0 60 




5.6 4 01 


6 






6 7 20 


7 


-256 



Loop stop If zero divisor 

b to 5 

2 - b to 5 

\ A to 6 

y 4 a/b to 7 

H Vs to 6 

Order pair to reset for -r or x to 7 

Plant reset order pair in 0+. 4 
and then return to 0+.5 

-i- reset order pair 
x reset order pair 
Optional P.P.01 = -2 n_I .2~ 38 (n - 9) 



Notes: 



1) The value of the parameter n, specifying the number of bits in the 
exponent, may not exceed 10 because the order 4 41 5 in 1+. 0+ 
assumes that there are not more than 10 bits in b - a. The maximum 
possible value of b - a is 2 n - 1. 

2) The order 5.6 6 03 in 1+. 3+ may very occasionally make the argument 
equal to 4% or ~ l A, which are not correctly normalized. See section 4 
of the Library Specification. 
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